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Th«  picture  which  emerges  from  our  investigations  is  that  vs  are  Indeed 
observing  the  behavior  of  a  two  dimensional  magnetic  system.  The  low 
damans ionellty  is  brought  about  by  the  shielding  of  the  intercalent  layers  by 
the  Intervening  graphite  layers.  The  shielding  also  depends  on  the  charge 
transfer  between  the  lntercalant  and  the  graphite.  In  graphite  there  Is  a  band 
near  K  ■  0  about  half  filled  with  electrons  (10).  The  intercalation  of 
an  acceptor  compound  will  draw  off  those  electrons  so  that  the  electrical 
conduction  Is  hole  .like.  Thus  low  stages  exhibit  a  positive  Hall  coefficient. 

The  shielding  of  the  lntercalant  layers  is  due  to  electromagnetic  shielding 
analogous  to  the  skin  effect,  irtiere  the  skin  depth  is  a  function  of  the  electrical 
conductivity*  Which  in  turn  is  a  function  of  the  number  of  charge  carriers. 

Stages  1  M  )  Almost  empty  the  band  of  electrons  and  thus  the  graphite  layers 
have  maximum  conductivity.  Above  stage  3,  the  quantity  of  acceptors  is  not 
large  enough  to  ebsorb  all  the  electrons  which  can  be  donated  by  the  graphite, 
and  the  conductivity  of  the  graphite  layers  decreases.'  Moreover,  the  conduct¬ 
ivity,  according  to  our  Hall  measurements,  is  due  to  both  positive  and  negative 
carrlers.«rthia  gives  rise  to  decreased  magnetic  shielding  of  the  lntercalant 
layers,  vluka  corresponding  decrease  in  the  peak  height  in  the  susceptibility 
at  1.7K.  TVi\,  the  transition  at  1.7K  can  become  a  very  sensitive  probe  of  the 
internal  merauc  field,  the  charge  transfer,  and  the  electrical  conductivity 
of  the  graphics \ayers. 
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Review  of  Accomplished  Work  and  Work  in  Progreso 

I.  Introduction 

Moet  of  the  work  to  date  was  performed  on  FeCls  intercalated  graphite,  where  the 
FeCls  is  the  magnetic  substance.  Work  on  FeCls  intercalated  graphite  at  Boston  University 
was  begun  in  the  late  1970’s  and  was  initially  confined  to  Mossbauer  investigation  of  the 
compound  [lj.  In  those  investigations  it  was  shown  that  the  FeCls  had  a  magnetic  phase 
transition  with  long  range  order  in  stage  1  and  stage  2  compounds,  while  no  such  transition 
took  place  in  stages  4  and  6.  Subsequent  measurements  of  the  magnetic  susceptibility, 
which  are  described  in  more  detail  in  Appendix  I,  have  shown  that  besides  the  susceptibility 
which  denote  a  ferromagnetic  transition  in  the  out  of  plane  direction  and  an 
antiferromagnetic  transition  in  the  in-plane  direction  in  stage  1  and  an  antiferromagnetic 
transition  both  in  the  in-plane  and  out-of-plane  directions  in  stage  2  [2,3]  There  exists 
an  anomaly  in  the  magnetic  susceptibility  which  occurs  in  all  stages  at  about  the  same 
temperature,  1.7K  to  1.8K. 

Graphite  intercalation  compounds  (GIC)  are  naturally  layered  [4]  and  are  of  interest 
because  it  was  believed  that  electronically,  only  the  graphite  layers  next  to  the  intercalant 
change  significantly  upon  intercalation  [5].  From  our  measurements  we  now  have  evidence 
that  the  charge  transfer  affects  layers  other  than  those  next  to  the  intercalant.  The  in¬ 
tercalant  forms  a  monomolecular  layer  which  is  boarded  by  graphite.  Fig.l  of  appendix 
II  shows  a  schematic  of  stage  2  intercalated  graphite  compound,  where  the  striped  layers 
represent  the  intercalant,  (stage  index  refers  to  the  number  of  graphite  layers  between 
two  consecutive  intercalant  layers).  For  a  magnetic  intercalant,  the  graphite  layers  can 
be  a  controlling  factor  in  the  magnetic  interactions  between  different  magnetic  intercalant 
layers.  The  graphite  layers  can  control  the  magnetic  interaction  by  separating  successive 
intercalant  layers  from  eacn  other.  Increasing  the  distance  will  generally  decrease  the  ex- 


change  and  dipolar  interactions.  They  can  also  actively  shield  the  intercalant  layers  from 
each  other  because  the  graphite  planes  are  good  electrical  conductors.  Thus  there  will 
be  electromagnetic  shielding.  The  systematic  change  of  the  number  of  graphite  layers 
between  two  adjacent  intercalant  magnetic  layers  can  provide  two-dii&ensional  magnetic 
systems  by  shielding  the  intercalant  layers  from  each  other  until  their  magnetic  interaction 
is  insignificant.  In  our  investigations  we  found  that  in  low  stages  this  shielding  increased 
exponentially  with  stage.  Thus  intercalation  of  magnetic  species  into  graphite  provides  us 
with  very  useful  model  systems  in  which  to  study  two  dimensional  mugnotium 

Magnetic  compounds  have  been  widely  used  in  the  intercalation  of  highly  ordered 
pyrolytic  graphite  (HOPG)  and  generally  one  observes  magnetic  anomalies  in  the  GIC’s, 
even  in  high  stages.  The  structure  of  the  intercalant  is  similar  to  that  of  the  bulk  material 
and  generally  magnetic  compounds  intercalated  into  graphite  also  undergo  a  phase  transi¬ 
tion  in  the  intercalated  state,  but  the  nature  of  the  transition  is  altered  by  the  enviroment. 
Additional  transitions  due  to  the  magnetic  two  dimensionality  of  the  systems  were  also 
reported  [6,7]. 

We  have  performed  a  systematic  study  of  the  magnetic  phase  transitions  of  FeCl3 
intercalated  graphite  compounds  and  their  behavior  as  a  function  of  the  stage  of  the 
sample,  (see  appendices).  We  also  investigated  the  dependence  of  these  transitions  on  an 
external  dc  magnetic  field. 

The  crystal  structure  of  FeCl3  is  a  repeated  sequence  of  chlorine,  iron,  and  chlorine 
hexagonally  arranged  atoms  [8],  such  that  every  iron  atom  is  surrounded  by  an  octahedron 
of  chlorine  atoms.  The  iron  atoms  in  subsequent  layers  are  displaced  with  respect  to  each 
other  by  1/3  of  the  unit  cell  and  thus  three  FeCl3  layers  are  encompassed  by  the  cell  vector 
in  the  C  direction.  Upon  intercalation  of  FeCls  into  the  graphite,  one  electron  is  donated 
by  the  graphite  host  for  every  four  iron  atoms  (9).  It  is  not  yet  clear  what  site  these  donated 
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ciectrona  occupy.  At  room  temperatures,  Mossbauer  data  [10,11]  and  Raman  spectroscopy 
[12]  do  not  support  the  existence  of  FeCla  as  a  result  of  the  acceptance  by  FeCl3  of  the 
donated  electron  [13].  However,  in  the  low  temperature  Mossbauer  data  for  FeCU,  Millman 
and  Kirczenow  [14],  at  Boston  University,  have  seen  Fe2+  ions  at  temperatures  as  high  as 
100K. 

II.  Magnetic  Susceptibility  (Appendix  I,  II  and  HI) 

In  its  pristine  state,  FeCls  undergoes  a  magnetic  phase  transition  at  8K  [15].  Measure¬ 
ments  of  the  magnetic  properties  of  stages  1  and  2  [la, 3],  prepared  from  HOPG,  revealed 
that  the  samples  undergo  a  magnetic  phase  transition  with  long  range  order,  however  the 
transition  temperatures  are  depressed  from  those  of  the  pristine  FeCls.  The  measurement 
of  Millman  et  al  [la]  was  performed  by  means  of  the  Mossbauer  effect,  while  Millman 
et  al  [3]  (see  appendix  I)  measured  the  magnetic  susceptibility  of  the  system,  the  latter 
method  has  a  better  temperature  resolution.  Both  measurements  found  that  stage  1  un¬ 
dergoes  a  phase  transition  at  4.3K  while  stage  2  undergoes  a  transition  at  1.3K.  Moreover, 
the  Curie- Weiss  $  indicated  that  stage  1  orders  antiferromagnetically  in-plane  and  ferro- 
magnetically  between  planes  while  stage  2  orders  antiferromagnetically  both  in-plane  and 
between  planes. 

A  Mossbauer  investigation  of  higher  stages  [1]  revealed  that  no  long-range  order  ex¬ 
isted  in  stages  4  and  6  down  to  temperatures  of  65mK.  Subsequent  susceptibility  mea¬ 
surements  on  those  samples  and  others  revealed  a  magnetic  susceptibility  m«rimnm  which 
occured  in  all  the  samples  at  a  temperature  between  1.7K  and  1.8K,  depending  on  stage 
[2a].  Although  the  maxima  occured  seemingly  at  the  same  temperature  for  all  stages,  the 
size  of  the  maximum  varied  widely  from  stage  to  stage  becoming  generally  greater  as  the 
stage  increased.  In  fact  our  measurements  reveal  that  the  size  of  the  maximum  increases 
by  a  factor  of  55  from  stage  1  to  stage  5,  (see  appendix  Q).  Moreover  the  susceptibility 


maximum  wu  supressed  by  relatively  low  magnetic  fields,  of  the  order  of  10G.  Since  the 
maximum  increased  in  intensity  with  stage,  we  assume  that  it  is  a  consequence  of  the 
two-dimensionality  of  the  system.  Within  the  same  stage,  the  height  of  the  susceptibility 
peak  increased  with  the  number  of  iron  atoms  next  to  iron  vacancies  [16], (see  appendix 

m). 

Phase  transitions  in  other  GIC  such  as  C0CI3,  and  NiClj  were  also  reported  by  Ka¬ 
rimov  et  al  [17],  recently  these  transitions  were  also  confirmed  by  Elahy  et  al  [2b].  The 
authors  of  reference  2b  have  concluded  that  the  magnetic  interc slants  exhibit  a  very  gen¬ 
eral  type  of  magnetic  behavior,  independent  of  species  or  of  stage.  More  recently  low 
temperature  magnetic  phase  transitions  were  reported  for  MnCla  [18],  and  also  for  stage 
1  C«Eu  [10].  Qualitatively  moet  of  the  reported  phase  transitions  have  shown  similarties 
in  both  magnetic  field  and  temperature  dependence,  differences  exist  only  in  the  loca¬ 
tion  of  the  transition  points  and  the  dc  magnetic  field  needed  to  quench  the  peaks  in  the 
measured  magnetic  susceptibility.  Most  of  the  reported  experimental  data  have  been  anal¬ 
ysed  on  the  basis  of  a  two-dimensional  XY  model  with  a  low  temperature  transition  to  a 
two-dimensional  antiferromagnetic  phase  [19,20]. 

Although  the  details  are  given  in  appendix  Q,  the  highlights  of  our  measurements 
were: 

1)  The  size  of  the  magnetic  susceptibility  peak  increases  exponentially  with  stage  up  to 
stage  5. 

2)  The  temperature  of  the  maxima  decrease  from  stage  1  to  stage  4  and  then  increase  after 
that.  All  the  maxima  occur  between  1.7K  and  1.8K. 

3)  At  any  stage  the  size  of  the  maximum  in  the  susceptibility  decreases  with  the  application 
of  a  small  magnetic  field,  the  size  of  the  decrease  depends  on  the  direction  of  the  field. 

4)  The  magnetic  field  shifts  the  temperature  at  which  the  maximum  occurs  to  a  higher 
temperature. 
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5)  At  high  fields,  the  susceptibility  starts  to  increase. 

6)  The  susceptibility  as  a  function  of  temperature  obeys  the  general  power  law  of  a  second 
order  phase  transition. 

• 

m.  Spin-lattice  Relaxation  (Appendix  IV) 

The  technique  used  to  measure  the  magnetic  susceptibility  allowed  us  to  measure 
both  the  inphase  and  out-of-phase  components  of  the  susceptibility.  The  relation  of  the 
in-phase  to  the  out-of-phase  component  at  various  frequencies  is  related  to  the  spin  lattice 
relaxation  time  of  the  system.  Appendix  IV  gives  the  details  of  the  measurement  and 
their  analysis.  The  main  conclusion  is  that  the  maximum  in  the  magnetic  susceptibility  is 
accompanied  by  a  corresponding  maximum  in  the  relaxation  times  and  in  the  specific  heat. 
The  maxima  in  the  relaxation  times  come  at  a  temperature  which  is  lower  than  that  of  the 
susceptibility  maximum  by  as  much  as  25mK.  Although  we  first  suspected  the  electrical 
conductivity  to  be  the  culprit  in  the  relaxation  maximum,  direct  measurements  of  the 
conductivity  revealed  a  smooth  and  continuous  behavior  with  no  evidence  of  anomalous 
behavior  in  the  temperature  region  of  the  susceptibility  maximum,  (see  following  section). 

IV.  Electrical  Resistivity  (Appendix  V) 

We  have  measured  the  c-axis  electrical  resistivity  of  the  various  stages  of  FeCl3  GIC. 
Appendix  V  gives  the  details  of  the  measurements.  We  have  measured  the  electrical 
resistivity  by  means  of  the  conventional  four  probe  technique.  The  temperature  range 
of  these  measurements  was  between  IK  and  300K.  At  room  temperature  stage  5  has  the 
highest  resistivity.  Stage  5  is  also  the  pivotal  stage  where  in  stages  higher  than  that  the 
c-axis  resistivity  decreases  with  temperature  (this  includes  unintercalated  HOPG),  while 
in  the  lower  stages  it  increases.  Thus  stage  5  seems  to  be  the  dividing  stage  between  the 
metalic  behavior,  low  stages,  and  activated  behavior,  high  stages.  The  behavior  of  the  c- 


axis  resistivity  is  in  qualitative  agreement  with  the  theory  of  Sugihara[2l]  which  attributes 
the  c-axis  conductivity  to  a  phonon  assisted  hopping  mechanism.  We  have  also  measured 
the  in-plane  resistivity  of  selected  samples  of  FeCl3  GIC’s  at  selected  temperatures  in 
connection  with  our  galvanomagnetic  measurements,  (see  following  section). 

V.  Galvanomagnetic  Measurements,  (Hall  Effect,  Magnetoresistance  and  the  Shubnikov-de 
Haas  Effect)  (Appendix  VI  and  VH) 

Again,  the  details  of  these  measurements  are  given  in  Appendix  VI  and  VII.  The  main 
results  of  these  measurements  are: 

1)  The  conduction  in  stage  2  and  5  is  dominated  by  hole  carriers,  while  stage  9  and  HOPG 
exhibit  both  electron  and  hole  conduction  with  the  dominant  mechanism  depending  on 
temperature  and  field. 

2)  The  magnetoresistance  is  positive  in  all  stages  and  at  low  temperatures  exhibits  quantum 
oscillations. 

3)  The  quantum  oscillations  in  the  magnetoresistance  persist  up  to  temperatures  as  high 
as  30K. 

4)  The  Fourier  analysis  of  the  data  indicates  that  the  oscillation  instage  5  is  dominated  by 
a  single  frequency  contributed  by  holes,  while  in  stage  9  several  frequencies  were  observed. 

VI.  Theoretical  Calculations  (Appendix  VHI  and  IX) 

A  theoretical  effort  was  started  after  the  experimental  investigation  had  yielded  in¬ 
teresting  results.  At  the  beginning  it  was  aimed  at  the  explanation  of  some  of  our  results. 
After  a  while,  however,  this  investigation  became  of  interest  in  itself.  The  model  explored 
was  that  of  a  two  dimensional  array  of  dipoles  situated  on  a  honeycomb  lattice.  A  clas¬ 
sical  dipolar  potential  was  assumed,  and  an  exact  ground-state  energy  was  calculated.  A 
mean-field  calculation  was  performed  to  determine  the  dependence  of  the  system  on  tem¬ 
perature.  So  far  only  the  planar  lattice  was  considered  and  calculations  of  the  influence  of 


boundaries  and  defects  in  the  array  are  currently  in  progress. 

The  details  of  the  calculations  and  their  results  are  given  in  Appendix  VUI  and  DC. 
The  major  results  of  the  calculations  are: 

1)  There  is  an  infinitely  degenerate  ground  state. 

2)  If  one  uses  the  spacing  and  the  dipole  moments  measured  for  the  iron  ions  in  FeCl3  the 
energy  of  the  ground  state  becomes  -0.98K. 

3)  There  is  a  critical  point  near  IK  in  zero  field. 

4)  We  are  able  to  map  out  a  phase  diagram  in  the  magnetic  field  -  temperature  plane  with 
three  phases.  At  low  temperature  and  field  there  is  the  ordered  phase,  at  low  temperature 
and  high  field  there  is  a  field  induced  ferromagnetic  phase,  and  at  high  temperature  we 
obtain  a  paramagnetic  phase. 

5)  In  the  low  and  moderate  field  region,  and  below  the  transition  temperature,  the  magnetic 
isothermal  susceptibility  increases  with  an  applied  external  field. 

6)  From  our  calculation  we  can  determine  the  exact  orientation  of  the  dipoles  as  a  function 
the  applied  magnetic  field  and  the  angle  of  the  field  with  respect  to  the  lattice. 

Further  calculations  to  determine  the  interaction  of  this  lattice  with  similar  lattices 
in  three  dimensions  above  and  below  it  are  in  progress.  We  also  propose  to  calculate 
the  influence  of  defects,  and  superlattices  on  the  behavior  of  the  system.  As  shown  in 
Appendix  DC,  some  of  this  work  has  already  been  started. 

VII.  Instrumentation  (Appendix  X) 

In  connection  with  our  experimental  work,  we  have  developed  a  self-balancing  re¬ 
sistance  bridge  which  uses  a  novel  feedback  system.  The  bridge  is  appropriate  for  low 
temperature  measurements  because  of  its  low  power  dissipation.  The  details  of  the  bridge 
are  given  in  appendix  X.  We  intend  to  develop  more  self  balancing  instrumentd,  applying 
this  feedback  system.  This  system  is  fully  compatible  with  computer  usage. 


8 


Characterization  of  Samples 


The  samples  which  were  used  in  the  above  measurements  were  characterized  for  stag¬ 
ing  by  means  of  001  x-ray  analysis  and  mixed  stage  samples  were  not  considered.  In  many 
cases  the  samples  were  also  examined  by  means  of  the  Mossbauer  effect  which  could  shed 
light  on  the  in-plane  density  of  the  intercalant  and  the  ratio  of  Fe2+  to  Fe3+. 


Summary 

The  picture  which  emerges  from  our  investigations  is  that  we  are  indeed  observing 
the  behavior  of  a  two  dimensional  magnetic  system.  The  low  dimensionality  is  brought 
about  by  the  shielding  of  the  intercalant  layers  by  the  intervening  graphite  layers.  The 
shielding  also  depends  on  the  charge  transfer  between  the  intercalant  and  the  graphite.  In 
graphite  there  is  a  band  near  jfc  =  0  about  half  filled  with  electrons [10].  The  intercalation 
of  an  acceptor  compound  will  draw  off  those  electrons  so  that  the  electrical  conduction 
is  hole  like.  Thus  low  stages  exhibit  a  positive  Hall  coefficient.  The  shielding  of  the 
intercalant  layers  is  due  to  electromagnetic  shielding  analogous  to  the  skin  effect,  where 
the  skin  depth  is  a  function  of  the  electrical  conductivity,  which  in  turn  is  a  function  of 
the  number  of  charge  carriers.  Stages  1  to  5  almost  empty  the  band  of  electrons  and  thus 
the  graphite  layers  have  maximum  conductivity.  Above  stage  5,  the  quantity  of  acceptors 
is  not  large  enough  to  absorb  all  the  electrons  which  can  be  donated  by  the  graphite,  and 
the  conductivity  of  the  graphite  layers  decreases.  Moreover,  the  conductivity,  according 
to  our  Hall  measurements,  is  due  to  both  positive  and  negative  carriers.  This  gives  rise 
to  decreased  magnetic  shielding  of  the  intercalant  layers,  with  a  corresponding  decrease  in 
the  peak  height  in  the  susceptibility  at  1.7K-  Thus,  the  transition  at  1.7K  can  become  a 
very  sensitive  probe  of  the  internal  magnetic  field,  the  charge  transfer,  and  the  electrical 
conductivity  of  the  graphite  layers. 

The  nature  of  the  magnetic  phase  transition  is  still  not  definitely  determined,  although 


some  facts  have  been  established:  1)  The  transition  is  definitely  due  to  the  two  dimensional 
nature  of  the  system;  2)  Dipolar  interactions  play  a  significant  role  in  the  transition.  Our 
theoretical  efforts  are  directed  towards  the  explanation  of  the  nature  of  the  low  temperature 
phase  and  so  far  we  have  managed  to  theoretically  reproduce  some  of  the  effects  observed 
experimentally.  Those  effects  are  the  increase  of  the  susceptibility  in  the  magnetic  field 
at  fields  above  40G,  which  agrees  with  our  measurements  which  are  given  in  appendix  II 
and  calculated  in  appendix  Vm,  and  the  decrease  at  extremely  low  fields  when  defects  are 
present  as  seen  in  appendix  II  and  calculated  in  appendix  DC.  It  is  evident  that  dipolar 
interactions  do  not  tell  the  whole  story,  and  that  exchange  interactions  will  also  play  a 
role.  However,  the  degeneracy  in  the  ground  state  of  our  dipolar  model  as  well  as  the 
extreme  sensitivity  of  the  dipole  configuration  to  the  direction  of  the  magnetic  field  give 
us  a  hint  that  the  material  may  resemble  a  two  dimensional  spin-glass.  This  will  be  the 
subject  of  further  investigations. 
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Measurement  of  eha  ae  auacapcibllity  of  wall  charaetarlzad  farrlc 
chloride  Intercalated  graphito  compounds  both  parallal  and  parpandl- 
cular  to  tha  baaal  plana,  ahowa  that  the  stage  1  and  2  compounds 
undergo  a  magnetic  phase  transition  at  4.3  *  0.2  K  and  1.3  t  0.2K 
respectively. 


1.  INTRODUCTION 

Tha  structure  of  graphite  Intercalation 
compounds  indicates  thnt  both  the  graphitic  and 
intercalate  layers  retain  the  identity  nf  their 
parent  oeterlale. 1  For  magnetic  Intercalants 
this  obeervat  Ion  suggeats  that  thaae  compounds 
should  show  magnetic  phase  transitions.  Thus. 
Intercalation  can  be  used  to  naps rate  the  mag¬ 
netic  layers  by  a  controlled  number  of  d  1  nomi¬ 
nee  Ic  layers,  gy  comparing  the  magnetic  proper¬ 
ties  of  tho  parent  magnetic  material  sad  the 
different  stage  compounds,  one  can  lnsna>lgate 
the  Influence  nf  interlayer  interactions  on  the 
■agnatic  properties  of  tho  system. 

Magnetic  susceptibility,  neutron  diffrac¬ 
tion.  and  Moaabauer  spectroscopy  data  shorn  that 
pristine  anhydrous  FeCl ,  undergoes  a  nagnatlc 
phase  transit  too  at  about  8.8  K.  However,  the 
reported  magnetir  properties  in  graph! te-FeCl j 
have  shown  discrepancies  which  nay  possibly  be 
sample  dependent.  ~ 1 1  For  thn  nagnatlc  suscep¬ 
tibility  neanursnents,  there  are  discrepancies 
nor  only  In  the  transition  temperature  and  char¬ 
acter  of  a  nagnatlc  phase  transition, •“ 18  but 
also  In  tho  very  existence  of  a  magnetic  state. 
However,  it  Is  not  surprising  to  find  such  dts- 
rreuanc toe  in  the  literature  due  to  the  uncher- 
arterized  differences  in  iron  vacancies,  chlor¬ 
ine  vacancies,  FeCl y  vseancios,  and  VeClj  im¬ 
purities  which  are  likely  to  vary  from  sample 
to  scmplo. 

tn  this  work,  wo  find  a  susceptibility  peak 
in  both  the  s  to  go  1  and  2  graphlte-FeCl j  com¬ 
pounds  similar  in  functional  fora  to  that  ob¬ 
served  In  pristine  anhydrous  FeClv*  Wo  show 
that  the  stage  i  and  2  compounds  undergo  a  mag¬ 
netic  phase  transition  at  4.3  1  0.2  K  and  1.3  t 
0.2  K  respectively.  Thus  the  Interactions  be¬ 
tween  mngnatlc  planes  seem  to  etrongly  Influence 
the  magnetic  order  observed  tn  these  compounds. 
In  addition  to  x-ray  analysts  we  have  also  char¬ 
acterised  our  samples  by  Mossbsuer  spectroscopy 
to  determine  tho  amount  of  Iron  vacancies,  chlo¬ 
rine  vacancies,  FeClj  vacancies,  and  FeClj  im¬ 
purities. 

2.  EXPERIMENTAL 

The  samples  studied  were  prepared  by  inter¬ 
calating  anhydrous  FeClj  into  highly  ordered  py¬ 


rolytic  graphite  (H0PC)  using  the  two-zone  vapor 
transport  method.  Employing  this  growth  process, 
the  graphite  temperature  wee  maintained  at  350°C 
while  the  temperature  of  the  ferric  chloride  was 
varied  to  achteva  different  stnges.  The  Inter¬ 
calated  compounds  were  characterized  for  Identi¬ 
ty  and  uniformity  of  staging  using  (OOt)  x-ray 
dtf fraction  both  before  and  after  our  magnetic 
susceptibility  measurements.  Both  samples  were 
essentially  single  etaged,  end  showed  less  than 
20X  end  )1  of  secondary  diffraction  poeks  from 
both  othor  stages  end  free  graphite  in  the  stage 
1  end  2  compounds  respectively. 

These  samples12  ware  further  checked  for 
fidelity  of  tn-plane  density  by  taking  their 
Moeabauar  spoctra  at  100  K.  For  the  stage  I 
sump Is  st  this  temperature  the  Mnsabauer  peak 
can  be  fit  with  only  a  single  Lorentztan  line 
whose  llnevldth,  0.38  rnm/sec,  and  lsoaar  shift, 
0.3b  ims/sec,  ars  essentially  the  same  ns  those 
of  pristine  anhydrous  FeClj*’.  This  indicates 
that  there  te  neither  a  change  In  site  aymmetry 
at  the  iron  lone  nor  is  there  any  appreciable 
format  ton,  l.e.  less  than  42,  of  either  iron 
vacancies,  chlorine  vacancies,  FeClj  vacancies, 
or  FoCl;  Impurities  In  this  sample.  An  analy¬ 
sis  of  tha  MDaabauer  spectrum  of  the  stage  2 
sample  indicates  that  9T  of  tha  sites  are  near¬ 
est  neighbors  to  Iron  vacancies.  Tills  identi¬ 
fication  will  be  discussed  in  detail  elsewhere.  *'• 
Assuming  these  vacancies  are  rnndomly  situated, 
this  implies  that  only  approximately  32  of  tha 
Iron  lattice  aitee  contain  vacancies. 

The  susceptibility  measurements  were  made 
by  a  standard  ac  bridge  method15  operating  at 
about  40  hz.  The  amplitude  of  the  magnetic 
field  at  the  sample  was  kept  below  1  Oe.  The 
coils  were  kept  at  a  constant  temperature  while 
tha  temperature  of  the  samples  wee  varied  inside 
an  isolated  chamber.  The  temperatures  above  4.2 
K  were  measured  hy  a  calibrated  carbon  resistor 
and  •  calibrated  metal  resistance  thermometer. 
Below  4.2  X  the  sample  was  in  direct  contact 
with  the  helium  bath  and  temperatures  were  de¬ 
termined  fro*  the  vapor  pressure  of  the  helium 
bath. 

3.  RESULTS 

In  Figure  i,  we  show  the  susceptibility 
curves  of  the  stage  I  and  2  camples  with  the 
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lipin  I.  iwc«ytibili(y  (I)  versus  toptntui 
(T),  plotted  on  a  sani-logarithnic 
scale,  for  tha  stage  I  and  2  coapounda. 
Tha  naasuring  field  waa  fixed  both  par¬ 
allel  (e-ax ia)  and  perpendicular 
(c-axis)  to  the  basal  plane.  Notice 
the  snail  ensue  t  of  stage  I  (undetec¬ 
table  ir  the  x-ray  diffractograne) 
which  can  be  detected  in  the  atage 
2  curve. 

neaeuring  field  both  parallel  (a-axla)  end  per¬ 
pendicular  (c-axia)  to  tha  baaal  plane  over  the 
tenperatura  range  of  1  K  to  200  K.  The  auacap- 
tibillty  ia  glean  in  arbitrary  unlta  which  have 
been  no rnall zed  for  the  anount  of  iron  la  each 
eaaple  by  aeawing  the  etolchlonatry  C(.M  PeCl  . 
However,  we  note  that  a  coeparlsoe  of  areas  1 
under  the  MBaebauer  curves  tap  lies  that  as 
auch  aa  SOX  of  tha  bulk  stags  1  sanpla  nay 
have  consisted  of  un Intercalated  rex ions. 

This  estlaate  could  be  nade  to  irithln  13Z. 

Coupled  with  Mosabaoer  inf ornatlon  that  tha 
electronic  spine  for  both  stag-s  are  point las 
in  the  basal  plans.11 

several  laportaat  results  can  be  dodeced  free 
the  data.  Por  tha  stage  1  conpound  an  antl- 
ferronagnetlc  susceptibility  peak  is  observed 
whose  nsxlnun  occurs  at  0.3  K  in  both  tha  par¬ 
allel  and  perpendicular  susceptibilities.  As 
expected  for  the  electronic  spina  pointing  in 
the  basal  plane,  tha  susceptibility  is  greater 
in  the  direction  perpendicular  to  the  planes. 

Also,  at  tenperatures  lass  than  X,  the  para¬ 
llel  susceptibility  is  clearly  decreasing  while 
the  perpendicular  susceptibility  is  approx lante- 
ly  constant.  Por  the  stage  2  conpound  the  naxl- 
nun  of  the  antlfarrouagwatlc  peak  has  shifted 
down  to  1.72  R  and  tha  peak  ia  tbs  perpendicu¬ 
lar  direction  ia  drastically  reduced.  This 
indicates  that  the  nagnatie  nonants  are  nore 
strongly  bound  ia  tha  basal  plane  In  stags  2 
than  in  atage  1.  A  eonparisen  of  the  parallel 
susceptibility  of  bath  stagee,  platted  on  a 
logarlthnlc  scale  ia  Figure  1  shews  that  the 
curves  are  qualitatively  similar  ia  reduced 
units,  T/T-. 

The  Mel  tenperatura,  TH,  for  taro  field 
susceptibility,  is  given  by  a  navi aim  in 
|d(Tx)/dT|  . 11  The  value  of  Tv,  obtained  for 
tha  stage  1  conpound  is  4.3  i  "6.2  K  for  ori¬ 
entations  both  parallel  and  perpendicular  to 
the  basal  plane.  Por  the  stags  2  conpound,  we 
obtain  1.3  t  0.2  R  for  both  parallel  and 
perpendicular  to  tha  baaal  plane.  Thus  the 
transition  tenperatura  decreases  by  approxl- 
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Figure  2.  Inverse  susceptibility  (if1)  versus  tear 
perature  (T)  for  the  curves  shown  in 
Figure  I.  Tha  solid  lines  represent 
least  squares  fits  to  the  data. 

nately  a  factor  of  three  when  the  lnterpla- 
nar  spacing  is  Increased  fraa  9.4  X  (stage  1) 
to  12.6  X  (stage  2). 

Aa  shown  in  Figure  2,  the  stage  1  con- 
pound  obeys  the  Curla-Welss  law  In  the  tenp- 
arature  region  above  13  K.  We  obtain  the 
values  of  6  ■  -3.E  t  0.5  K  for  tha  direction 
parallel  to  the  planes  and  9  *  >3.8  ♦  0.3  K 
for  tha  direction  perpendicular  to  the  planes. 
Although,  because  of  tha  large  anisotropy 
energies  it  is  difficult  to  relate  9  to  the 
exchanged  interactions,  these  results  are 
consistent  with  the  transition  tenperatura. 

Tha  values  of  the  Curla-Welss  constants,  9, 
nay  Indicate  that  tbs  donlnant  exchange 
interaction  is  antlferronagnetlc  within  the 
layers  while  it  is  ferronagnetlc  between  tha 
layers.  Previous  results  have  Indicated  that 
the  Interlayer  Interaction  was  antlferronag- 
netlc.1 ’*  Por  tha  stage  2  conpound,  we 
obtain  tha  result  9  “  -7.6  t  1.0  E  and 
9  -  -33  t  2.0  R  in  the  direction  parallel 
and  perpendicular  to  the  layers  respectively, 
and  the  Curla-Welss  law  is  followed  above  5  K 
as  shown  in  Figure  2.  These  values  of  9  are 
Markedly  different  Iron  the  TH  values  obtained 
by  finding  the  naxinun  in  |d(fx)/dT|.  The 
transition  tenperatures,  T_,  values  of  9  and  the 
Curia  constants  are  llsted"ln  Table  1. 

In  order  to  corroborate  our  results  for 
the  stage  2  conpound,  we  Measured  the  a-axls 
susceptibility  of  two  other  stage  2  compounds . 

The  susceptibility  of  one  of  the  easplea  ex¬ 
hibited  the  sane  qualitative  behavior  as  the 
stage  2  conpound  reported  here.  Its  MBsabauar 
spec t run  was  also  essentially  the  sane  as  the 
stage  2  conpound  reported  here,  except  that 
slightly  nore,  UZ,  of  the  sites  were  nearest 
neighbors  to  iron  vacancies.  The  susceptibil¬ 
ity  curve  of  the  other  stage  2  sanpla  had  both 
a  peak  at  the  sane  tenperatura,  7.1  K,  and  the 
sans  qualitative  behavior  as  the  stage  2  con¬ 
pound  of  Karlnov  at  al.*  However,  its  MBse- 
bauar  spec true  indicated  that  at  least  nost, 
and  probably  all  of  this  sanpla  was  composed 
of  FeCl}.  We  note  that  this  sanpla  was  pro¬ 
duced  In  the  sane  apparatus  and  under  the  sane 
growth  conditions  as  all  the  other  samples. 

These  results  underscore  the  need  for  both  a 
structural  and  chenieal  characterisation  of 
these  coapounds  In  addition  to  the  standard 
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TABLE  1  This  table  Hat  a  the  magnetic  properties  of  the  compounds. 

CURIE 

CONSTANT 


STAGE 

DIRECTION 

TRANSITION 

TEMPERATURE 

SUSCEPTIBILITY 

MAXIMUM  0 

(ARBITRARY 

UNITS) 

NEAREST  NEIGHBORS 
TO  IRON  VACANCIES 

1 

a-axla 

4.3 

6.5 

-3.8 

222 

S4X 

1 

c-axis 

4.3 

6.5 

+3.8 

132 

*4X 

2 

a- axis 

1.3 

1.72 

-7.6 

242 

91 

2 

c-axis 

1.3 

1.72 

-33.0 

575 

9X 

characterisation  by  (OOt)  x-ray  dlf fractocraas . 

The  data  at  temperatures  below  2.2  K,  in  the 
curves  of  Figure  1,  were  all  aeasuted  in  a  dc 
field  of  80  Oe.  As  seen  in  Fig..  3,  a  small  narrow 
peak  is  observed  at  1.7  K  when  no  external  dc 
field  is  applied.  This  peak  also  appears  at 
this  temperature  in  higher  stages  but  is  re¬ 
duced  by  approximately  a  factor  of  SO  in  the 
low  stages  (stage  1  nnd  2)  relative  to  the 
high  stages  (stage  4  and  6). *  It  can  be  al¬ 
most  completely  quenched  by  fields  as  small 
as  20  Oe.  Similar  behavior  was  observed  in 
both  NlClj  and  C0CI2  intercalated  graphite 
compounds.17  Because  of  the  coincidence  of 
this  peak  and  that  due  to  antiferromagnetism 
in  the  stage  2  compound,  wa  suppressed  it  in 
an  80  Oe  field  In  order  to  observe  the  full 
features  of  the  anti ferromagnetic  transition. 

The  occurence  of  the  small  peek  at  a  lower 
temperature  in  the  stage  1  sample  than  the 
three-dimensional  magnetic  transition  lndl- 


Flgure  3.  Subtraction  of  the  susceptibility  at 
80  Oe  from  that  at  sero  field.  For 
clarity,  the  susceptibility  scale  has 
been  expanded  by  a  factor  of  3  from 
that  of  Figure  I.  Mote,  within  the 
accuracy  of  our  experiment  the  stags 
I  c-sxis  susceptibility  was  the  same 
both  at  sero  field  and  at  80  Os  and 
is  therefore  not  plotted. 


cates  that  it  any  be  caused  by  high  stage 
impurities  la  the  predominately  low  stage 
samples. 

4.  DISCUSSION 

Previous  susceptibility  results  have  been 
interpreted  as  showing  that  even  for  stage  i 
graphlta-FeCl)  the  Interlayer  interactions  are 
negligible  compared  with  the  Intralayer  inter¬ 
actions.  However,  our  results  Indicate 
that  the  magnetic  order  seen,  at  least  In  the 
stags  I  compound,  in  strongly  influenced  hy 
the  Interlayer  interactions.  In  fact,  we  have 
shown  that  the  transition  temperature  decreases 
by  a  factor  of  2  la  going  from  pristine  anhy¬ 
drous  FeCl),  T_  «  8.8  E2-1*  to  stage  1  and  by 
an  additional  Factor  of  ]  in  going  from  stage 
1  to  stage  2.  Both  the  parallel  and  perpen¬ 
dicular  susceptibilities  follow  the  Cur  le¬ 
ws  lss  law,  and  an  Increase  in  either  the  ac 
measuring  frequency  or  the  amplitude  of  the 
sc  field  by  about  an  order  of  magnitude,  or 
an  applied  dc  field  of  160  Oe  causes  no  change 
in  the  transition  temperature. 

The  small  susceptibility  and  the  quite 
largs  value  of  8  in  the  direction  porpendl- 
culsr  to  the  basal  plana  in  the  stage  2  com¬ 
pound  nay  be  due  to  the  presence  of  iron 
vacancies  in  this  compound.  Iron  vacancies 
would  crests  an  electric  field  gradient  in 
the  basal  plane  which  would  causa  anisotropy 
between  the  ln-plane  and  out-of-plane  direc¬ 
tions. 

We  are  happy  to  acknowledge  many  help¬ 
ful  discussions  with  Professors  C.  K lrc renow, 
and  S.  Rednsr.  Our  deep  gratitude  is  expressed 
to  Dr.  R.B.  Frankel  and  Dr.  C.  Papaefthymlou 
for  the  use  of  their  MOssbauer  equipment,  to 
Dr.  A.W.  Moore  of  Union  Carbide  for  donating 
the  HO PC,  and  to  Dr.  H.  Hasurek,  Dr.  C. 
Nlcollnl,  and  A.  Thompson  for  growing  and 
characterizing  one  of  the  samples.  This  work 
was  supported  in  part  by  AP0SR  contract 
#F4 96 20-8 1-0006. 
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APPENDIX  II 

MAGNETIC  PROPERTIES  OF  GRAPHITE  INTERCALATION  COMPOUNDS* 

A.K.  Ibrahim  and  George  O.  Zimmerman 
Physics  Department,  Boston  University 
Boston,  MA  02215 

Comprehensive  measurements  of  the  magnetic  properties  of  FeCU  graphite  interca¬ 
lation  compounds  are  presented  with  an  emphasis  on  the  low  temperature  region  where 
a  susceptibility  murimum  is  observed  in  all  stages.  This  maximum,  which  varies  in  size 
according  to  stage,  occurs  In  a  very  narrow  temperature  range  and  is  attributed  to  the 
two-dimensionality  of  the  intercalate  system.  It  obeys  the  power  law  of  a  second  order 
phase  transition  with  an  exponent  7  which  has  a  value  between  1.8  and  2.0.  The  max¬ 
imum  occurs  only  in  the  in-plane  direction  with  no  corresponding  e-axis  susceptibility 
response.  The  application  of  an  external  magnetic  field  drastically  supresses  the  suscep¬ 
tibility  m«ritntnw  and  shifts  it  to  higher  temperatures.  Both  in  plane  and  out  of  plane 
measurements  are  presented  and  the  magnetic  properties  of  stages  1  through  6  as  well  as 
stage  0  are  compared. 

: 

•  Supported  by  Air  Force  Office  of  Scientific  Research  Grant  AFOSR  82-0288.  j 
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INTRODUCTION 


Graphite  intercalation  compounds  (GIC)  are  naturally  layered  [1]  and  are  of  inter¬ 
est  because  electronically,  only  the  graphite  layers  next  to  the  intercaiant  change  signif¬ 
icantly  upon  intercalation  [2j.  For  a  magnetic  intercaiant,  the  inner  layers  can  be  a  con¬ 
trolling  factor  in  the  magnetic  interactions  between  different  magnetic  intercaiant  layers. 
The  systematic  change  of  the  number  of  graphite  layers  between  two  adjecent  intercaiant 
magnetic  layers  can  provide  two-dimensional  magnetic  systems  by  shielding  the  inter¬ 
caiant  layers  from  each  other  until  their  interaction  is  insignificant. 

Magnetic  compounds  have  been  widely  used  in  the  intercalation  of  highly  ordered 
pyrolytic  graphite  (HOPG)  and  generally  one  observes  magnetic  anomalies  in  the  GIC, 
even  in  high  stages  (stage  index  refers  to  the  number  of  graphite  layers  between  two  con¬ 
secutive  intercaiant  layers).  It  has  been  reported  [3]  that  for  all  the  magnetic  compounds 
which  have  a  low  temperature  phase  transition  there  exists  a  corresponding  transition 
when  these  compounds  are  intercalated  into  graphite.  Additional  transitions  due  to  the 
magnetic  two  dimensionality  of  the  systems  were  also  reported  [4,5].  Thus  intercalation 
of  magnetic  species  into  graphite  provides  very  useful  model  systems  for  the  study  of 
two-dimensional  magnetism. 

In  this  work  we  present  a  systematic  study  of  the  low  temperature  phase  transition 
of  FeCl3  intercalated  graphite  compounds  and  its  behavior  as  a  function  of  the  stage  of 
the  sample.  We  also  investigate  the  dependence  of  these  transitions  on  an  external  dc 
magnetic  field.  Our  data  for  FeCl*  GIC  are  compared  with  reported  data  for  similar 
GIC  sytems.  The  standard  low  frequency  ac  susceptibility  bridge  technique  [A]  is  used 
to  obtain  the  experimental  data  . 

The  crystal  structure  of  FeCl*  is  a  repeated  sequence  of  chlorine,  iron,  and  chlorine 
hexagonally  arranged  atoms  [7],  such  that  every  iron  atom  is  surrounded  by  an  octa- 
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hedron  of  chlorine  atoms.  The  iron  atoms  in  subsequent  layers  are  displaced  with  re¬ 
spect  to  each  other  by  1/3  of  the  unit  cell  and  thus  three  FeCl3  layers  are  encompossed 
by  the  cell  vector  in  the  C  direction.  Upon  intercalation  of  FeCl3  into  the  graphite,  one 
electron  is  donated  by  the  graphite  host  for  every  four  iron  atoms  [8].  It  is  not  yet  clear 
what  site  these  donated  electrons  occupy,  room  temperatures  Mossbauer  data  [0,10]  and 
Raman  spectroscopy  [11]  do  not  support  the  existence  of  FeCl3  as  a  result  of  the  accep¬ 
tance  by  FeCl3  of  the  donated  electron  [12].  However,  there  is  a  discrepancy  in  the  low 
temperature  Mossbauer  data  for  FeCI3.  Miliman  and  Kirczenow  [13]  have  reported  the 
existence  of  Fe++  ions  at  temperatures  as  high  as  100K  which  is  contradictory  to  the 
work  of  Ohhashi  et  al  [10]. 

Magnetic  susceptibility  measurements  on  FeCl3  graphite  intercalated  compound 
were  performed  by  Karimov  et  al  [14]  and  Hohlwein  et  al  [15]  in  the  70’s.  The  interest 
there  was  the  nature  of  magnetic  behavior  of  the  FeCls  which  in  pristine  state  under¬ 
goes  a  magnetic  phase  transition  at  8K  [16].  Those  measurements  concentrated  mainly 
on  the  properties  of  stages  1  and  2.  Although  the  results  of  those  measurements  dis¬ 
agreed  on  the  details  the  authors  reported  a  phase  transition  in  the  temperature  re¬ 
gion  of  their  measurements .  The  method  used  in  those  measure  ments  was  the  Faraday- 
ballance  supplemented  by  an  ac  method,  the  samples  consisted  of  small  crystallites. 
More  recent  measurements  of  the  magnetic  properties  of  stages  1  and  2  [17,18],  pre¬ 
pared  from  HOFG,  revealed  that  the  samples  undergo  a  magnetic  phase  transition  with 
long  range  order,  however  the  transition  temperatures  are  depressed  from  those  of  the 
pristine  FeCls-  The  measurement  of  Miliman  et  al  [17]  was  performed  by  means  of  the 
Mossbauer  effect,  while  Miliman  et  al  [18]  measured  the  magnetic  susceptibility  of  the 
system,  the  latter  method  has  a  better  temperature  resolution.  Both  measurements 
found  that  stage  1  undergoes  a  phase  transition  at  4.3K  while  stage  2  undergoes  a  tran¬ 
sition  at  1.3K.  Moreover,  the  Curie- Weiss  #  indicated  that  stage  1  orders  antiferromag- 
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netically  in-plane  and  ferromagnetically  between  planes  while  stage  2  orders  antiferro- 
magnetically  both  in-plane  and  between  planes. 

A  Mossbauer  investigation  of  higher  stages  [19]  revealed  that  no  long-range  order 
existed  in  stages  4  and  6  down  to  temperatures  of  65mK.  Subsequent  susceptibility  mea¬ 
surements  on  those  samples  and  others  revealed  a  magnetic  susceptibility  maximum 
which  oceured  in  all  the  samples  at  a  temperature  between  1.7K  and  1.8K,  depending 
on  stage  [20].  Although  the  maxima  oceured  seemingly  at  the  same  temperature  for  all 
stages,  the  size  of  the  maximum  varied  widely  from  stage  to  stage  becoming  generally 
greater  as  the  stage  increased.  In  fact  this  investigation  reveals  that  the  size  of  the  max¬ 
imum  increases  by  a  factor  of  55  from  stage  1  to  stage  5.  Moreover  the  susceptibility 
maximum  was  supressed  by  relatively  low  magnetic  fields,  of  the  order  of  10G,  which 
may  have  accounted  for  the  fact  that  it  was  not  seen  by  other  investigators  prior  to  that 
time.  Since  the  maximum  increased  in  intensity  with  stage,  it  was  assumed  that  it  was  a 
consequence  of  the  two-dimensionality  of  the  system.  Within  the  same  stage,  the  height 
of  the  susceptibility  peak  increased  with  the  number  of  iron  atoms  next  to  iron  vacancies 
[21]. 


Phase  transitions  in  other  GIC  such  as  CoCIa,  and  NiCla  were  also  reported  by  Ka¬ 
rimov  et  al  [22,23],  recently  these  transitions  ware  also  confirmed  by  Elahy  et  el  [24]. 

The  authors  of  [24]  have  concluded  that  the  magnetic  intercalsnta  exhibit  a  very  general 
type  of  magnetic  behavior,  independent  of  species  or  of  stage.  More  recently  low  tem¬ 
perature  magnetic  phase  transitions  were  reported  for  MnCla  [25],  and  also  for  stage  1 
C«Eu  [25].  Qualitatively  most  of  the  reported  phase  transitions  have  shown  similartics 
in  both  magnetic  field  and  temperature  dependence,  differences  exist  only  in  the  loca¬ 
tion  of  the  transition  points  and  the  dc  magnetic  field  needed  to  quench  the  peaks  in 
the  measured  magnetic  susceptibility.  Most  of  the  reported  experimental  data  have  been 
analysed  on  the  basis  of  a  two-dimensional  XY  model  with  a  low  temperature  transition 
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EXPERIMENTAL 


a)  Characterisation  of  samples 

The  FeCl3-GIC  samples  were  prepared  using  a  standard  two-zone  furnace  technique 
where  stage  index  was  controlled  by  the  temperature  difference  between  the  graphite 
host  (HOPG)  and  the  FeCls  powder.  The  samples  were  in  the  form  of  thin  rectangular 
plates  of  dimensions  1.5x0. 5x0.1  cm3.  Well-staged  samples  were  achieved  by  controlling 
the  pressure  of  Cla  gas  inside  the  intercalation  tube,  as  well  as  the  partial  pressure  of 
FeGb  through  rigid  temperature  control!  After  intercalation,  the  samples  were  charac¬ 
terized  for  identity  and  uniformity  of  staging  using  x-ray  (001)  diffraction.  Only  single 
stage,  well  staged  samples  were  used  in  the  magnetic  measurements.  The  x-ray  diffrac- 
tograms  were  also  used  to  determine  the  c-aads  repeat  distance  I«  after  cycling  the  sam¬ 
ples  from  room  to  liquid  helium  temperature,  and  showed  that  the  cycling  did  not  affect 
this  staging  distance. 

Most  of  the  samples  measured  were  characterised  by  means  of  the  MSssbauer  effect, 
which  can  reveal  the  ratio  of  iron  atoms  next  to  iron  vacancies  to  that  of  the  total  num¬ 
ber  of  iron  atoms  in  the  intercalant  layer  as  well  as  that  of  Fa3*  to  that  of  Fe3+  [17,19]. 
The  initial  measurements  were  performed  on  the  same  samples  used  in  the  Mdssbauer 
investigations.  The  same  measured  samples  were  characterised  periodically,  before  and 
after  the  measurements,  over  a  time  spaa  of  over  a  year,  having  been  stored  at  room 
temperature  la  a  dry  aitrogea  enviromeat.  No  deterioration  was  observed  in  the  sam¬ 
ples  during  the  period  of  meesuremeats  [27].  Once  deterioration  was  observed  the  sam¬ 
ples  were  discarded.  Ia  higher  stage  samples,  where  stage-disorder  is  expected  [28],  the 
Headricks-Teller  [29],  and  Mats  and  Hohlwcin  [30]  analysis  teehaique  were  used  to  cal¬ 
culate  the  intensity,  width  aad  location  of  the  x-ray  reflections.  We  find  that  our  exper¬ 
imental  x-ray  data  on  the  stage  9  samples  reported  ia  this  article  are  in  good  agreement 
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with  that  calculated  for  the  pure  and  well-staged  stage  9.  Although  a  small  admixture  of 
stage  10  can  not  be  excluded,  this  stage  serves  as  an  example  of  a  high  stage  sample. 


As  shown  in  Table  I  of  reference  (19),  the  Mossbauer  analysis  of  our  initial  samples 
showed  that  17±3%  of  the  iron  sites  were  next  to  vacancies  (ISNV)  in  stages  1  and  2 
while  this  number  was  19±3%  and  13±3%  for  stages  4  and  6,  respectively.  The  Fea+ 
sites  comprised  22±Z%  in  stages  1  and  2  at  low  temperatures  while  only  3±1%  in  stages 
4  and  0.  In  a  similar  set  of  measurements  it  was  also  shown  that  the  amount  of  Fe2+  in¬ 
creased  systematically  with  a  decrease  in  temperature  from  about  100K  to  10K,  and  was 
constant  below  that  temperature  [13].  In  some  of  the  subsequent  samples  the  percentage 
of  iron  atoms  next  to  vacancies  was  reduced  to  7%  or  lower. 

If  one  assumes  that  most  of  the  sites  near  to  vacancies  come  from  the  atoms  at  the 
boundary  of  the  intercalate  islands,  one  can  calculate  the  minimum  average  diameter  of 
an  island  to  be  94A  for  a  1996  ISNV  and  ITSA  for  796  ISNV  corresponding  to  clusters  of 
200  and  1500  atoms  respectively.  If,  however,  most  of  the  ISNV's  were  to  come  from  the 
island  boundaries,  one  would  observe  two  ISNV  sites.  One  would  come  from  the  ISNV 
next  to  one  vacancy  and  another  next  to  two  vacancies.  Only  the  ISNV  next  to  one  va¬ 
cancy  was  observed  [31].  It  is  therefore  argued  that  since  only  one  kind  of  ISNV  was  ob¬ 
served,  the  ISNV’s  were  within  the  islands.  Because  on  a  honeycomb  lattice  there  would 
be  three  ISNV's  near  to  one  vacancy,  the  vacancy  percentage  would  be  1/3  that  of  the 
ISVN  percentage  and  the  island  sise  would  be  much  bigger  than  the  above  limits.  On 
the  whole,  there  was  consistency  in  the  magnetic  measurements  between  the  samplee 
within  one  stage.  This  is  true  with  the  exception  of  several  samples  which  significantly 
deviated  from  the  norm  and  upon  further  examination  turned  out  to  be  FeCl*  interca¬ 
lated  compounds.  In  each  case,  more  than  one  sample  was  used  in  our  magnetic  mea¬ 
surements,  and  in  some  eases  as  many  as  5  or  5. 
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b)  Measurement  techniques 


A  standard  ac  bridge  technique  [6]  was  employed  to  probe  the  magnetic 
susceptibility^*)  of  the  system.  The  signal  was  picked-up  by  a  two-phase  lock-in  ana¬ 
lyzer  which  can  detect  signals  down  to  1  m  volt.  The  data  were  taken  at  several  frequen¬ 
cies  ranging  between  40  Hz  to  1000  Hz. 

A  computer-controlled  system,  via  A/D  and  D/A  converters,  was  used  to  operate 
the  apparatus  at  all  desirable  conditions.  The  temperatures  of  the  samples  above  2K 
were  measured  by  a  calibrated  silicon  diode  thermometer,  while  temperatures  below  2K 
were  determined  by  means  of  the  He4  vapor  pressure.  This  allowed  for  an  accuracy  of 
±lmK  at  2K.  When  greater  resolution  was  required,  a  carbon  resistor  was  employed. 

The  resistance  of  certain  carbon  resistors,  in  our  case  an  IRC  brand,  is  very  temperature 
sensitive  at  low  temperatures.  This  method  enabled  us  to  detect  temperature  changes  of 
less  than  O.lmK 

The  susceptibility  coils  were  always  kept  in  a  cryogenic  bath,  thus  changing  the  the 
temperature  of  the  sample  did  not  change  the  temperature  of  the  coils.  At  high  temper¬ 
atures  (nitrogen  to  room  temperature),  the  coils  were  immersed  in  liquid  nitrogen. 

An  ac  current  in  the  primary  circuit,  of  magnitude  4  mA  and  below,  was  used  to 
keep  the  amplitude  of  the  satcitfng  ac  field  below  0.1G,  thus  non-linear  susceptibility  ef¬ 
fects  were  excluded.  To  Investigate  the  susceptibility  as  a  function  of  the  magnetic  field, 
a  dc  magentk  field  la  the  range  of  (0-150G)  was  applied  to  the  samples.  Fig(l)  shows 
schematically  several  layers  of  GIC  with  the  intercalant  (striped  layers).  Generally  the 
probing  ac  and  the  external  dc  magnetic  fields  are  parallel  to  the  a-b  plane,  thus  the  bi¬ 
plane  "**1"^  susceptibility  component  was  measured  in  this  configuration.  Other  con¬ 
figurations  were  used  in  the  various  measurements  and  will  be  mentioned  as  appropriate. 
Our  sample  holder  was  so  that  the  orientation  of  the  sample  can  be  adjusted  to 


any  desired  configuration  relative  to  the  exciting  field  direction.  The  c-axia  susceptibility 
data  were  taken  when  the  a-b  plane  aligned  perpendicular  to  the  exciting  field.  Config¬ 
urations  where  the  dc  field  was  perpendicular  to  the  measuring  ac  field,  either  along  the 
c-axis  or  along  the  a-b  plane  were  also  used. 

Mechanical  vibrations  can  cause  serious  problems  in  this  kind  of  experiment,  thus 
careful  attention  has  been  paid  to  ensure  that  the  sample  was  firmly  attached  to  the 
sample  holder  and  in  a  rigid  configuration  with  the  susceptibility  coils.  A  special  com¬ 
puter  program  was  made  to  transfer  the  data  from  a  MACSYM-350,  ANALOG  DE¬ 
VICES,  Inc.,  computer,  which  monitors  the  apparatus  and  collects  the  data,  to  a  Digital 
Equipment  Corporation  VAXll  minicomputer  for  routine  analysis. 


RESULTS  AND  ANALYSIS 


a)  In- Plana  Susceptibility 

The  susceptibility  results  in  zero  external  dc  magnetic  field  of  stage  1  through  6, 
and  9  samples  are  presented  in  Fig.  (2).  In  this  configuration  the  measuring  ac  field  was 
perpendicular  to  the  c-axis,  thus  in-plane  susceptibility  data  are  shown  in  the  figure. 

The  susceptibility  is  giren  in  arbitrary  units,  however  it  is  normalized  for  the  amount 
of  iron  in  each  sample  and  at  a  frequency  of  40  Hz.  Also,  all  the  susceptibility  data  are 
subtracted  from  the  value  at  the  lowest  measured  temperature  1.1K.  At  low  frequencies 
and  away  from  the  transition,  (at  the  transition  the  relaxation  times  can  be  very  long, 
r  =  flx  10'*3s  [32])  the  arbitrary  unit  is  equal  to  1.0  emu/ mole  of  iron.  However,  the  sus¬ 
ceptibility  varies  with  the  measuring  frequency  and  thus  we  give  our  data  in  "arbitrary” 
units.  For  a  more  detailed  analysis  of  this  we  refer  the  reader  to  reference  32.  As  shown 
in  the  figure,  all  the  stages  have  low  temperature  susceptibility  maxima  at  about  the 
same  temperature  but  the  relative  size  of  the  peaks  is  different,  stage  5  has  the  largest 
peak. 

As  shown  in  the  insert  to  figure  2,  where  the  logarithm  of  the  peak  size,  Xm,  is  plot¬ 
ted  as  a  function  of  stage,  n,  the  peak  sise  for  stages  1  through  5  follows  an  exponential 
law.  More  exactly, 

Xm  *  0.028exp(n)  (1) 

for  1  <  it  £  5.  This  suggests  that  the  graphite  layers  shield  the  intercalant  layer  interac¬ 
tion  exponentially,  to  the  shielding  of  electromagnetic  radiation  from  the  interior 

of  a  conductor  due  to  the  skin  depth.  The  conducting  graphite  layers  between  the  inter- 
calaat  layers  shield  the  interplanar  intercalant  interaction,  and  thus  controll  the  dimen¬ 
sionality  of  the  magnetic  system.  In  stage  5  the  magnetic  intercalant  layers  are  22. 7  A 
apart.  The  sharp  increase  of  the  peak  sise  from  stage  1  to  stage  5  indicates  that  five 
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stages  of  the  graphite  layers  between  the  interealant  layers  are  sufficient  to  completely 
reduce  the  interlayer  interactions  and  produce  a  two-dimensional  magnetic  system.  As 
shown  in  Fig.(2)t  the  peak  sizes  started  to  decrease  for  samples  of  stages  higher  than 
stage  5. 

We  surmise  that  the  decrease  in  the  relative  peak  size  beyound  stage  5  is  related 
to  the  in-plane  density  of  the  magnetic  ions,  because  the  in-plane  interealant  density  of 
graphite-FeClj  compounds  decreases  as  the  stage  is  increased  [31].  As  shown  in  reference 
31,  the  in-plane  density  varies  within  stage  at  low  stages  from  100%  to  95%.  This  vari- 
tion  was  not  observed  in  our  susceptibility  measurements.  The  variation  of  the  in-plane 
density  can  be  attributed  to  the  existence  of  large  islands  [33,34]  in  the  interealant  com¬ 
pound,  a  decrease  in  the  number  of  iron  ions  in  these  islands  is  expected  when  the  stage 
index  is  increased.  Therefore  we  suggest  that  the  maximum  in  the  relative  peak  size  of 
the  susceptibility  versus  the  stage  index,  presented  in  Table  I,  is  due  to  the  competition 
between  the  dimensionality  of  the  system  and  the  interealant  in-plane  density.  Moreover 
it  is  interesting  to  note  that,  as  shown  in  Table  I,  the  transition  temperatures  also  vary 
with  the  stage  index  and  stage  4  has  the  lowest  transition  temperature.  This  indicates 
an  additional  aspect  of  the  competition  between  dimensionality  and  in  plane  density, 
since  otherwise  all  the  transitions  would  have  occured  at  the  same  temperature  for  all 
the  stages. 

The  anomaly  in  this  system  depends  strogiy  on  the  external  applied  dc 

magnetic  field.  The  susceptibility  data  versus  the  external  applied  dc  magnetic  field  for 
stage  6  and  stags  I.  as  examples,  am  presented  in  Fig.(3)  and  7ig.(4)  respectively.  Note 
that  the  scales  of  Fig. (3)  and  Fig.(4)  are  different  from  those  of  Fig.(2)  at  the  maximum. 
The  reason  is  that,  because  of  the  relaxation  times  [32]  the  susceptibility  varies  with  fre¬ 
quency.  Fig. (2)  represents  the  peak  sizes  at  397Hs  while  Fig.(3)  is  at  39Hs  and  Fig.(4) 
is  at  197Hs.  In  these  figures  the  exciting  field  is  perpendicular  to  the  c-axis,  thus  the 
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in-plane  susceptibilities  are  the  measured  data.  In  this  paper  we  have  adopted  the  con¬ 
vention  that  B  stands  for  the  enternal  exciting  ac  magnetic  field  and  H  stands  for  the 
external  applied  dc  magnetic  field.  In  this  case  H  is  applied  parallel  to  B.  As  shown  in 
the  figures,  the  peak  size  is  drastically  reduced  when  a  magnetic  field  of  the  order  of  5G 
is  applied  to  the  sample.  One  would  expect  that  in  three  dimensions  a  phase  transition 
would  be  quenched  by  a  magnetic  filed  at  least  in  the  order  of  k*T  or  l.AKG.  Therefore 
we  suspect  that  the  increase  in  the  number  of  graphite  layers  between  two  successive  in- 
tercalant  layers,  which  reduce  the  in  tore  slant  interaction,  contribute  to  the  strong  field 
dependence.  Spin  glasses  also  exhibit  strong  low  field  dependence  and  there  may  be  a 
relation  of  our  system  to  that  of  a  two-dimensional  spin  glass. 

In  Fig.(S)  and  Fig.  (A)  the  in-plane  susceptibility  at  constant  temperature,  normal¬ 
ized  to  the  aero  field  value,  is  shown  as  a  function  of  the  external  applied  dc  magnetic 
field  H  for  different  temperatures  for  stage  A  and  stage  3  respectively.  Here  H  is  again 
parallel  to  B.  The  remarkable  features  in  these  figures  are  the  low  field  local  minima  in 
the  susceptibility,  more  pronounced  for  stage  3,  which  steeply  decrease  on  the  left  shoul¬ 
der  near  the  transition  temperature.  This  rapid  turnover  of  the  susceptibility  at  a  very 
small  magnetic  fields  and  its  subsequent  increase  is  in  qualitative  agreement  with  the  re¬ 
sults  of  mean  field  calculations  for  two  dimensional  dipoles  on  a  honeycomb  lattice  [35]. 

In  these  mean  field  calculations,  we  considered  a  system  of  dipoles  on  a  two- 
dimensional  honeycomb  lattice.  Dipolar  interactions  are  inherently  anisotropic  since 
they  depend  on  both  in  •  and  ft  *  where  ft  denotes  the  magnetic  moment  of  the 
dipole  and  fy  is  the  distance  betweer  the  two  interacting  dipoles.  The  dipoles  were  con¬ 
strained  to  rotate  about  an  axis  normal  to  the  plane  and  thus  to  point  along  the  plane. 
That  system  has  a  continuously  degenerate  ground  state  whose  energy  is  -1.95K  per 
dipole  if  the  dipole  moment  and  separation  of  the  dipoles  is  that  of  the  iron  in  FeCU. 
There  is  a  phase  transition  in  zero  field  which  takes  place  at  IK  and  is  lowered  in  an 
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applied  field.  In  a  perfect  lattice,  the  magnetic  susceptibility  increases  with  an  applied 
magnetic  field  at  low  fields.  If,  however,  defects  are  taken  into  account,  the  suscepti¬ 
bility  has  a  minimum  with  the  application  of  an  external  field  at.  a  constant  tempera¬ 
ture.  Since  the  results  calculated  for  this  model  system  are  of  the  same  order  of  magni¬ 
tude  and  qualitatively  similar  to  those  measured  for  FeCfe  IGC,  we  assume  that  in-plane 
dipolar  interactions  play  a  significant  role  in  the  behavior  of  our  measured  systems.  As 
will  be  seen  below,  in  the  c-axis  measurements  at  low  temperatures,  the  spins  tend  to 
rotate  in  the  plane,  thus  approximating  the  configuration  of  the  model  system  above. 

It  is  of  interest  to  note  the  relative  large  increase  in  the  suaceptibiliy  on  the  high 
field  side  of  the  for  stage  3,  compared  with  that  for  stage  0,  where  the  suscep¬ 

tibility  exceeds  the  zero  field  value.  This  emphasises  the  still  existing  interplanar  inter¬ 
actions  in  the  lower  stages.  Stage  2  exhibits  similar  behaviour  to  that  of  stage  3  [18]. 

As  we  have  shown  the  interplanar  interactions  decrease  as  the  stages  increase,  thus  the 
in-plane  magnetic  phase  of  stage  3,  where  3-dimensional  effects  are  still  significant,  is  ef¬ 
fected  by  the  applied  field  in  a  way  as  to  reduce  the  interplanar  interaction  and  increase 
the  in-plane  interactions  of  the  system.  In  stage  6  the  interplanar  interactions  are  al¬ 
ready  quenched  by  the  graphite  layers  and  the  external  field  is  not  expected  to  increase 
the  in-plane significantly  in  this  low  magnetic  field. 

The  explanations  of  other  susceptibility  peaks  [24]  were  based  on  a  two-dimensional 
XY  model  [30],  however  the  existence  of  a  spin  glass  transition  was  also  suggested  [21]. 
The  model  calculations  mentioned  above  [36]  have  the  XY  symmetry.  Lundgren  et  al 
[37]  have  suggested  a  three  dimensional  spin  glass  model  to  relate  the  in-phase  and  out- 
of-phase  parts  of  the  magnetic  susceptibility.  We  have  tried  their  approach  to  the  analy¬ 
sis  of  the  magnetic  in  our  system  in  the  frequency  range  of  our  data  and  found 

no  quantitave  agreement  between  the  suggested  model  and  our  data.  Therefore  we  again 
conclude  that  the  magnetic  anomaly  in  this  system  is  two  dimensional  in  nature.  The 


spina  are  locally  locked  in  layered  planes  and  contibute  to  the  two-dimensional  anomaly. 
The  fact  that  the  spins  are  confined  to  the  the  a-b  plane  is  confirmed  by  our  e-axis  sus¬ 
ceptibility  measurements. 

b)  C-Axis  Susceptibility 

The  c-axis  susceptibility  was  measured  over  wide  temperature  range.  For  stage  3 
and  stage  6  we  did  not  observe  any  magnetic  anomaly  at  low  temperatures,  however  the 
susceptibility  showed  a  steep  decline  as  the  temperature  decreased  thus  indicating  that 
the  spins  are  locked  in  the  planes  as  temperature  is  lowered.  Appearently  the  Curie  law 
is  not  applicable  for  these  data  in  any  temperature  range. 

The  low  temperature  data  are  presented  in  Fig.(7).  The  figure  shows  the  c-axis  sus¬ 
ceptibility  of  stage  0  versus  the  temperature  at  different  external  dc  magnetic  fields.  As 
shown  in  the  figure,  the  susceptibility  in  this  direction  is  down  by  a  factor  of  100  from 
the  in-plane  value.  No  significant  anomaly  exists  in  this  direction.  Experimentally,  per¬ 
fect  orientation  of  the  sample  is  very  difficult  and  in  addition  the  graphite  planes  may 
be  out  of  parallel  by  as  much  as  a  few  degress,  thus  we  expect  a  small  contribution  from 
the  in-plane  anomaly  to  the  c-axis  data.  This  contribution  due  to  misorientation  appears 
as  a  small  peak  in  the  c-axis  susceptibility  data  near  the  transition  temperature. 

As  shown  in  Fig.(7),  the  applied  dc  magnetic  field  does  not  change  the  the  c-axis 
susceptibility  significantly  as  it  does  in  the  in-plane  ease  which  is  shown  in  Fig.  (3).  In 
this  high  stags  low  temperature  limit  the  spins  are  expected  to  be  confined  to  a  plane 
parallel  to  the  graphite  layers,  and  thus  are  free  to  realign  themselves  along  the  applied 
field  only  in  the  in-plane  directions.  Therefore,  in  this  small  magnetic  field  limit  and  for 
the  higher  stage  samples,  the  e-axis  susceptibility  Is  less  sensitive  to  the  applied  mag¬ 
netic  field  than  the  in-plane  data. 


Since  the  c-axis  susceptibility  data  at  low  temperature  were  continously  increasing 
as  a  function  of  temperature,  we  decided  to  extend  our  measurements  to  higher  tempera* 
tures  where  a  high  temperature  magnetic  susceptibility  maximum  might  exist.  As  shown 
in  Fig.(8),  the  stage  3  susceptibility  increases  monotonically  while  the  stage  6  suscep¬ 
tibility  has  a  shoulder  at  about  65K  and  a  small  &  temperature  of  about 

140K.  At  65K  the  e-axis  susceptibility  of  stage  8  is  a  factor  of  200  greater  than  that  at 
2K  and  below.  We  take  this  as  evidence  that  the  spins  are  locked  in  the  plane  at  low 
temperatures.  This  is  contrasted  with  stages  1  and  2  where  above  15K  the  magneti  sus¬ 
ceptibility  obeys  the  Curie* Weiss  law  both  along  the  a*b  and  along  the  e-direction. 

The  large  size  of  stage  6  c-axis  susceptibility  versus  the  temperatures  relative  to 
that  for  stage  3,  shown  in  Fig.  (8),  is  another  manifestation  of  the  two  dimensionality  at 
low  temperatures.  In  stage  6  the  spins  are  locked  in  the  plane  at  low  temperatures  and 
as  the  temperature  is  raised  the  spins  are  released  from  the  planar  orientation  and  are 
free  to  respond  to  the  field  in  the  e  direction.  In  stage  3,  the  interaction  between  planes 
hinders  the  response  of  the  spins  in  the  e  direction,  and  thus  the  susceptibility  is  smaller 
Stage  3  is  not  completely  two  dimensional,  and  interplane  interactions  play  a  role  in  this 
stage. 
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c)  Critical  Behavior 


A  maximum  in  the  magnetic  susceptibility  often  denotes  a  magnetic  phase  transi¬ 
tion  which,  if  it  is  a  second  order  transition,  can  be  analyzed  in  terms  of  critical  expo¬ 
nents  [38].  As  mentioned  above,  similar  transitions  were  anlyzed  in  terms  of  the  Koster- 
lit*  and  Thoules  XY  model  [30,40]  and  thus  we  decided  to  subject  our  data  to  a  similar 
analysis.  Before  this  analysis,  however,  the  measured  susceptibility  has  to  be  modified  so 
as  to  account  for  the  shape  of  the  sample  [41|. 

Because  of  the  high  susceptibility  at  the  maximum,  each  spin  does  not  see  the  exter¬ 
nally  applied  field,  but  a  field  which  is  modified  by  the  field  of  the  surrounding  dipoles. 
This  modification  or  shielding  depends  on  the  shape  of  the  sample.  In  what  follows  we 
designate  the  externally  measured  susceptibility  by  x««*.  and  the  actual  susceptibility  of 
the  spin  as  x*nt*  The  relation  between  the  two  susceptibilities  is 

Xint  *  X«st/(1  -  'Xfst)  (2) 

where  c  is  the  shape  factor.  The  shape  factor  can  be  exactly  calculated  only  for  ho- 
mogenious  substances  of  ellipsoidal  shape.  Neither  of  these  conditions  applies  to  our 
substance  since  it  is  an  inhomogenous  substance  in  the  form  of  thin  rectangular  plates. 
Because  of  that  on*  has  to  resort  to  experiment. 

If  one  assumes  that  at  the  maximum  this  Xwu  is  infinite  then 

•  *  1/X«at  (3) 

at  that  point.  A  finite  but  very  large  susceptibility  does  not  change  this  result  signifi¬ 
cantly. 

If  one  uses  the  value  obtained  for  t  from  equation  3  and  applies  it  to  our  data  by 
use  of  Eq.(2)  one  obtains  susceptibility  values  shown  in  Fig.(9).  This  figure  shows  the 


natural  logarithm  of  x.*i  a a  a  function  of  the  reduced  temperature  t  =  (T  -  Te)/Te 
where  Tc  is  the  temperature  at  the  susceptibility  maximum.  Although  this  procedure  in¬ 
troduces  an  infinity  where  only  rounded  maxima  exist,  it  is  justified  by  the  fact  that  our 
data  agree  well  with  equation  (4),  below,  over  a  wide  range  of  temperature.  Moreover, 
the  best  fit  was  found  with  x»»*  being  finite  and  of  the  order  of  400  rather  than  infinity. 

Figure  (10)  shows  the  natural  logarithm  of  the  susceptibility  as  a  function  of  the 
natural  logarithm  of  the  absolute  value  of  the  reduced  temperature  for  stage  0.  Simi¬ 
lar  data  reduction  was  performed  on  stages  4,  5  and  9  with  similar  results.  Stages  1  ,  2 
and  3  were  not  suited  for  this  analysis  because  of  interference  from  the  high  susceptibil¬ 
ity  due  to  3-dimensional  interactions  [18]  which  necessitated  the  subtraction  of  the  field 
sensitive  component  from  the  total  suseptibility,  with  an  accompanying  loss  in  accuracy. 
The  squares  are  for  T  >  Te  while  X  represent  data  for  2*  <  Te.  This  graph  suggests  the 
usual  power  law  behavior  of  the  susceptibility 

x«|(r-r«)/r.p  (4) 

with  7  =s  1.97  for  T  >  Tc  and  7  =  1.85  for  T  <  Tc.  The  error  in  7  is  ±0.1.  The  slopes 
of  the  drawn  lines  denote  7.  The  value  for  7  is  unusually  large.  For  a  three  dimensional 
system  7  has  the  value  between  1  and  1.25  while  the  calculated  value  for  the  two  dimen¬ 
sional  Ising  model  is  1.75  (see  reference  38  p.  47).  Our  values  appear  to  be  higher  than 
that  and  thus  it  appears  that  we  are  dealing  wih  a  somewhat  different  phenomenon.  Our 
data  were  also  fitted  to  an  equation  derived  by  Kosterlits  [40] 

X  «  ezp(fc“*)  (5) 

where  b  ia  a  constant.  Our  data  fit  Eq.(5)  only  in  a  very  limited  temperature  region. 

The  data  shown  in  Fig. (10)  are  well  within  the  accuracy  of  our  experiment  since  in|(T  - 
Te)/Te\  =  -7.5  represents  a  temperature  difference  of  one  millidegree. 
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Figures  (11)  and  (12)  show  the  magnetic  field  dependence  of  the  susceptibility.  Two 
configurations  were  used.  The  measuring  field  was  along  the  *»b  plane  while  the  applied 
dc  magnetic  field  was  normal  to  the  measuring  field  either  along  the  a-b,  in  plane  direc¬ 
tion  Ha,  or  along  the  c-axis  direction  H«.  Both  the  temperature  shift  of  the  maximum 
and  the  size  show  a  high  field  and  low  field  behavior.  The  horizontal  axis  is  calibrated  in 
terms  of  H«  where  H*  has  a  value  of  17G  for  He  and  7.5G  for  Ha. 

Figure  (11)  shows  the  temperature  at  which  the  peak  occurs,  in  terms  of  the  re¬ 
duced  temperature  (T-Te)/Te  as  a  function  of  the  applied  field.  Here  Te  is  the  trasition 
temperature  at  zero  field.  One  sees  that  the  maximum  is  shifted  to  higher  temperatures 
as  the  magnetic  field  is  applied.  Moreover,  two  straight  lines  can  be  drawn  through  the 
points.  One  at  low  field  with  a  steep  slope  and  one  at  high  field  with  a  shallower  slope. 
The  slope  at  low  field  is  a  factor  of  3  greater  than  that  at  high  field.  In  addition,  H« 
is  more  effective  in  shifting  the  temperature  of  the  maximum  than  H«  by  a  factor  of 
17/7.5,  the  ratio  of  H«  in  the  respective  directions.  A  similar  low  and  high  field  behav¬ 
ior  is  seen  if  one  plots  the  logarithm  of  the  susceptibility  peak  size  as  a  function  of  the 
applied  magnetic  field  normalized  to  H,,  as  shown  in  Fig.(12).  One  can  fit  the  data  in 
Fig.(12)  with  the  expression 

XwesW  *  / H»)  (6) 

with  Xmu  denoting  the  susceptibility  at  the  maximum.  For  H«  the  value  of  ^  =  1.3  for 
the  low  field  and  4  *  0.6  at  high  field,  while  for  H«  the  value  of  ^  a*  1  for  low  field  and 
4  »  0.4  for  high  field.  Similar  behavior  was  observed  in  other  GIC  samples  [24]. 

To  summarise,  our  data  obey  the  power  law  behavior  of  general  second  order  phase 
transitions,  with  an  exponent  7  greater  than  that  expected  for  a  two-dimensional  Ising 
model.  In  an  externally  applied  magnetic  field  the  susceptibility  maximum  exhibits  a 
distinct  low  field  and  high  field  behavior. 


CONCLUSION 


The  in-plane  magnetic  susceptibility  data  of  FeClj  GIC  indicate  that  all  the  stages 
of  this  intercaiant  compound  possess  a  low  temperature  magnetic  phase  transition  at 
temperature  between  1.7K  and  1.8K.  This  transition  obeys  the  power  law  dependence 
of  second  order  transitions.  The  data  suggest  that  the  transition  is  two-dimensional  in 
nature.  Similarities  in  the  behavior  of  the  transition  and  calculations  of  the  behavior  of 
a  two-dimensional  dipolar  system  on  a  honeycomb  lattice  suggest  that  dipolar  interac¬ 
tions  play  a  significant  role.  The  relative  size  of  the  anomaly  in  the  magnetic  suscepti¬ 
bility  is  maximum  for  stage  5.  If  the  size  of  the  low  temperature  susceptibility  maximum 
is  taken  as  an  indication  of  two-dimensionality,  then  the  shielding  of  the  magnetic  in¬ 
teraction  by  graphite  layers  is  exponential.  The  variation  in  the  size  of  the  *n«rfm«tn 
is  described  as  a  competition  between  the  staging  mechanism  which  contributes  to  the 
two-dimensionality  and  the  in-plane  density. 

As  shown  in  GIC,  the  graphite  layers  between  the  intercaiant  layers  screen  the  in¬ 
tercaiant  interplane  interactions  in  a  way  as  to  provide  a  two-dimensional  magnetic  sys¬ 
tem.  The  low  dimensional  nature  of  this  phase  trasition  is  supported  by  the  dependence 
of  the  susceptibility  peaks  on  the  magnetic  field,  the  application  of  a  small  dc  magnetic 
field  in  the  order  of  10G  smears  out  the  peak  completly,  whereas  one  would  expect  this 
to  happen  at  fields  of  the  order  of  l.OkG  In  throe  dimensions.  The  c-axis  susceptibility 
does  not  show  any  anomaly  at  low  temperature.  Moreover,  the  susceptibility  in  this  di¬ 
rection  is  not  sanative  to  the  small  dc  external  magnetic  fields  which  implies  the  confine¬ 
ment  of  the  intercaiant  magnetic  spins  to  the  planes  parallel  to  the  graphite  layers. 
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Table  I 


Stage  index  versus  the  relative  peak  size  and  the  peak 
temperature  of  the  susceptibility  for  FeCl3  IGC. 


Stage 

Peak  Size 

Temperature 

1 

0.0846 

1.793 

2 

0.2286 

1.745 

3 

0.4381 

1.746 

4 

1.4381 

1.710 

5 

5.1877 

1.735 

6 

1.8320 

1.761 

1  9 

1.4706 

1.768 
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FIGURE  CAPTIONS 


Fig.  1.  Layered  structure  for  stage  2  FeCl3-GIC,  the  striped  layers  are  the  intercalant. 
Fig.  2.  The  in-plane  magnetic  susceptibility  vs  temperature  for  different  stages  of  FeCb- 
GIC  near  the  transition  temperature.  The  insert  shows  the  logarithm  of  the  peak  height 
(Xm)  es  a  function  of  the  stage  index  (n)  to  emphasise  the  exponential  behavior  for 
1  <  n  <  5. 

Fig.  3.  The  in-plane  magnetic  susceptibility  vs  temperature  for  stage  6  FeClyGIC.  H 
denotes  the  applied  dc  magnetic  fields  in  Gauss  in  the  a-b  direction. 

Fig.  4.  The  in-plane  magnetic  susceptibility  vs  temperature  for  stage  3  FeClyGIC.  H 
denotes  the  applied  dc  magnetic  fields  in  Gauss  in  the  a-b  direction. 

Fig.  5.  The  in-plane  magnetic  susceptibility  vs  applied  magnetic  field  for  stage  6  FeCly 
GIC  near  the  transition  temperature. 

Fig.  6.  The  in-plane  magnetic  susceptibility  vs  applied  magnetic  field  for  stage  3  FeCl3- 
GIC  near  the  transition  temperature. 

Fig.  7.  The  low  temperature  c-axis  magnetic  susceptibility  vs  temperature  near  the 
transition  temperature  for  stage  A  FeClyGIC.  H  denotes  the  applied  dc  magnetic  fields 
in  Gauss. 

Fig.  8.  The  c-axis  susceptibility  vs  temperature  for  stage  3  and  stage  0  FeClyGIC  ex¬ 
tended  to  high  temperatures. 

Fig.  0.  Natural  logarithm  of  the  shape-corrected  in-plane  susceptibility  at  zero  dc  mag¬ 
netic  field  vs  the  reduced  temperature  for  stage  8  FeClygraphite  intercalation  com¬ 
pound. 

Fig.  10.  Natural  logarithm  of  the  in-plane  susceptibility  at  zero  dc  magnetic  field  vs 
natural  logarithm  of  the  reduced  temperature  showing  the  universal  power  law  behav¬ 
ior.  The  squares  denote  data  for  T>Te  and  the  x  denote  data  for  T<Te. 

Fig.  11.  Temperature  dependence  of  the  ir-plane  susceptibility  maximum  vs  the  applied 
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de  magnetic  field  for  stage  6  FeClj-GIC.  The  squares  are  for  H  along  the  c-axis  (He)  and 
the  x  are  for  H  along  the  a-b  plane  (Ha). 

Fig.  12.  Natural  logarithm  of  the  in-plane  susceptibility  maximum  vs  the  applied  de 
magnetic  filed  for  stage  6  FeClrGIC.  The  squares  are  for  H  along  the  c-axis  (He)  and 
the  x  are  for  H  along  the  a-b  plane  (H«). 
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Observation  of  spin  glass  state  in  FeCI3:  intercalated 
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Abstract.  FeClj  intercalated  graphite  compounds  exhibit  a  sharp  maximum  in  magnetic 
susceptibility  at  1.7  K.  By  using  Mdssbauer  spectroscopy  as  a  characterisation  tool,  we  show 
that  Hus  peak  is  dramatically  enhanced  as  the  number  of  iron  vacancies  in  the  intercalant 
lattice  is  increased.  These  results  strongly  suggest  that  the  sharp  susceptibility  maximum  is 
due  to  a  spin  glass  transitioo  where  the  iron  vacancies  inhibit  long-range  magnetic  order, 
thus  providing  a  new  series  of  compounds  to  help  elucidate  the  many  unknown  properties 
of  the  spin  glass  sum. 


The  onset  of  a  spin  glass  state  is  frequently  signalled  by  a  rather  sharp  peak  in  the 
susceptibility  at  low  magnetic  fields  (Mydosh  1981,  Maletta  1980).  Although  this  state 
from  a  theoretical  point  of  view  is  still  a  puzzle,  one  frequently  invokes  strong  magnetic 
short-range  order  as  a  qualitative  description  of  the  sharp  peak.  Chemical  short-range 
order  would  strongly  affect  such  magnetic  short-range  order  (Binder  1982).  In  this 
Letter  we  use  a  very  recently  developed  characterisation  technique  for  graphite — FeCI  j 
compounds.  Mdssbauer  spectroscopy,  to  analyte  the  microscopic  structure  of  each 
sample.  In  particular,  we  find  sites  which  are  nearest  neighbours  to  iron  vacancies  whose 
numbers  vary  from  sample  to  sample.  By  comparing  samples  having  x-ray  diffractograms 
which  are  almost  identical,  we  find  that  the  amplitude  of  the  sharp  susceptibility  peak 
shows  a  dramatic  increase  as  the  number  of  sites  which  are  nearest  neighbours  to  iron 
vacancies  is  only  slightly  increased,  as  if  approaching  a  percolation  threshold.  In  addi¬ 
tion,  we  provide  a  new  series  of  compounds,  the  magnetic  graphite  intercalation  com¬ 
pounds  on  which  the  spin  glass  state  can  be  studied,  which  has  quite  different  physical 
properties  from  the  usual  metallic  alloys. 

Previous  experiments  by  Elahy  et  d  (1982)  have  shown  that  the  existence  of  a  sharp 
peak  in  susceptibility  is  a  rather  universal  property  among  many  magnetic  compounds 
which  may  be  intercalated  into  graphite.  However,  these  experiments  have  been  inter¬ 
preted  as  being  due  to  the  existence  of  a  two-dimensional  XY  model  transition  and  have 
been  compared  to  the  theory  of  Kosteriitz  and  Thouless  (Suzuki  and  Ikeda  1981).  The 
sharp  peak  in  susceptibility  wit  proposed  by  Elahy  et  d  (1982)  as  being  the  universal 
magnetic  signature  for  magnetic  intercalants  of  stage  »  2.  In  this  work  we  show  that,  in 
the  limit  of  no  iron  vacancies,  the  susceptibility  curves  tend  to  an  antiferromagnetic  peak 
for  stage  2  FeGi  intercalants.  Even  for  samples  exhibiting  a  large  sharp  peak  in  the 
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susceptibility  curves  the  bulk  of  the  iron  atoms  still  enter  into  a  long-range  magnetically 
ordered  state. 

The  samples  studied  were  prepared  by  intercalating  anhydrous  FeClj  into  highly 
ordered  pyrolytic  graphite  using  the  two-zone  vapour  transport  method.  Employing  this 
growth  process,  the  graphite  temperature  is  maintained  at  350  °C  while  the  temperature 
of  the  ferric  chloride  is  varied  to  achieve  different  stages.  An  overpressure  of  about  250 
Torr  of  chlorine  gas  was  used  to  encourage  staging.  The  intercalated  compounds  were 
characterised  for  identity  and  uniformity  of  staging  using  (000  x-ray  diffraction  both 
before  and  after  our  susceptibility  measurements.  All  samples  were  essentially  single 
staged,  showing  only  very  weak  secondary  diffraction  peaks  due  to  other  stages. 

The  molecular  structure  of  grephite-FeCI)  has  been  determined  by  x-ray  and  electron 
diffraction  studies  and  by  electron  microscopy  studies  (Cowley  and  Ibers  1956,  Evans 
and  Thomas  1975).  These  studies  have  conducted  that  the  molecular  structure  of 
graphite-FeClj  is  essentially  unchanged  from  that  of  its  two  constituents,  as  in  pristine 
anhydrous  FeClj,  the  iron  ions  form  hexagonal  planes  where  each  iron  atom  is  sur¬ 
rounded  by  an  octahedron  of  chlorine  ions. 

The  susceptibility  measurements  were  made  by  a  standard  ac  bridge  method  (Max¬ 
well  1965)  operating  at  about  40  Hz.  The  amplitude  of  the  magnetic  field  at  the  sample 
was  kept  below  1  Oe.  The  susceptibility  coils  were  kept  at  a  constant  temperature  while 
the  temperature  of  the  samples  was  varied  inside  an  isolated  chamber.  The  temperatures 
above  4.2  K  were  measured  by  a  calibrated  carbon  resistor  and  a  calibrated  metal 
resistance  thermometer  (Millman  et  at  1982).  Below  4.2  K  the  sample  was  in  direct 
contact  with  the  helium  bath  and  temperatures  were  determined  from  the  vapour 
pressure  of  the  helium  bath. 

Mdssbauer  spectra  were  taken  in  transmission  geometry  with  the  y  ray  direction 
parallel  to  the  graphite  c  axis.  Temperatures  were  controlled  to  within  0. 1 K  and  were 
measured  by  a  calibrated  diode. 

In  figure  1  we  show  the  low-temperature  portion  of  our  susceptibility  curves  with  the 
measuring  field  parallel  to  the  basal  plane  for  three  stage  2  samples.  Although  the 
samples  have  almost  identical  x-ray  diffractograms,  the  susceptibility  curves  show  some 
quite  dramatic  differences.  In  particular,  the  sharp  peak  present  at  1.7  K  has  a  much 
different  amplitude  for  the  three  samples.  We  will  label  the  samples  in  order  of  increasing 
peak  amplitude  as  samples  1  to  3.  The  susceptibility  is  given  in  arbitrary  units  but  each 
curve  has  been  normalised  to  the  amount  of  iron  in  each  sample.  The  relative  amount 
of  iron  in  each  sample  was  calculated  by  comparing  the  areas  under  the  Mdssbauer 
curves  of  each  of  the  samples  at  the  same  temperature.  The  only  error  arising  in  this 
calculation  is  that  the  background  radiation  is  assumed  to  be  sample-independent. 
However,  since  the  same  Mdssbauer  apparatus  was  used  for  this  measurement  we  can 
assign  at  most  a  5%  error  bar  to  possible  slightly  different  sample  alignments.  We  will 
show  later  by  comparing  the  backgrounds  of  the  susceptibility  curves  that  this  method 
of  calculating  the  relative  amount  of  iron  is  quite  sensitive .  Although  one  usually  assumes 
that  the  stoichiometry  Q.*,  FeClj  (where  n  is  the  stage)  is  valid,  we  have  found  that  this 
method  may  be  quite  inaccurate,  due  either  to  formation  of  intercalant  islands,  a  gradient 
in  concentration  through  the  bulk  sample,  or  a  loss  of  some  of  the  intercalant  over  long 
periods  of  time  (Millman  1982a). 

The  Mdssbauer  spectra  of  these  samples  at  10.0  K  are  shown  in  figure  2.  Notice  the 
small  peaks  which  are  superimposed  on  both  sides  of  the  main  peak,  as  shown  in  the 
figure.  Although  we  will  discuss  the  interpretation  of  this  site  in  more  detail  elsewhere 
(Corson  etal  1982,  Millman  1982b),  briefly  it  has  an  isomer  shift  indicative  of  an  Fe}* 
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Flpn  1.  Susceptibility  venal  tempera tur*  Cor  three  diflenat  *U§*  2  graph) te-FeCli  com- 
pounda.  The  susceptibility  if  plotted  tn  arbitrary  uniti  but  each  simple  haa  been  normalised 
for  the  relative  amount  of  iron  it  contate.  A,  sample  1;  □.sample  2;  O.  sample  3. 
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Figure  2.  Corresponding  MOssbauer  spectra  for  the  three  samples whore  susceptibility  curves 
are  shown  in  figure  1.  The  position  of  the  two  peaks  which  comprise  the  iron  sites  nearest 
neighbour  to  iron  vacancies  are  indicated  by  the  straight  lines.  Zero  velocity  is  measured 
relative  to  the  centre  of  gravity  of  an  iron  foil  spectrum  at  room  temperature. 


ion  but  a  quite  large  quadrupole  splitting  for  an  Fe3*  ion  of  1.3  mm  s'1.  The  relative 
ratio  of  peak  intensities  is  3:3,  which  indicates  that  the  major  axis  of  the  electric  field 
gradient  tensor  is  pointing  in  the  basal  plane.  The  large  quadrupole  splitting  along  with 
the  direction  of  the  major  axis  of  the  electric  held  gradient  tensor  indicates  that  these 
sites  are  nearest  neighbours  to  iron  vacancies  in  the  intercalant  planes.  A  model  for 
iron  vacancies  in  graphite  FeCb  has  already  been  proposed  by  Wertheim  et  al  (1980). 
We  will  show  later  that  the  susceptibility  curves  are  also  consistent  with  the  existence 
of  an  in-plane  anisotropy  axis.  The  main  result  of  these  Mdssbauer  spectra  is  that  as 
the  number  of  sites  which  are  nearest  neighbours  to  iron  vacancies  is  increased  the 
sharp  peak  in  the  susceptibility  correspondingly  increases.  By  comparing  the  relative 
areas  of  the  Mdssbauer  spectra  our  least-squares  fits  to  the  spectrum  show  that 
7  ±  1%,  9  ±  1%  and  11  ±  1%  of  the  total  iron  sites  are  nearest  neighbours  to  iron 
vacancies. 
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The  height  of  the  sharp  susceptibility  peak  does  not  seem  to  depend  linearly  on  the 
number  of  sites  nearest  neighbour  to  iron  vacancies.  In  fact,  the  relative  magnitude  of 
this  peak  increases  for  the  three  samples  in  the  sequence  1,  l.S,  12.S  as  the  number  of 
sites  nearest  neighbour  to  vacancies  increases  as  7. 9, 11 .  This  nonlinear  dependence  on 
vacancies  may  indicate  that  a  percolation  threshold  is  being  approached.  We  also  find 
that  the  susceptibility  curves  show  very  anisotropic  behaviour.  If  we  measure  the  sus¬ 
ceptibility  perpendicular  to  the  basal  plane  we  find  only  a  very  tiny  contribution  of  the 
sharp  peak  at  1.7  K,  which  can  be  completely  explained  as  arising  from  a  parallel 
susceptibility  component  due  to  the  small  crystalline  misorientation  known  to  exist  in 
these  compounds.  Therefore,  the  susceptibility  curves  corroborate  the  Mossbauer 
measurements,  which  indicates  a  large  in-plane  anisotropy  axis. 

We  remark  that  even  for  the  sample  exhibiting  the  largest  peak  in  the  susceptibility, 
the  bulk  of  the  spins  still  exhibit  long-range  magnetic  order  as  we  have  reported  elsewhere 
(Millman  et  at  1982).  Since  the  sharp  peak  in  susceptibility  is  shown  to  be  quenched  by 
oc  fields  as  small  as  20  Oe  (Zimmerman  et  al  1981,  Eiahy  et  al  1982)  it  is  quite  easy  to 
deconvolve  the  two  contributions  to  the  susceptibility.  We  show  in  figure  3  the  curve  of 
one  of  the  samples,  sample  2,  for  which  the  sharp  peak  has  been  quenched  by  a  DC  field 
of  80  Oe,  along  with  that  with  no  external  DC  field.  We  note  that  the  backgrounds  of  the 
curves  shown  in  figure  1  are  almost  identical  in  the  paramagnetic  regime  (above  2.3  K) 
for  samples  1  and  2  and  only  slightly  smaller  for  sample  3.  The  slightly  smaller  back¬ 
ground  for  sample  3  can  be  explained  as  arising  from  the  increased  number  of  spins  for 
this  sample  which  are  contained  in  the  spin  glass  peak.  The  consistency  of  the  back- 


Flgas«3.  Susceptibility  versus  temperature  over  a  wide  rang*  of  temperature  tor  sample  2 
ta  figures  1  and  2.  The  ttne  rcptaaaata  the  susceptibility  curve  in  aero  external  oc  field  while 
the  trianglea  represent  the  susceptibility  cum  at  80  Oe,  which  has  been  applied  to  quench 
the  gain  glam  peak.  Notice  how  tmafl  the  sharp  epin  glass  peak  is  compared  to  the  rest  of  the 
susceptibility  cam  even  when  drawn  on  a  logarithmic  scale.  The  susceptibility  units  are  the 
tame  as  thoee  in  figure  1. 
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grounds  indicates  that  the  normalisation  of  the  susceptibility  curves  by  Mossbauer 
spectroscopy  is  quite  accurate. 

These  results  indicate  that  the  sharp  peak  observed  in  the  susceptibility  is  not  due  to 
the  two-dimensionality  of  the  compound  but  rather  to  vacancies  inhibiting  the  long- 
range  order  of  some  of  the  spins.  Anhydrous  FeCI}  is  known  to  contain  competing 
interactions  since  its  spin  structure  has  been  shown  to  be  a  spiral  through  neutron 
diffraction  (Cable  et  at  1962).  Therefore,  icreasing  the  number  of  vacancies  is  akin  to 
putting  non-magnetic  impurities  into  a  sys.em  containing  competing  interactions  which 
for  sufficient  non-magnetic  impurities  will  cause  the  system  to  go  over  to  the  spin  glass 
state  (Mydosh  1981,  Maletta  1980).  Although  these  results  have  only  been  shown  for 
the  FeG]  intercalant,  the  existence  of  similar  sharp  peaks  in  other  magnetic  intercalants 
suggests  that  intercalant  vacancies  also  occur  in  these  compounds.  In  fact,  small  inter¬ 
calant  islands  of  only  about  100  A  in  diameter  have  been  observed  in  graphite-NiOi 
(Flandrois  etal  1981).  It  seems  rather  likely  that  if  such  small  islands  are  occurring,  then 
vacancies  would  also  exist  in  these  islands. 

In  conclusion,  by  using  Mdssbauer  spectroscopy  as  a  tool  for  calibrating  the  micro¬ 
scopic  structure  of  FeClj  intercalated  graphite  compounds  we  have  found  a  strong 
dependence  in  the  magnetic  susceptibility  curves  with  the  amount  of  sites  which  are 
nearest  neighbours  to  iron  vacancies.  The  sensitivity  of  Mdssbauer  spectroscopy  is  even 
more  surprising  when  we  consider  that,  assuming  uncorrelated  vacancies,  the  threefold 
honeycomb  lattice  implies  that  only  4%  of  the  local  number  of  intercalant  atoms  are 
vacant,  even  for  sample  3.  These  results  may  also  help  in  understanding  the  spin  glass 
state  better,  since  they  introduce  a  new  series  of  compounds  for  which  it  can  be  studied. 

We  are  happy  to  acknowledge  many  helpful  discussions  with  G  Dresselhaus,  M  S 
Dresselhaus,  R  B  Frankel,  G  Kirczenow,  G  Papaefthymiou  and  S  Redner.  We  would 
also  like  to  thank  R  B  Frankel  and  G  Papaefthymiou  for  making  their  Mdssbauer 
spectrometer  available  to  us.  This  work  was  supported  by  AFOSR  subcontract 
F49620-81-C-0006  and  by  the  Francis  Bitter  National  Magnet  Laboratory  which  is 
supported  by  the  NSF. 
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The  ic  susceptibility  bridge  technique  is  applied  to  the  measurement  of  a  long  spin-lattice  relaxa¬ 
tion  time  with  the  magnetic  anomaly  which  hai  been  reported  for  FeOj  graphite  inter¬ 

calation  compounds  at  1.73  K.  Investigations  of  this  relaxation  mechanism  show  that  the  specific 
heat  of  the  spin  system  at  low  temperatures  contributes  to  an  anomaly  in  the  relaxation  time  near 
the  temperature  of  the  magnetic  anomaly.  Our  results  indicate  that  the  graphite  layers  between  the 
magnetic  intsrealant  layers  act  as  nonmagnetic  spacing  layers  allowing  the  c-axis  spin-spin  coupling 
to  be  varied. 


INTRODUCTION 


In  graphite  intercalation  compounds  (GlCs)  the  electri¬ 
cal  resistivity  is  one  of  the  properties  most  drastically 
changed  by  intercalation  of  both  acceptor  and  donor  ma¬ 
terials,1'2  it  is  also  one  which  holds  the  greatest  promise 
for  technological  applications.1  If  one  wants  to  use  the 
conventional  four-probe  technique  to  measure  the  resis¬ 
tivity  of  these  materials,  the  high  anisotropy  of  GICs  is 
fspatfil  to  introduce  difficulties  in  the  measurements. 
As  a  result  of  these  difficulties,  experimentalists  had 
thought  of  alternative  means  to  measure  the  resistivity  for 
these  highly  anisotropic  compounds. 

One  of  the  widely- used  approaches  in  the  measurements 
of  the  electrical  resistivity  in  GICs  is  the  low-frequency 
eddy-current  technique.4**  The  basic  concept  of  this 
method  is  that  when  a  sample  is  inserted  into  an  ac  field, 
the  induced  currents  in  the  conductor  reduce  the  magnetic 
flux  which  penetrates  the  conductor.  If  the  signal  is 
detected  by  a  standard  susceptibility  bridge7  the  induced 
eddy  current,  which  is  proportional  to  the  resistivity  of 
the  sample,  is  related  to  the  out-of-phase  or  quadrature 
component  of  the  n  qnetic  susceptibility.* 

Because  relaxatioi  Tecta  which  indues  an  out-of-phase 
signal  are  also  associated  with  low-temperature  magnetic 
anomalies,  which  have  been  repotted  for  a  variety  of 
GICs,*’10  one  has  to  be  careful  in  relating  the  out-of¬ 
phase  component  of  the  susceptibility  to  the  resistivity  of 
the  system.  Relaxation  effects,  therefore,  which  might 
contribute  to  the  out-of-phase  signal  should  not  be  ig¬ 
nored  near  a  magnetic  anomaly. 

In  our  system,  near  the  temperature  of  the  magnetic 
anomaly,  a  spin-lattice  relaxation  mechanism  is  suggested 
as  the  contributor  to  the  quadrature  component  of  the 
magnetic  susceptibility.  When  a  spin  system  (magnetic 
ions)  is  disturbed  from  its  thermal  equilibrium  by  an 
external  source,  such  as  an  applied  magnetic  field,  the 
heat  developed  in  this  spin  system  has  to  be  given  to  the 
lattice.  We  assume  that  the  spins  are  in  thermal  equilibri¬ 
um  with  each  other  and  at  a  temperature  slightly  different 
from  the  lattice  temperature.  The  transfer  of  heat  be¬ 
tween  the  spin  system  and  the  lattice  is  associated  with  a 
relaxation  time  (r),  which  depends  on  the  thermal  con¬ 
duction  between  the  spins  and  the  lattice  and  the  specific 


heat  of  the  spin  system.  Relaxation  effects  for  GICs  in 
which  low-temperature  magnetic  phase  transitions  were 
reported  have  not  yet  been  explored.  As  shown  in  the  fol¬ 
lowing  sections,  investigation  of  such  effects  can  provide 
information  about  the  heat  capacity  of  the  system  and 
also  about  the  nature  of  the  mentioned  magnetic  anomaly 
in  these  compounds. 

In  this  work  we  have  measured  the  long  relaxation  time 
in  the  FeClj-GIC  system,  using  the  susceptibility  bridge 
twhniqu*  We  have  also  used  the  conventional  four-probe 
technique  to  investigate  the  electrical  resistivity  at  the 
low-temperature  anomaly  which  has  been  observed  in  the 
out-of-phase  component  of  FeClrintercalated  graphite.’ 
In  fact,  no  anomaly  in  the  resistivity  along  the  c  axis  was 
observed  at  any  temperature.  These  resistivity  results, 
which  will  be  shown  in  this  work  for  comparison  and 
whose  details  are  presented  elsewhere,  are  in  qualitative 
agreement  with  the  theory  reported  by  Sugihara.11  There¬ 
fore,  it  is  dear  that  the  anomaly  of  the  out-of-phase  com¬ 
ponent  of  the  susceptibility  of  FeClj-intercalated  graphite 
compounds  is  due  to  the  relaxation  effects  in  the  system. 
Although  we  have  concentrated  on  the  measurements  of 
the  relaxation  time  of  stage-6  FeClj  GIC,  because  of  the 
similarity  in  the  «Mg«— anomaly  between  the  various 
stages,  we  can  draw  analogies  from  samples  other  than 
that  stags.  Thus  we  have  chosen  stages  5  and  9  for  the 
comparison  between  the  maximum  in  the  out  of  phase 
susceptibility  and  the  c-axis  resistivity. 

The  crystal  structure  of  FeClj  is  a  repeated  sequence  of 
three  layers  of  hexagonalty  arranged  atoms,12  these  layers 
are  displaced  relative  to  each  other  so  that  every  iron 
atom  is  surrounded  by  an  octahedron  of  chlorine  atoms. 
Upon  intercalation  of  FeClj  into  the  graphite,  one  elec¬ 
tron  is  donate  by  the  graphite  host  for  every  four  iron 
atoms.12  It  is  not  yet  clear  what  site  these  donated  elec¬ 
trons  occupy,  room  temperature  Mossbauer  data141’  and 
Raman  spectroscopy1*  do  not  support  the  existence  of 
FeClj  as  a  result  of  the  acceptance  by  FeClj  of  the  donat¬ 
ed  electron.17  However,  there  is  a  discrepancy  in  the  low 
temperature  Mossbauer  data  for  FeClj.  Millman  and 


Kircxenow11  have  reported  the  existence  of  Fe 


ly 


ions  at 


temperatures  as  high  as  100  K  which  is  contradictory  to 
the  work  of  Ohhashi  et  oi.1’  More  theoretical  and  experi¬ 
mental  work  is  necessary  to  understand  the  intercalation 
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mechanism  and  the  associated  charge  transfer  in  this  sys¬ 
tem. 

Studies  of  the  magnetic  properties  of  low  stage  FeClj 
GIC,  as  well  as  other  systems,  have  shown  a  low- 
temperature  phase  transition  and  different  mechanisms 
were  proposed9' 10, 19  to  describe  the  origin  of  this  transi¬ 
tion.  Limited  transport  measurements  have  been  reported 
for  low-stage  FeClj  GIC.  In-plane  thermal  conductivity 
data20  show  a  dominant  electronic  contribution  at  low 
temperatures,  while  high-temperature  data  indicate  a  pho¬ 
non  contribution.  The  c-axis  thermal  conductivity,  how¬ 
ever,  is  driven  by  phonons  at  all  temperatures.  The  au¬ 
thors20  have  also  shown  that  the  anisotropy  of  the  thermal 
conductivity  is  much  smaller  than  that  of  the  electrical 
conductivity.  It  has  been  reported21  that  the  in-plane 
resistivity  data  of  stage- 1  -and  -2  FeClj  GIC  exhibit  a 
low-temperature  phase  transition  while  stage  3  or  4  did 
not  show  any  anomaly.  It  is  the  purpose  of  this  work, 
which  was  initiated  a  few  yean  ago  at  Boston  University, 
to  investigate  in  a  full  and  comprehensive  way  the  proper¬ 
ties  of  the  FeClj-GIC  system. 


EXPERIMENTAL 


ing  ac  field  below  0.1  G,  thus  nonlinear  suscepubi|jtv  t 


feets  were  excluded.  To  investigate  the  relaxation 


j 

function  of  the  magnetic  field,  an  external  dc  field  jn 
range  of  (0—50  G)  was  applied  to  the  samples.  The 
netic  dc  and  ac  fields  configurations  are  shown  in  Fig  { 

It  shows  schematically  several  layers  of  GIC's  with  the 
probing  ac  and  the  external  dc  magnetic  fields  parallel  to 
the  a-b  plane,  thus  the  induced  eddy  currents  would  be 
normal  to  that  plane  along  the  c  axis.  Mechanical  vibra. 
tions  can  cause  serious  problems  in  this  kind  of  experi. 
ment,  thus  careful  attention  has  been  paid  to  ensure  that 
the  sample  was  firmly  attached  to  the  sample  holder  and 
in  a  rigid  configuration  with  the  susceptibility  coils. 

The  four-probe  method  was  used  to  probe  the  out-of. 
plane  resistivity,  and  an  on-off  current  technique  enabled 
us  to  subtract  out  any  ambient  or  spurious  voltage  from 
that  created  by  the  measuring  current.  A  special  comput¬ 
er  program  was  made  to  transfer  the  data  from  a  Symbol¬ 
ics,  Inc.  MACSYM-350  (ANALOG  DEVICES)  comput- 
er,  which  monitors  the  apparatus  and  collects  the  data,  to 
a  Digital  Equipment  Corporation  VAX  11  minicomputer 
for  routine  analysis. 


The  FeClj-GIC  samples  were  prepared  using  a  standard 
two-zone  furnace  technique2  where  stage  index  was  con¬ 
trolled  by  the  temperature  difference  between  the  graphite 
host  [highly-oriented  pyrolytic  graphite  (HOPG)]  and  the 
FeClj  powder.  The  graphite  samples  were  in  the  form  of 
thin  rectangular  plates  of  dimensions  1.5x0. 5x0.1  cm2. 
Well-staged  samples  were  achieved  by  controlling  the 
pressure  of  Cla  gas  inside  the  intercalation  tube,  as  well  as 
the  partial  pressure  of  FeClj  through  rigid  temperature 
control.  After  intercalation,  the  samples  were  character¬ 
ized  for  identity  and  uniformity  of  staging  using  x-ray 
(001)  diffraction.  The  x-ray  diffractograms  were  also 
used  to  determine  the  c-axis  repeat  distance  /,  after  cy¬ 
cling  the  samples  from  room  to  liquid  helium  temperature 
and  showed  that  the  cycling  did  not  affect  this  staging 
distance. 

A  standard  ac  bridge  technique7  was  employed  to  probe 
the  signal  which  is  in-phase  with  die  exciting  field  (related 
to  the  relaxation  of  the  system  or  the  out-of-phase  suscep¬ 
tibility  X")  and  the  signal  which  is  out-of-phase  with  the 
exciting  field  (related  to  the  in-phase  magnetic  susceptibil¬ 
ity  of  the  system  X').  These  components  were  picked  up 
by  a  two-phase  lock-in  analyzer  which  can  detect  signals 
down  to  1  n  V.  The  data  were  taken  at  several  frequencies 
ranging  between  40—1000  Hz. 

A  computer-controlled  system,  via  analog-to-digital 
(A/D)  and  digital-to-analog  (D/A)  converters,  was  used 
to  operate  the  apparatus  at  all  desirable  conditions.  The 
temperatures  of  the  samples  above  2  K  were  measured  by 
a  calibrated  silicon  diode  thermometer,  while  tempera¬ 
tures  below  2  K  were  determined  by  means  of  the  He  va¬ 
por  pressure.  The  susceptibility  coils,  were  always  kept  in 
a  cryogenic  bath,  thus  changing  the  temperature  of  the 
sample  did  not  change  the  temperature  of  the  coils.  At 
high  temperatures  (room  to  nitrogen  temperature),  the 
coils  were  immersed  in  liquid  nitrogen. 

An  ary  current  in  the  primary  circuit  of  magnitude 
about  VmA  was  used  to  keep  the  amplitude  of  the  excit¬ 


RESULTS  AND  ANALYSIS 


The  most  striking  result  in  this  work  is  the  temperature 
dependence  of  the  out-of-phase  component  of  the  magnet¬ 
ic  susceptibility  which  exhibits  an  anomaly  in  the  form  of 
a  sharp  peak  at  temperatures  near  1.75+0.05  K  in  zero  dc 
magnetic  field.  This  anomalous  behavior  is  correlated 
with  the  same  anomaly  which  we  have  seen  in  the  in- 
phase  component  of  the  magnetic  susceptibility.'*  Figure 
2  shows  the  correlation  between  the  in-phase  and  the  out- 
of-phase  component  as  a  function  of  temperatures  for 
stage-6  FeClj  GIC  at  /  =  39.7  Hz  and  in  zero  dc  field. 
As  shown  in  the  figure,  both  the  in-phase  part  UT'  and 
the  out-of-phase  part  (X")  exhibit  anomalies  in  the  form 
of  sharp  peaks  with  the  maximum  in  X"  shifted  to  lower 
temperatures. 

The  samples  were  oriented  inside  the  coils  m  such  a 
way  as  to  probe  the  in-plane  susceptibility  and  the  out-of- 
plane  resistivity  as  explained  in  the  preceding  section.  To 
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STAGE-2 

FIG.  1.  Layered  structure  for  stage-2  FeCL  GIC.  S,.  and 
Bjc  are  ihe  applied  ac  and  dc  magnetic  fields  along  (he  b  a  si* 
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FIG.  2.  The  in-phase  and  out-of-phase  components  of  the 
magnetic  susceptibility  versus  temperatures  near  the  transition 
point  for  stage-6  FeClj  GIC. 


determine  whether  the  correlation  between  X'  and  X"  is 
due  to  a  resistivity  anomaly  or  due  to  the  long  relaxation 
time  in  the  system,  the  four-probe  resistance  measurement 
technique  was  engaged  simultaneously  with  the  suscepti¬ 
bility  measurements.  Figure  3  shows  the  out-of-phase 
components  X"  as  a  function  of  temperature  for  stages  3 
and  9,  and  Fig.  4  shows  the  out-of-plane  resistivity  for  the 
same  samples  and  in  the  same  temperature  range.  As 
shown  in  the  figures,  there  is  no  anomaly  in  the  electrical 
resistivity  measured  by  the  four-probe  method  and  the 
resistivity  basically  behaves  as  predicted  by  Sugihar*." 
Therefore,  we  can  conclude  that  the  low-temperature 
maximum  in  X"  cannot  be  ascribed  to  an  anomaly  in  the 
resistivity,  but  is  related  to  the  relaxation  of  the  spin  sys¬ 
tem. 

Based  on  the  above  results  and  the  theoretical  descrip¬ 
tions  of  Casimir  and  Du  Prt,22  we  have  made  a  detailed 
investigation  of  the  relaxation  mechanism  for  stage-6 
FeClj  OIC.  The  real  and  imaginary  parts  of  the  complex 
susceptibility  are  related  to  the  relaxation  time  of  the  sys¬ 
tem  and  the  frequency  by  the  following  equations:11 


FIO.  3.  The  out-of-phase  susceptibility  component  versus 
temperature  for  stage-3  and  -9  FeClj  QIC's  near  the  transition 
temperature. 


FIG.  4.  The  c-axis  resistivity  versus  temperature  in  the  vicin¬ 
ity  of  the  magnetic  anomaly  for  stage-5  and.  -9  FeClj  G1CY 

x'mTT+h?)+x°-a  ’ 

r-7lSp?<!-r°“-c"/c*’-  121 

where  Cm  and  C*  are  the  heat  capacities  at  constant  mag¬ 
netization  and  at  constant  field  respectively,  and  T0  is  the 
static  susceptibility,  r  is  the  relaxation  time  of  the  spin 
system,  and  a  is  the  angular  frequency  of  the  exciting 
field. 

The  out-of-phase  magnetic  susceptibility  component 
data  were  taken  as  function  of  frequency  and  at  tempera¬ 
tures  between  1.1  and  2.4  K.  Since  that  component  varied 
significantly  with  frequency  only  in  the  temperature  range 
between  1.63  and  1.83  K,  constant  temperature  cuts  were 
made  in  this  temperature  range  and  the  measured  values 
of  the  quadrature  component  of  the  magnetic  susceptibili¬ 
ty  at  constant  temperature  fit  to  the  function  given  in  Eq. 
(2).  We  used  a  nonlinear  least-squares  computer  fitting 
program  allowing  G  and  r  to  be  varied  to  achieve  the  best 


FIG.  5.  The  out-of-phase  component  of  the  magnetic  'u-.c 
tibtlity  versus  the  angular  frequency  for  stage-0  FeCI  t<K 
The  continuous  lines  are  the  theoretical  functions  and  the  "  n 
bolt  are  the  experimental  data. 
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fit.  Figure  S  shows  the  fits  of  the  measured  data  at  tern* 
peratures  1.68,  1.72,  and  1.78  K.  As  shown  in  the  figure, 
the  data  indicate  a  good  fit  to  the  Casimir— -Du  Pre  rela¬ 
tion  in  the  low-frequency  region  and  a  somewhat  worse  fit 
at  frequencies  near  1  kHz.  This  is  due  to  the  fact  that  in 
the  high-frequency  region,  it  is  difficult  to  decouple  the 
two  components  of  the  measured  susceptibility  from  each 
other  because  capacitive  and  skin  effects  impose  an  upper 
limit  on  the  sensitivity  of  the  bridge.  The  high-frequency 
points,  however,  were  given  a  lower  weight. 

The  fitting  parameters  r  and  G  were  found  to  be  tem¬ 
perature  dependent  and,  in  fact,  both  of  them  exhibit 
maxima  near  the  temperature  of  the  susceptibility  anoma¬ 
ly.  In  Fig.  6  the  values  of  these  parameters  are  represent¬ 
ed  as  a  function  of  the  temperature  and,  as  shown  in  the 
figure,  the  peak  of  the  relaxation  time  is  shifted  towards 
lower  temperature  from  that  in  G.  It  is  interesting  to  note 
that  the  out-of-phase  component  of  the  susceptibility,  as 
shown  in  Fig.  2,  shifts  towards  lower  temperatures  as  the 
relaxation  time  does,  while  the  in-phase  component  as 
well  as  the  fitting  parameter  G  are  shifted  towards  higher 
temperatures. 

The  parameter  G  can  be  also  written  as 
<7  •  <3> 

where  X,  and  X,  are  the  susceptibility  at  constant  tem¬ 
perature  and  entropy,  respectively.  In  the  low-frequency 
region  and  at  low  temperature,  one  can  assume  that  X,  is 
small  and  expect  G  to  have  a  maximum  which  coincides 
with  that  of  the  in-planc  susceptibility  which  is  propor¬ 
tional  to  Xq.  This  explains  the  correlation  between  the 
graphs  in  Fig.  2  and  Fig.  6  which  showed  a  shift  towards 
higher  temperatures  for  X'  and  G  relative  to  X"  and  r, 
respectively.  Measurements  at  low-frequency  provide  an 
isothermal  spin  system  so  that  the  spin-spin  relaxation 
can  be  ignored.  The  relaxation  time,  on  the  other  hand, 
depends  on  the  temperature  through  the  specific  heat  at 
constant  field  of  the  spin  system  and  the  beat  conduction 
between  the  spin  system  ami  the  lattice  system.  Thus,  in 
order  to  extract  useful  information,  careful  examination 
has  to  be  carried  out  in  this  case. 


FIG.  6.  The  fitting  parameters  G  and  r  versus  temperatures 
for  stage-6  FeClj  GIC.  The  continuous  lines  are  used  to  guide 
the  eye. 
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We  have  also  measured  the  spin-lattice  relaxation  tim 
in  stage-3  FeClj  GIC  and  found  that  the  relaxation 
in  stage  3  were  shorter  by  a  factor  of  4  compared  to  those 
of  stage  6.  According  to  the  Casimir-  Du  Pre  picture,  re 
taxation  occurs  when  heat  is  transferred  from  the  spin 
system,  in  our  case  FeClj  whose  temperature  changes  in  > 
sinusoidal  fashion,  to  a  lattice  which  is  at  a  constant  tent 
perature  and  serves  as  a  heat  reservoir: 

r»C/K ,  (4) 

where  C  is  the  spin  specific  heat  and  K  is  the  thermal 
conductivity.  If  heat  transfer  were  due  to  the  electrons  in 
the  bounding  layers,  the  retaxation  time  r  would  be  too 
short  for  us  to  measure  by  this  method.  It  was  also  shown 
that  the  spacing  between  adjacent  graphite  layers  is  nearly 
independent  of  staging.24  23 

Because  of  this,  and  the  fact  that  near  a  magnetic  tran¬ 
sition  the  specific  heat  is  related  to  the  magnetic  suscepti. 
bility  X  by2i 


d(XT) 

dT 


(5) 


we  can  attribute  the  difference  between  the  stage-3  and  -6 
relaxation  times  to  the  differences  in  their  respective 
specific  heats  at  the  anomaly.  Table  I  gives  the  relative 
magnitudes  of  the  susceptibility,  normalized  per  iron 
atom,  of  various  stages  of  FeClj  GIC  at  the  susceptibility 
maximum.  As  shown  in  Table  l,  stage  3  has  a  maximum 
whose  size  is  approximately  a  factor  of  4  smaller  than 
that  of  stage  6. 

Another  possible  mechanism  for  the  difference  in  the 
relaxation  times  between  the  various  stages  could  be  the 
intercalate  bounding  layers  of  graphite.  One  imagines 
that  the  Fe3*  spins  relax  to  the  FeClj  lattice  which  then 
relaxes  to  the  graphite  layers  mainly  through  the  ch tonne 
graphite  interaction  due  to  phonons.  An  enhancement  in 
the  low-frequency  phonon  spectrum  of  intercalated  gra¬ 
phite  over  that  of  HOPG  was  actually  calculated27  and  ob¬ 
served.2*  Since  high  stages  have  a  smaller  density  of  low- 
frequency  phonons  than  low  stages,  according  to  equation 
(4)  the  conductivity  K  in  the  low  stages  would  be 
enhanced  and  thus  again  lead  to  a  shorter  relaxation  lime. 
However,  the  changes  in  the  relaxation  times  can  be  ade¬ 
quately  accounted  for  by  the  relative  size  of  the  specific 
heats  without  considering  the  changes  in  the  thermal  con¬ 
ductivity.  Thus  we  conclude  that  the  change  in  thermal 
conductivity  between  different  stages  of  FeClj  is  small 
compared  to  the  specific  heat  effect.  * 

Therefore  the  maximum  in  the  relaxation  time  is  attn- 
buted  to  the  specific  heat  of  the  spin  system  which  is  con- 

TABLE  I.  Stage  index  versus  the  relative  peak  sue  of  (he  in- 
phase  magnetic  susceptibility  for  Fed)  QIC 


Stage 

Peak  sue 

1 

0.0846 

2 

0.2256 

3 

0.4311 

4 

1  4381 

J 

5.1877 

6 

1  8329 

4 

1.4706 
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siiHfnt  with  theoretical  calculations  based  on  a  model  of 
two-dimensional  interacting  dipoles.29  This  model  shows 
an  increase  in  the  specific  heat  of  the  system  near  the 
transition  point  which  is  similar  to  that  exhibited  by 
many  other  systems  near  their  critical  points  and  is  attri¬ 
buted  to  critical  fluctuations.  Mossbauer  studies20  of 
anhydrous  FeClj  describe  the  broadness  of  the  resonance 
tin**  near  the  critical  temperature  as  a  result  of  the  in¬ 
crease  in  fluctuations  of  the  internal  field  direction  which 
are  caused  by  the  spin  relaxation  in  the  system. 

If  the  quantity  or  is  eliminated  between  Eqs.  (1)  and 
(2),  then  another  representation  of  X'  and  X"  can  be  given 
by  the  following  equation: 

(X'-*o+t G)l+iX")l=iGl .  (6) 


T  (K) 


As  long  as  the  relaxation  is  controlled  by  one  relaxation 
time,  Eq.  (6)  predicts  that  when  X'  is  plotted  versus  X"  the 
points  should  lie  on  a  circle.  Figure  7  shows  the  experi¬ 
mental  data  of  X'  versus  X"  at  different  frequencies  at  the 
peak  temperatures  and  1.72  K,  which  is  about  30  mK 
below  the  peak  temperature. 

The  existence  of  the  spin  relaxation  at  the  critical  point 
for  both  anhydrous  FeClj  and  FeClj  GIC  supports  the 
correlation  between  the  magnetic  properties  of  these  two 
systems.  Thus,  qualitatively,  we  suggest  that  the  three- 
dimensional  anomaly  which  is  observed  at  about  8  K  in 
anhydrous  FeClj  is  the  source  of  the  low-temperature 
anomaly  in  FeClj  GIC.  The  staging  process  decreases  the 
interactions  between  the  intercalant  layers  in  such  a  way 
as  to  have  a  two-dimensional  interacting  spin  system. 


FIG.  8.  The  out-of-phase  component  of  the  magnetic  suscep¬ 
tibility  versus  temperatures  for  stage-6  FeClj  GIC.  The  0's  are 
the  applied  dc  magnetic  fields  in  Gauss. 

in  Fig.  8;  as  shown  in  the  figure  at  a  field  of  4.2  G  the 
anomaly  has  disappeared. 

This  anomaly,  if  it  is  a  three-dimensional  transition  be¬ 
tween  two  magnetic  phases,  is  expected  to  exist  in  fields  as 
high  as  1  kG;  however,  the  disappearance  in  such  a  small 
field  emphasizes  the  two-dimensional  nature  of  this  phase 
transition.  In  addition,  the  stage  dependence  of  this 
anomaly  shows  more  pronounced  peaks  for  higher  stage 
samples  in  both  components  of  the  susceptibility. 

CONCLUSION 


I 


Our  data  indicate  that  the  graphite  layers  between  the 
magnetic  intercalant  layers  act  as  nonmagnetic  spacing 
layers  allowing  the  c-axis  spin-spin  coupling  to  be  varied. 

The  application  of  a  small  dc  magnetic  field  of  the  or¬ 
der  of  5  G  smears  out  the  anomaly  in  the  susceptibility21 
and  thus  no  relaxation  was  observed  at  this  field.  The 
out-of-phase  component  (*")  as  a  function  of  temperature 
and  in  different  external  dc  magnetic  fields  is  represented 


We  have  shown  that  for  systems  which  possess  a  mag¬ 
netic  anomaly  the  out-of-phase  component  of  the  suscep¬ 
tibility  is  dominated  by  spin-lattice  relaxations,  thus  it 
would  not  be  reliable  to  relate  the  electrical  resistivity  of 
these  systems  to  the  out-of-phase  component.  Fitting  our 
data  to  a  temperature-dependent  relaxation  shows  that  the 
relaxation  time  itself  exhibits  an  anomaly  near  the  transi¬ 
tion  temperature  of  the  magnetic  anomaly.  On  a  qualita¬ 


tive  basis,  and  in  agreement  with  a  two-dimensional 
mean-field  calculation,21  we  have  found  that  the  relaxa¬ 
tion  anomaly  is  related  to  a  similar  one  in  the  specific 
heat  at  constant  field. 

In  agreement  with  repotted  thermal  conductivity  mea¬ 
surements20  and  phonon  density  calculations27  at  low  tem¬ 
peratures,  the  out-of-plane  phonon  vibrations  are  deter¬ 
mined  to  be  the  carriers  of  the  heat  from  the  spin  system 
to  the  graphite  lattice.  The  contribution  of  the  conduc¬ 
tion  electrons  is  not  significant  in  the  temperature  range 
over  which  our  measurements  were  carried  out.  Finally, 
we  have  shown  that  the  intercalant  FeClj  controls  the 
magnetic  properties  of  the  FeClj-GIC  system  and  be¬ 
comes  more  two  dimensional  as  the  stage  is  increased. 
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The  four  probe  technique  is  used  to  measure  the  c-axis  resistivity  of  various  stages 
of  FeCla  acceptor  graphite  intercalated  compounds.  Comprehensive  investigation  of  the 
resistivity  behavior  between  room  and  liquid  helium  temperatures  for  stages  2-5,  and  9 
is  presented.  The  experimental  results  are  compared  with  one  of  the  existing  theoretical 
models.  There  is  good  agreement  between  the  experimental  data  and  the  predictions 
of  one  of  the  theoretical  models.  The  results  of  this  work  reveal  that  the  conduction 
along  the  c-axis  in  acceptor  compounds  is  a  hopping  mechanism  assisted  by  phonons 
and  impurities.  In  lower  stage  compounds,  the  conduction  is  dominated  by  impurity- 
assited  hopping,  and  the  resistivity  shows  metallic  temperature  behavior  at  high  temper¬ 
atures  and  tends  to  saturate  at  low  temperatures.  For  higher  stage  compounds,  phonon- 
assisted  hopping  dominates  the  conduction.  The  resistivity  approaches  constant  values 
at  high  temperatures  and  exhibits  an  activated  behavior  which  is  proportional  to  T-2  at 
low  temperatures.  Stage  5  seems  to  be  an  intermediate  stage  where  the  conduction  has 
the  contribution  of  both  phonon  and  impurity  hopping  mechanisms. 
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INTRODUCTION 


In  graphite  intercalation  compounds  (GIC’s)j_the  electieal_conductivity  is  one  of  the 
most  promising  physical  quantities  for  which  real  applications  might  take  place  [1].  The 
in-plane  and  c-axis  conductivities  are  among  the  properties  most  drastically  changed  by 
intercalation  of  acceptor  and  donor  GIC’s  [2,3].  This  unexpected  change  of  the  conduc¬ 
tivity  attracts  numerous  theoreticians  and  experimentalists  to  investigate  these  GIC’s. 
Inspite  of  that,  there  is  a  lack  of  a  coprehensive  theory  and  experimental  results  which 
can  describe  the  electrical  conduction  in  these  materials.  The  high  anistropy  of  the  con¬ 
ductivity,  which  can  be  a  factor  of  10®  for  some  compounds  [4],  and  the  staging  property 
introduce  difficulties  to  the  current  theories  which  attempt  to  formulate  a  complete  and 
realistic  picture  about  the  conduction  in  these  compounds. 

The  conduction  in  GIG's  has  been  described  by  Sugihara  et  al  [5]  and  Sugihara  [6] 
in  terms  of  a  scattering  mechanism  caused  by  in-plane  and  out-of-plane  vibrations.  The 
interaction  between  phonons  and  these  out-of-plans  vibrations  is  the  contributor  to  the 
c-axis  conduction.  The  authors  of  reference  [6]  suggest  that  this  conduction  is  a  hop¬ 
ping  mechanism  assisted  by  phonons  and  impurities.  The  in-plane  conduction  is  not  yet 
based  on  any  satisfaction  theory,  however  recent  papers  [7,8]  concluded  that  a  strong 
charge  screening  mechanism  takes  place  in  the  e-direction  and  most  of  the  free  carriers 
are  confined  to  the  carbon  planes  adjacent  to  the  intercalated  layers.  It  is  the  purpose 
of  this  work  to  investigate  the  c-axis  conduction  in  the  FeCly  graphite  intercalated  sys¬ 
tem  and  provide  experimental  results  which  can  be  compared  with  existing  theoretical 
models. 


The  FeCls-GIC  systenfls  one  of  the  earliest  systems  to  have  been  investigated  by 
numerous  researchers.  However  the  understanding  the  low  temperature  magnetic  phase 
transition,  the  unexpected  behavior  of  the  transport  and  electronic  properties,  and  the 
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staging  mechanism  in  this  system,  as  well  as  other  GIC  systems,  still  are  open  questions. 
Met*  et  al  [0]  have  made  a  structural  investigation  using  x-ray  diffraction  data  and  con¬ 
cluded  that  the  stacking  of  layers  is  disordered  and  under  special  circumstances,  depend¬ 
ing  on  the  preparation  condition,  nearly  periodic  compounds  can  be  produced.  IR  and 
Raman  spectroscopy  of  FeCl3-graphite  [10]  indicate  in-plane  compression  within  the 
graphite  layers  and  most  of  the  effective  charge  is  in  the  graphite  bounding  layers,  with 
the  rest  of  the  charge  distributed  among  the  graphite  interior  layers.  Thermal  conduc¬ 
tivity  measurements  [11]  showed  that  the  electronic  thermal  conductivity  of  the  graphite 
bounding  layers  increases  upon  intercalation,  whereas  the  lattice  thermal  conductivity 
decreases.  Magnetic  susceptibility  measurements  [12,13]  on  FeCfe  indicate  the  existence 
of  a  low  temperature  magnetic  phase  transition. 

In  this  work,  the  c-axis  resistivity  measurements  of  the  FeCl3-graphite  system  us¬ 
ing  the  convetional  four-probe  technique  were  stimulated  by  the  investigation  of  the  low 
temperature  anomaly  of  the  magnetic  susceptibility  [13].  Contactless  bridge  measure¬ 
ments  had  been  made  on  various  stages  and  a  peak  in  the  susceptibility  was  observed  at 
around  1.7SK.  This  peak  was  also  noted  in  the  out-of-phase  component  of  the  suscepti- 
bilty  which  indicated  a  possible  contribution  from  the  resistivity  of  the  sample  [14].  For 
the  chosen  sample  orientation  in  the  measurements,  this  contribution  would  have  been 
primarily  from  the  c-axis  resistivity.  Thus  it  was  decided  to  investigate  the  resistivity  by 
means  of  independent  measurements  concurrently  with  the  contactless  bridge  method. 

Although  numerous  theoretical  and  experimental  investigations  of  the  in-plane  con¬ 
ductivity  have  been  reported,  less  attention  was  paid  to  the  c-axis  conductivity  of  GIC’s 
[15].  Theoretically,  very  limitted  work  has  been  reported,  Markiewics  [16]  has  proposed 
a  simple  band  model  for  calculating  the  c-axis  conductivity  in  a  series  of  GIC’s.  In  this 
model,  it  was  assumed  that  for  acceptor  compounds  the  bandwidth  is  small  enough  so 
one  can  treat  these  compounds  ss  two-dimensional  systems.  Although  there  is  a  small 


contribution  duo  to  band  conduction,  it  is  generally  overpowered  by  other  mechanisms. 
The  hopping  mechanism  suggested  by  Sugihara  [6]  excludes  the  contribution  of  the  band 
conduction  along  the  c-axis  in  acceptor  compounds.  Experimentally,  Ohts  et  al  [17] 
have  measured  the  temperature  dependence  of  the  c-axis  resistivity  for  stages  1  and  2 
of  ICl-graphite  compounds.  Their  data  exhibit  order-disorder  transition  at  room  temper¬ 
ature  and  anistropy  of  the  order  of  10s  for  stage-1  and  104  for  stage-2.  This  anistropy 
was  attributed  to  the  hopping  mechanism  proposed  by  Sugihara  [6].  The  investiga¬ 
tion  of  the  c-axis  resistivity  of  SbClg  by  Uher  and  Morelli  [18]  indicates  that  low  stage 
compounds  exhibit  metallic-like  conduction,  which  was  described  by  a  model  based  on 
defect-mediated  short-circuiting  channels  along  the  c-axis.  In  in  higher  stages,  however, 
a  metal-insulator  transition  was  observed  at  low  temperature  and  the  conduction  is  es¬ 
sentially  caused  by  a  thermally  activated  hopping  mechanism.  Measurements  of  the  c- 
axis  resistivity  for  stage  2  FeCls-GIC's  by  Issi  et  al  [10]  indicate  that  resistivity,  after 
subtracted  the  residual  part,  has  Ta  dependence  at  low  temperature  and  a  T  dependence 
at  high  temperature.  A  compendium  of  various  measurements  of  the  c-axis  conductivity 
for  several  GIC’s,  both  acceptor  and  donor,  was  presented  recently  by  McRae  et  al  [20]. 
It  shows  that  there  are  sosna  •‘nmmnn  characteristic  signatures  in  the  c-axis  conductivity 
in  both  acceptor  and  donor  GIC's.  Nevertheless,  those  characteristics  are  not  fully  un¬ 
derstood  and,  to  reach  a  comprehensive  picture  about  this  conduction,  more  theoretical 
and  experimental  works  on  several  different  systems  should  be  done. 
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The  FeCU-GIC  —  mples  were  prepared  us  Inf  a  standard  two- sons  furnace  technique 
[l]  where  state  index  was  controlled  by  the  temperature  difference  between  the  graphite 
hoet  highly  ordered  pyrolytic  graphite  (HOPO)  and  the  FeCls  powder.  The  samplee 
were  in  the  form  of  thin  rectangular  plates  of  dimensions  1.5x0. 5x0.1  cm9.  Well-staged 
samples  were  achieved  by  controlling  the  pressure  of  Clj  gas  Inside  the  Intercalation 
tube,  as  well  as  the  partial  pressure  of  FeCU  through  rigid  temperature  controll.  After 
intercalation,  the  samples  were  characterised  for  identity  and  uniformity  of  staging  us¬ 
ing  x-ray  (001)  diffraction.  Only  single  stage,  well  staged  samples  were  used  in  the  mag¬ 
netic  measurements.  The  x-ray  diffractograms  were  also  used  to  determine  the  c-axis 
repeat  distance  I«  after  cycling  the  samples  from  room  to  liquid  helium  temperature,  and 
showed  that  the  cycling  did  not  affect  this  staging  distance. 

Moet  of  the  samples  measured  were  characterised  by  means  of  the  Mdesbauer  effect, 
details  are  reported  in  references  (13,21).  In  higher  stage  samples,  where  stage-disorder 
is  expected  (22),  the  Hendricks- Teller  (23),  and  Mats  and  Hohlwein  [0]  analysis  technique 
were  used  to  calculate  the  intensity,  width  and  location  of  the  x-ray  reflections.  We  find 
that  our  experimental  x-ray  data  on  the  stage  0  samples  reported  in  this  article  are  in 
good  agreement  with  that  calculated  for  the  pure  and  well-staged  stage  0.  Although  a 
small  admixture  of  stage  10  can  not  be  excluded,  this  stage  serves  as  an  example  of  a 
high  stage  sample. 

Samplee,  for  resistivity  measurements,  of  approximate  dimensions  3  mm  by  5x5 
mm9  were  peeled  out  from  the  intercalated  samples  using  scotch  tape.  These  dimensions 
were  chosen  primarily  for  convenience  when  mounting  the  samples  in  our  sample  holder, 
whereas  optimal  signal  would  have  been  obtained  by  minimising  the  in-plane  area  while 
comfortable  lead  separation  on  the  faces  and  maximising  the  c-axis  thick- 
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ness.  Because  of  the  high  in-plane  conductivity  relative  to  that  of  the  c-axis,  the  non- 
uniform  current  injection  would  not  be  a  problem.  The  leads  were  attached  to  the  sam¬ 
ple  in  the  conventional  four-prope  configuration.  Contacts  were  made  on  the  sample  us¬ 
ing  GC  conductive  silver  print  (GC  ELECTRONICS).  The  samples  were  inserted  in  the 
axial  slot  of  a  phenolic  rod  which  was  then  wrapped  with  mylar  tape  to  insure  that  the 
sample  did  not  move  in  the  slot  and  that  subsequent  stress  to  the  leads  was  minimized. 
The  leads  were  attached  to  the  main  leads  leading  from  the  cryostat  by  wrapping  and 
then  glueing  with  silver  print. 

The  sample  holder  was  mounted  in  a  vacuum  case  which  was  surrounded  by  a  he¬ 
lium  dewar  with  an  outer  nitrogen  dewar  jacket.  The  vacuum  case  contained  an  ex¬ 
change  gas  of  air  or  helium,  depending  on  the  temperature  region  being  investigated. 
Temperatures  between  room  and  liquid  nitrogen  temperature  woe  reached  by  means 
of  liquid  nitrogen  in  the  outer  dewar.  Liquid  helium  was  then  transfered  with  a  mini¬ 
mal  amount  of  exchange  gas  to  cool  down  to  helium  temperature  in  a  controlled  fashion. 
Pumping  on  the  exchange  gas  enabled  us  to  eontroll  the  temperatures  for  both  liquid 
nitrogen  and  liquid  helium  ranges.  The  temperature  was  measured  using  a  calibrated  sil¬ 
icon  diode  thermometer  mounted  close  to  the  sample,  and  wired  in  a  four  point  contact 
configuration. 

Measurement  of  the  resistance  of  the  sample  was  controlled  by  a  MACSYM  350 
(ANALOG  DEVICES,INC)  based  computer  system.  A  computer  program  was  used 
to  collect,  simultaneously,  the  susceptibility  and  the  c-axis  voltages.  The  sample  volt¬ 
age  was  first  read  with  sero  current  through  the  sample  by  a  channel  on  an  analog  in¬ 
put  card  in  the  MACSYM.  Then  the  sample  current  was  set  to  5  ma  by  a  channel  on  an 
analog  output  card  in  the  MACSYM.  This  value  of  current  was  chosen  as  a  tradeoff  be¬ 
tween  low  power  dissipation  in  the  dewar  and  low  noise  in  the  signal.  The  voltage  at  5 
ma  was  taken  and  then  the  current  was  switched  off.  The  sample  voltage  was  then  read 
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again.  The  average  of  the  zero  current  voltages,  representing  the  thermal  voltage  in  the 
junctions  of  the  leads,  was  subtracted  from  the  voltage  at  5  ma.  Finally  the  sample  re- 

( 

sistance  was  calculated  from  this  corrected  voltage.  The  data  then  were  trasfered  to  a 
Digital  Equipment  Corporation  VAXll  minicomputer  for  routine  analysis. 
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Karaite  and  Analysis 


The  data  plotted  in  Fig(l)  are  the  absolute  resistivity  versus  temperature  for  stages 
2, 3, 4, 5,0  of  FeCl3  GIC,  and  HOPG.  Several  features  can  be  observed  in  the  data:  1)  For 
stage  index  less  than  5  the  temperature  behavior  of  the  resistivity  is  metallic-like,  and^ 
the  resistivity  saturates  at  low  temperatures.  The  metallic  behavior  is  most  pronounced 
for  stages  2  and  3  and  less  so  for  stage  4.  2)  Note  that  for  stage  5  the  resistivity  ex¬ 
hibits  metallic-like  behavior  at  high  temperatures  and  crosses  to  activated  behavior  at 
low  temperatures.  3)  Observe  that  the  highest  stage  sample,  stage  9,  exhibits  activated 
behavior  throughout  the  entire  range  of  temperatures  shown.  These  behaviors  are  more 
easily  seen  on  a  plot  of  resistivity  normalized  to  room  temperature  versus  temperature, 
shown  in  Fig(2).  In  this  figure,  stage  number  increases  in  the  positive  y  direction  (ex¬ 
cept  for  HOPG).  This  presentation  supports  more  clearly  the  identification  of  stage  5  as 
the  stage  at  which  a  definite  transition  in  behavior  occurs  [13,24].  Stage  2  and  4  sam¬ 
ples  behave  qualitively  alike  whereas  stage  9  behaves  fundamentally  differently  and  is 
similar  to  HOPG  in  character  as  one  would  expect  for  higher  stages.  As  the  number  of 
pristine  graphite  layers  between  the  intercalate  layers  becomes  large  enough  to  screen 
out  the  out-of-plane  interactions,  the  c-axis  conduction  is  expected  to  be  controlled  by 
the  graphite  layers. 

The  absolute  resistivity  versus  the  inverse  of  the  stage  index  is  presented  in  Fig(3) 
at  293  K,  77  K,  4.2  K.  At  293  K  and  77  K  the  peak  in  resistivity  is  at  stage  5.  This  is, 
qualitatively,  in  agreement  with  the  results  of  the  magnetic  properties  measurements  for 
the  FeCls-graphite  system  [13].  In  those  measurements,  it  was  shown  that  the  system 
behaves  as  a  two-dimensional  system  and  stage  5  exhibits  maximum  characterstics  of  the 
two-dimensional  nature  relative  to  the  other  stages.  If  the  charge  carriers  are  confined 
between  the  planes,  one  would  expect  the  conduction  along  the  c-axis  to  be  minimum 
and  thus  the  resisitivity  in  this  direction  should  rise. 
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The  resistivity  temperature  coefficient  versus  the  inverse  of  the  stage  index  for  the 
temperatures  4  ,77,  and  293  K  are  presented  in  Fig(4).  These  coefficients  are  the  slopes 
of  the  resistivity  versus  temperature  curves  at  the  particular  temperature  normalized  to 
the  room  temperature  resistivity  values.  As  it  turns  out,  the  labeling  of  the  y-axis  gives 
the  change  in  resistivity  per  degree  kelvin  as  a  percentage  of  the  resistivity  at  room  tem¬ 
perature.  All  the  values  are  seen  to  be  small  in  absolute  magnitude-  less  than  one  half  of 
a  percent  per  degree,  It  is  the  sign  of  the  coefficient  which  is  important.  This  gives  the 
type  of  temperature  dependence-  positive  coefficients  indicate  "metallic”  temperature 
behavior  and  negative  coefficients  indicate  "activated”  behavior.  Again  the  fundamental 
difference  in  behavior  between  stages  with  stage  index  less  than  five  and  greater  than  5 
is  readily  apparent.  Also  note  that  for  all  stages  the  resistivity  coefficient  at  room  tem¬ 
perature  is  very  nearly  equal  to  zero. 

Athough  we  tried  to  compare  our  data  with  various  theories,  we  found  the  best 
qualitative  agreement  with  the  model  proposed  by  Sugihara  [6j.  In  this  model  it  was 
concluded  that  the  c-axis  conduction  in  acceptor  GIC’s  is  a  hopping  conduction  as¬ 
sisted  by  phonons  and  impurities.  The  total  conductivity  in  a  particular  stage  sample 
is  the  sum  of  the  conductivity  due  to  phonon-assisted  hopping  and  that  due  to  impurity- 
assisted  hopping.  The  temperature  dependences  of  the  two  hopping  mechanisms  differ. 
Phonon-assisted  hopping  conductivity  is  expected,  theoretically,  to  vary  as  T 2  for  low 
temperatures  and  to  be  constant  at  high  temperatures.  Impurity-assisted  hopping,  how¬ 
ever,  is  constant  at  low  temperatures  and  decreases  linearly  with  T  at  higher  temper¬ 
atures.  Phonon-assisted  hopping  dominates  the  conductivity  in  high  stages,  whereas 
impurity-assisted  hopping,  as  a  result  of  the  larger  number  of  defects  introduced  by  in¬ 
tercalation,  should  dominate  in  low  stages.  These  qualitative  results  have  been  shown 
to  agree  with  data  taken  for  SbCU  acceptor  GIC’s  [18]  and  our  data  for  FeCl3  GIC’s,  as 
shown  below.  The  model  of  Sugihara  [6]  can  be  summarized  by  the  following  equations: 


a)  Phonon-hopping  Conduction 


<re  =  Oep  +  <*c »•  (3) 

Here  <rcp,  crc,,  and  <rc  are  the  c-axia  phonon-,  impurity-,  and  total-assisted  hopping  con¬ 
ductivity,  respectively. 

We  have  used  the  above  theory  to  analyze  the  data  of  two  stages,  stage  2  as  an 
example  of  a  low  stage  and  stage  0  as  an  example  of  a  high  stage  of  FeCh  GIC’s.  For 
stage  0,  we  find  an  interesting  quantitative  agreement  of  our  data  with  theory.  The  data 
presented  in  Fig(S)  is  the  c-axis  conductivity  of  stage  9  versus  temperature  from  4  K  to 
70  K.  The  data  points  are  shown  by  the  squares.  A  fit  has  been  made  to  a  T 2  depen¬ 
dence.  It  can  be  seen  that  the  fit  is  excellent  from  10  K  to  about  00  K.  This  is  the  max¬ 
imum  temperature  for  agreement  with  theory  found  in  an  exhaustive  series  of  fits  for  a 
large  range  of  low  temperatures.  Similarly  for  the  data  of  stage  2,  a  fit  has  been  made 
to  a  T  dependence.  The  result  of  the  fit  is  shown  in  Fig(6).  The  data  points  are  shown 
by  the  squares.  As  shown  in  the  figure  the  data  agree  well  with  a  linear  T  dependence 
in  a  temperature  range  between  300  K  and  100  K.  Therefore  the  data  of  the  two  stages, 
which  are  chosen  as  examples  of  low  and  high  stages,  are  in  excellent  agreement  with  the 
prediction  of  the  Sugihara's  theory.  The  results  of  the  these  fits  are  as  follows:  In  stage 
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9  the’ data  were  compared  to  a  fit  to  the  equation  ffep  =  A  +  BTa,  where  the  fitting  pa¬ 
rameters  A  and  B  have  the  values;  A  =  0.530  Q-cm,  B  =  65.3  x  10~6  H-cm K~2.  In  stage 
2,  however,  the  impurity  contributions  are  expected  to  dominate  the  c-axis  conductivity, 
and  thus  the  data  were  compared  to  a  fit  to  the  equation  erei  =  C  - 1-  DT.  The  fitting 
parameters  C,  and  D  have  the  values  C  =  0.346  fl-cm,  and  D  =  2.92  x  10-3  fl-cm K~l. 

Now  all  the  data  shown  in  Fig(2)  can  be  analyzed  on  the  basis  of  the  Sugihara’s 
model.  The  data  of  stages  2  and  4,  classified  as  low  stages,  shown  in  the  figure  exhibit 
metallic-like  behavior  in  the  high  temperature  region  and  nearly  temperature  indpen- 
dence  at  low  temperatures.  Stage  2  is  expected  to  have  a  relatively  large  number  of  im¬ 
purities  in  comparison  to  that  of  stage  4.  This  would  explain  the  difference  in  the  tem¬ 
perature  dependence  of  the  resistivity  for  the  two  stages.  Stages  2  and  3  have  a  similar 
behavior.  Therefore,  the  temperature  dependence  of  the  c-axis  resistivity  for  the  lower 
stages  behaves  consistently  with  the  Sugihara’s  theory.  In  the  high  temperature  region 
the  data  for  stage  2  of  this  work  and  those  of  Issi  et  al  [11]  are  in  quantitative  agree¬ 
ment.  At  low  temperatures,  however,  subtracting  the  residual  parts  from  the  total  re¬ 
sistivity  leads  to  a  complicated  temperature  dependence.  The  overall  behavior  of  the  c- 
axis  resistivity  in  this  system  is  controlled  by  the  competition  between  two  mechanisms, 
namely  the  phonon-  and  impuriy-assisted  hopping  conductions.  The  data  of  Fig(2)  indi¬ 
cate  that  both  of  these  mechanisms  contribute  to  the  c-axis  resistivity  of  stage  5.  There¬ 
fore  the  resistivity  of  stage  5  is  expected  to  exhibit  both  the  activated  and  metallic  be¬ 
haviors.  Clearly  the  data  of  stage  5  in  figure  2  are  consistent  with  this  prediction.  The 
large  c-axis  resistivity  of  stage  5  is  due  to  the  delocalization  of  the  electrons  in  the  a-b 
plane  as  a  result  of  its  two-dimensional  nature.  The  c-axis  resistivity  of  higher  stages, 
•fcowu  far  stage  9  in  Fig(2),  shows  semiconductor-like  behavior  of  dominant  phonon- 
aaeteied  bopping  conduction  which  is  characterized  by  the  Ta  dependence.  This  is  also 


it  with  the  prediction  of  the  theory  for  the  higher  stage  samples.  The  c-axis  re- 


!W5C 


sistivity  of  higher  stage  samples,  as  expected,  exihibits  similar  behavior  to  that  of  the 
pristine  graphite.  The  data  of  Uher  and  Morelli  [18)  are  also  consistent  with  the  predic¬ 
tion  of  the  Sugihara’s  theory,  however,  none  of  their  stages  exhibits  intermediate  transi¬ 
tion,  similar  to  that  of  stage  5  of  this  work,  between  phonon  and  impurity  contributions. 

One  of  the  remarkeable  features  about  the  FeCU-GIC’s  is  the  existence  of  stage  5 
as  a  boundary  stage  between  low  and  high  stages.  Moreover,  the  magnetic  susceptibility 
data[l3]  indicate  that  the  two-dimensional  nature  of  this  system  has  exponetial  depen¬ 
dence  on  the  stage  index  and  it  was  shown  that  it  is  maximum  at  stage  5.  In  fact,  the 
conductivity  anistropy  (c re/<?at)  versus  the  stage  index,  shown  in  Fig.  (7),  has  a  sim¬ 
ilar  exponential  behavior  and  between  stage  2  and  stage  5  the  anistropy  is  maximum 
at  stage  5.  Once  the  system  comes  completely  into  a  two  dimensional  state,  the  out-of¬ 
plane  phonon  contributions  to  the  conductivity  decrease  and  thus  the  anistropy  is  in¬ 
creased.  This  results  again  emphasizes  the  two-dimensinality  of  this  system.  In  figure 
(7),  the  natural  logarithm  of  the  anistropy  (<re/Cab)  is  plotted  versus  the  stage  index.  As 
shown  on  the  figure,  the  data  for  stage  2  through  5  follows  an  exponential  law.  On  the 
same  figure  the  natural  logarithm  of  the  susceptiblity  peak,  at  the  reported  phase  transi¬ 
tion^, 13],  is  plotted  versus  the  stage  index.  As  shown  in  the  figure,  the  data  from  two 
independent  measurements  follow  the  same  general  exponential  behavior. 

In  order  to  see  if  there  are  hysteretic  effect,  thermal  cycling  of  the  resistivity  was 
also  performed  for  all  stages.  No  hysteresis  was  observed  in  any  of  the  stages  except  an 
insignificant  one  for  stage  9.  The  data  plotted  in  Fig(8)  are  the  scan  up  and  scan  down 
of  the  c-axis  resistivity  versus  temperature  for  stage  0.  We  do  not  expect  that  the  hys¬ 
teresis  shown  in  the  data  below  100  degrees  to  be  related  to  the  sample  structure  and 
attribute  it  to  the  nature  of  the  experiment.  In  this  temperature  range  the  resistivity 
drastically  varies  with  temperature.  In  the  experiment  it  was  difficult  to  control  the 
cooling  rate  to  collect  the  data  at  the  same  points  as  during  the  warming  rate.  There- 
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fore  the  FeCls-GIC’s,  unlike  the  SbCls-GIC’a,  has  no  significant  hesteresis  which  would 
be  related  to  microscopic  transition  mechanisms  such  as  the  pinning  of  discommensura- 
tions  caused  by  defect  sites  in  HOPG  [18].  The  high  temperature  order-disorder  transi¬ 
tion  [25]  which  occurs  near  room  temperature  in  stages  1  and  2  ICl-graphite  intercalated 
compounds  [17]  does  not  exist  in  the  FeClr graphite  compounds.  The  absence  of  the  his- 
teresis  and  such  high  temperature  structural  phase  transitions  and  the  long  period  room 
temperature  stability  suggest  that  the  FeCU  sytem  is  a  suitable  candidate  among  the  se¬ 
ries  of  GIC’s  for  further  structural  studies. 

CONCLUSION 

The  result  of  experimental  measurements  indicates  that  the  temperature  depen¬ 
dence  of  c-axis  resistivity  has  a  universal  behavior,  consistent  with  the  existing  theo¬ 
ries,  for  acceptor  graphite  intercalated  compounds.  Detailed  investigations  in  terms  of 
the  model  reported  by  Sugihara  have  shown  that  the  lower  stages  of  FeCl3-  graphite  in¬ 
tercalated  compounds  exhibit  metallic-like  conduction  along  the  c-axis,  whereas  higher 
stages  show  an  activated  behavior.  The  conduction  along  the  c-axis  in  this  system  is  a 
result  of  the  competition  between  two  mechanisms,  the  phonon-assisted  hopping  and  the 
impurity-assisted  hopping.  The  first  contributes  dominatly  in  higher  stages  at  low  tem¬ 
peratures,  while  the  second  contributes  in  lower  stages  at  high  tempratures.  As  a  man¬ 
ifestation  of  the  competition  between  the  two  hopping  mechanisms,  stage  5  behaves  as 
intermediate  stage  and  exhibits  nearly  constant  temperature  dependence  between  room 
and  liquid  helium  temperatures.  The  results  of  this  work  and  those  of  the  magnetic 
measurements,  reported  elsewhere  [13],  independently  support  the  existence  of  a  two- 
dimensional  state  in  this  system  with  stage  5  exhibiting  the  greatest  two-dimensional 
characteristics. 
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FIGURE  CAPTIONS 


Fig.  1.  The  c-axis  resistivity  vs  temperature  for  stages  2  (□),  3  (OM  (A),  5  (+),  9  (x) 
FeCU-GIC’s,  and  HOPG  (❖). 

Fig.  2.  The  c-axis  resistivity  normalized  to  the  room  temperature  values  vs  temperature 
for  stages  2  (□),  3  (OM  (A),  5  (+),  9  (x)  FeCl3-GIC’s,  and  HOPG  (0). 

Fig.  3.  Room  temperature  c-axis  resistivity  vs  the  inverse  of  the  stage  index  for  stages  2 
(□),  3  (OM  (A),  5  (+),  9  (x)  FeCla-GIC’s,  and  HOPG  (<>). 

Fig.  4.  Resistivity  coefficient  of  temperature  at  room  (□),  liquid  nitrogen  (0)»  and  liq¬ 
uid  helium  (A)  temperatures  vs  the  inverse  of  the  stage  index  for  stages  2, 4,5, 9  FeCla- 
GIC’s,  and  HOPG. 

Fig.  5.  The  fits  of  the  c-axis  conductivity  vs  temperature  for  stage  9  FeCl3-GIC’s.  The 
(□)  are  the  data  points. 

Fig.  6.  The  fits  of  the  c-axis  conductivity  vs  temperature  for  stage  2  FeCl3-GIC’s.  The 
(□)  are  the  data  points. 

Fig.  7.  Natural  logarithm  of  the  anistropy  (<re/oab)  (Cl)and  the  peak  height  of  the  mag¬ 
netic  susceptibility  (A)  vs  stage  index  for  FeCl3-GIC’s. 

Fig.  8.  Scan  up  (□)  and  scan  down  (O)  of  resistivity  vs  temperature  for  stage  9  FeCh- 
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APPENDIX  VI 

HAIL  EFFECT  AND  MAGNETORESISTANCE  IN  LOW  DIMENSIONAL 

MAGNETIC  SYSTEMS* 

A.K.  Ibrahim,  E.  Powers,  G.O.  Zimmerman,  and  M.  Tahar 
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Boston,  MA  02215 

ABSTRACT 

We  have  measured  the  Hail  effect  and  magnetorssistance  of  stages  5  and  0  FeClj- 
graphite  intercalated  compounds  in  a  magnetic  field  up  to  20T  and  at  several  temper* 
aturss.  Our  preliminary  results  indicate  that  the  Hall  Coeffieient(Rff)  of  stage*5  is  pos¬ 
itive  at  all  field  values  and  has  insignificant  variation  between  room  and  liquid  helium 
temperature.  In  stage-0,  howevere,  R*  has  a  complicated  field  and  temperature  depen¬ 
dence.  The  magnetoresistance  of  the  two  stages  at  all  temperature  and  field  values  is 
positive  and  stage-5  exhibits  strong  Shubinkov-daHaas  oscillations  at  helium  tempera¬ 
ture.  Fourier  transforms  of  these  peculations  are  performed  and  the  results  indicate  that 
stage-5  has  a  dominant  fundamental  frequency  which  is  attributed  to  the  hole  carriers, 
whereas  stage-0  exhibits  a  modulation  of  several  frequencies  which  are  related  to  differ¬ 
ent  carrier  pockets. 

INTRODUCTION 

Because  of  the  drastic  in  the  physical  properties  upon  the  transition  from  a 

three-dimensional  state  to  a  two-dimensional  one,  low  dimensional  systems  have  recently 
been  the  subject  of  numerous  research  efforts.  The  extremely  high  anisotropy  in  most  of 
the  physical  properties  of  highly  oriented  pyrolytic  graphite  (HOPG)  makes  this  system 
an  interesting  candidate  for  low  dimensional  studies.  Moreover  this  anistropy  can  be  in¬ 
creased  by  the  intercalation  of  acceptor  or  donor  compounds  into  the  graphite  (IGS’»)[I). 
Because,  electronically,  only  the  graphite  layers  next  to  the  intercalant  change  signifi¬ 
cantly  upon  intercalation  1-2] ,  the  intercalation  of  magnetic  species  into  the  graphite  can 


thus  provide  a  layered  magnetic  compound  with  variable  spacing  between  the  magnetic 
layer*.  Therefore  IGC’s  are  useful  system*  for  two-dimensional  magnetic  and  structural 
studias[3-4). 

In  the  last  few  years,  IGC’s  have  been  extensively  investigated  by  numerous  re¬ 
search  laboratories.  However  less  attention  was  paid  to  galvanomagnetic  (GM)  effect 
studies  in  these  compounds,  and  only  a  limited  amount  of  theoretical  and  experimen¬ 
tal  work  has  been  reported[S-6].  As  an  extension  to  the  comperhensive  investigation  of 
the  FeCl;j~graphite  intercalated  system,  which  was  begun  few  years  ago  at  our  labora¬ 
tory,  we  have  measured  the  Hall  effect  and  the  magnetoresistance  of  stages  5  and  0  of 
this  system(stage  index  refers  to  the  number  of  graphite  layers  between  two  consecutive 
intercalant  layers).  In  higher  stages,  where  there  are  enough  inner  layers  to  screen  out 
the  interactions  between  the  magnetic  intercalant  layers,  one  can  expect  the  behavior  to 
approximate  a  two-dimensional  magnetic  system.  The  measurements  were  carried  out  in 
a  field  up  to  20T  and  at  several  temperatures. 

both  the  long  period  air  stability!  7}  and  the  ease  of  the  technique  of  preparation 
[2]  of  FeCls-graphite  compounds  attract  researchers  to  investigate  this  system.  The 
measurements  of  the  magnetic  properties  of  this  system  were  performed  by  several 
authors(&-10],  the  results  of  the  magnetic  susceptibility  measurements  indicate  the  ex¬ 
istence  of  a  low  temperature  phaee  transition  in  the  fora  of  a  sharp  peak  in  the  suseep- 
tibilty.  Other  investigations  such  as  x-ray  structural(ll],  thermal  conductivity!  12|,  Ram- 
man  spectroscopic  and  infraredjlS),  and  relaxation  effact(M)  measurements  have  been 
reported.  However  there  are  no  comprehensive  studfcs  of  GM  effects  reported  for  the 
FeCls-graphite  system.  Quantum  oedllations(QO),  which  provide  information  about  the 
Fermi  surfaces  and  density  of  carriers,  were  observed  in  some  IGC*»fl5-lflj.  In  addition 
to  the  Hall  offset  and  magnetoceaistaaca  We  also  investigate  the  quantum  oscillations  in 
stages  5  and  0  of  FeCls  IGC’s. 


EXPERIMENTAL 

The  FtCls-GIC  samples  were  prepared  using  a  standard  two- sons  furnace  technique 
where  stage  index  was  controlled  by  the  temperature  difference  between  the  graphite 
host  (HOPG)  and  the  FeCls  powder.  The  samples  were  in  the  form  of  thin  rectangular 
plates  of  dimensions  l.SxO.&xO.l  cm3.  Well-staged  samples  were  achieved  by  controlling 
the  pressure  of  Clj  gas  inside  the  intercalation  tube,  as  well  as  the  partial  pressure  of 
FeCls  through  rigid  temperature  control.  After  intercalation,  the  samples  were  char¬ 
acterised  for  identity  and  uniformity  of  staging  using  x-ray  (001)  diffraction.  The  x-ray 
diffractograms  were  also  used  to  determine  the  e-axis  repeat  distance  I«  after  cycling  the 
samples  from  room  to  liquid  helium  temperature,  and  showed  that  the  cycling  did  not 
affect  this  staging  distance. 

Most  of  the  samples  measured  were  characterised  by  means  of  the  Mdssbauer  effect, 
which  can  reveal  the  ratio  of  iron  atoms  next  to  iron  vacancies  to  that  of  the  total  num¬ 
ber  of  iron  atoms  in  the  intercalant  layer  as  well  as  that  of  Fea+  to  that  of  Fe3+  [17], 

The  initial  measurements  were  performed  on  the  same  samples  used  in  the  Mossbauer 
investigations.  The  same  samples  and  others  were  characterised  periodically  over  a  span 
of  several  years,  having  bean  stored  at  room  temperature  in  a  dry  nitrogen  enviroment, 
and  showed  no  deterioration  over  a  period  of  up  to  three  years  (7).  The  Hendrichs-Teller 
analysis  technique(lt]  was  also  used  to  check  the  stags  disorder  of  stage  9  sample.  In¬ 
deed,  our  x-ray  data  were  fat  good  agresment  with  the  data  calculated  from  the  theo¬ 
retical  analysis  of  a  pure  and  wall-staged  stags  9  sample.  Gahranomsgnetk  measure¬ 
ments  were  performed  in  high  magnetic  field  at  the  Francis  Bitter  National  Magnet  Lab¬ 
oratory  (MIT).  Layers  of  about  l.Bx.Bx^Olcm3  were  peeled  from  the  samples  and  then 
mounted  on  flat  substrates.  The  five  probe  dc  technique  was  used  to  measure  the  Hall 
voltage  and  the  transverse  magnetoresistance. 
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Eesults  and  Analysis 
A-  The  Hall  Effect 

The  Hall  resistivity  versus  the  applied  magnetic  field  at  room,  liquid  nitrogen,  and  liquid 
helium  temperatures  is  shown  in  Fig.(l)  and  Fig.(2)  for  stage  S  and  stage  0,  respectively. 
As  shown  in  Fig.(l)  the  Hall  voltage  is  positive  and  nearly  linear  with  the  field  at  all 
temperature  and  field  values*  indicating  the  dominance  of  hole  carriers  in  stage  5.  Stage 
9,  however,  exhibits  a  different  character  at  room  temperature  since  the  Hall  voltage 
fluctuates  between  positive  and  negative  values.  The  Hall  coefleient,  thus,  has  a  compli¬ 
cated  field  dependence  at  room  temperature  indicating  the  near  equality  of  the  electron 
and  hole  densities.  The  Hall  coefficient  of  HOPG(19*20J  at  low  temperatures  exhibits 
similar  behavior  to  that  of  stage  9  at  room  temperature.  As  shown  in  Fig.  (2),  at  liquid 
nitrogen  and  liquid  helium  temperatures,  the  Hall  voltage  is  positive  at  all  the  field  val¬ 
ues,  however,  it  is  not  aa  linear  in  the  field  as  that  of  stage  5.  Therefore,  the  Hall  coef¬ 
ficient  of  stage  5  has  the  least  complicated  field  dependence  relative  to  those  of  stages  2 
and  9,  and  HOPGfsee  ref.20  for  data  on  HOPG  and  stage-2). 

At  this  stage  of  our  work,  we  have  not  performed  detailed  investigation  of  the  tem¬ 
perature  dependence  of  the  Hall  coefficient  or  the  magnetoresistance.  However,  qualita¬ 
tive  information  can  be  extracted  from  our  preliminary  results.  The  Hall  resistivity  data 
shown  in  figures  1  and  3  for  stagss  5  and  9  and  those  of  HOPG  and  stage-2(20]  indicate 
that  the  Hall  coefficient  of  stage  5  is  loos  sensthre  to  the  variation  with  temperature  rela¬ 
tive  to  the  other  stages  and  HOPG  aa  well.  The  Hall  coefficient  Rg,  for  a  single-carrier 
and  two-carrier  models.  Is  related  to  the  carrier  mobilities  and  the  carrier  densi¬ 
ties  n*,n«  by  the  following  equations 


p  -1  _  1 

fig  *  —  or  — — , 

(1) 

.  i  i 

n«e  i  (n*/vt«  +  m«/ M*)a 

(2) 
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According  to  Eq.(l),  the  Hall  coefficient  R#  would  be  constant  between  liquid  nitrogen 
and  liquid  helium  temperatures  only  if  the  carrier  concentration  is  constant  in  this  tem¬ 
perature  range.  However  Eq.(2),  the  two-carrier  model,  indicates  that  the  constancy  of 
Rjr  is  controlled  by  the  temperature  dependence  of  the  four  quantitiee  n*.  n«,  mk.  and 
or  at  least  the  carrier  concentration  and  the  mobility  ratio.lt  would  be  unrealistic  to 
assume  the  constancy  of  these  four  quantities. 

• 

Therefore,  based  on  the  field  and  temperature  dependence  of  the  Hall  Coefficient, 
the  conduction  in  atage-6  is  dominated  by  one  kind  of  charge  carrier.  Our  meesurements 
of  the  c-axie  conductivity  in  thie  system{21]  indicate  that  stage-i  has  the c- 
ajds  conductivity  among  the  other  stages  of  FeClfgraphite  compounds.  Moreover,  the 
temperature  coefficient  of  the  resistance  is  nearly  aero  in  stage  6  while  the  other  stages 
exhibit  posthre  and  negative  reeletivity  temperature  dependence.  These  results  suggest 
that,  in  stage  S,  the  carriers  are  confined  into  layered  planes  and  the  electronic  transport 
has  the  behavior  of  a  two-dimensional  system.  The  above  quantitative  descriptions  are 
also  supported  by  the  results  of  the  QO,  as  shown  in  part  B  of  thie  work,  and  the  ob¬ 
served  low  temperature  magnetic  phase  transition  in  this  system(10].  It  has  been  shown 
that(22],  thie  phase  transition  has  the  nature  of  a  two-dimensional  magnetic  system  and 
atage-S  among  all  other  stages  of  TeClrgraphite  compounds  exhibits  the  maximum  tran¬ 
sition  [23  J. 

The  behavior  of  tho  Hall  data  of  stags-9  is  diffiwent  from  that  of  stage-8,  it  has 
a  complicated  field  and  temperature  dependence  ehnilar  to  that  of  HOPG.  At  room 
temperature  the  Hall  voltage  rhenges  sign  foom  negative  to  positive  values  at  a  field  of 
about  4T.  Thus  the  single  carrier  model  of  Eq.(l)  is  not  applicable  at  this  temperature 
and  one  has  to  use  the  two-carrier  model  of  Bq.(2)  to,  qualitatively,  analyse  these  data. 
Because  of  the  difireneo  in  the  effective  masom  of  the  electrons  and  hoks,  one  can  ex- 
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p«ct  the  high  field  to  drift  holes  more  than  electrons.  Thus  at  room  temperature  holes 
contribute  dominantly  to  the  Hall  voltage  in  the  high  field  region  .  In  terms  of  Eq.(2), 
this  result  can  be  described  as  a  change  of  the  inequality  <  1  in  low  field 

to  £  >  ( jjj)a  in  high  field.  Thus  it  is  possible  that  one  could  have  positive  Hall  Coeffi¬ 
cient  even  when  <  1.  At  low  temperatures  the  Hall  voltage  is  dominantly  contributed 
by  holes  and  then  exhibits  very  weak  QO  at  helium  temperature.  Thus  QO  oscillations 
exist  in  the  Hall  voltage  of  stage-0,  stage-2 (20),  and  HOPG[20),  however  less  significant 
oscillations  are  observed  in  stage-8. 

B-Magnetoresistaace 

The  tranverse  magnetoresistivity  fr*  as  a  function  of  the  applied  field  at  room,  nitro¬ 
gen,  and  helium  temperatures  is  presented  for  stage  5  and  stags  0  of  FeClj-graphite 
compounds  in  Fig.  (3)  and  Fig. (4),  respectively.  The  value  of  Ap  increases  with  the  field 
at  all  temperatures,  there  is  no  sign  of  saturation  up  to  a  field  of  20T  in  stage-9,  but  in 
stage-8  there  is  a  weak  indication  of  saturation  at  very  high  fields.  At  room  and  nitrogen 
temperatures  the  two  stages  exhibits  similar  field  and  temperature  dependence,  however 
at  helium  temperature,  although  the  over  all  behavior  is  still  the  same,  the  sise  of  the 
change  in  Ap  is  different.  The  magnetoresistance  of  both  stages  oscillates  at  helium  tem¬ 
perature,  but  the  magnitudes  of  the  oscillations  in  stage  8  are  relatively  larger  and  faster 
than  those  of  stage-9.  The  magnetoresistance  of  HOPG,  and  stages  2,8,  and  9  exhibit 
similar  field  and  temperature  dependence. 

The  similarities  between  the  msgnetoresistance  data  of  HOPG  and  the  three  stages 
of  FeCls  sugpst  that  the  transport  mechanism  in  this  system  can  be  described  in  terms 
of  a  two-dimensional  modeL  The  scattering  of  the  charge  carriers  is,  thus,  an  in-plane 
mechanism  controlled  by  the  graphite  lattice  and  the  density  of  free  carriers.  The  in¬ 
crease  of  the  number  of  free  carriers  due  to  the  charge  transfer  between  the  graphite 
and  the  intercalant  is  accounted  for  by  the  sharp  decrease  in  the  magnitude  of  Ap  of  the 


three  stage*  relative  to  that  of  HOPG.  Magnetoresistance  data  are,  generally,  analyzed 
on  the  bask  of  a  simple  two-carrier  model  in  which  an  average  Hall  mobility  could  be 
obtained  from  the  formula 


(3) 


For  conductors  that  obey  a  one-carrier  model,  the  carrier  density(n)  determined  by 
nen'po  3  1  should  be  equal  to  that  determined  by  Eq.(l).  However,  these  simple  one- 
carrier  and  two-carrier  models  fail  to  give  adequate  results  for  many  IGC’s.  Details  of 
these  calculation*  will  be  given  in  a  later  paper. 

The  QO  of  the  magnstorssistance  of  stags*  8  and  0  are  plotted  as  a  function  of  the 
applied  field  in  Fig.(5)(note  that  the  data  of  stage-8  are  multiplied  by  a  factor  of  70). 

As  shown  in  the  figure,  the  sise  of  the  oscillations  relative  to  the  base  line  in  stage  5  is 
significantly  larger  than  that  of  stage-0.  The  oscillations  in  the  electric  resistivity  arise 
from  the  oscillations  in  the  relaxation  time  for  the  scattering  carriers.  Because  these 
oscillations  are  a  manifestation  of  the  periodic  variation  in  the  density  of  states  at  the 
Fermi  energy,  the  large  amplitude  of  the  oscillations  in  stage-8  can  be  attributed  to  a 
longer  in-plans  relaxation  time.  The  results  of  the  in-plane  conductivity  measurements 
also  indicate  that  stage-8  has  a  relatively  large  value  compared  with  that  of  stage-0.  In 
Fig.(6)  the  peculations  in  at  ago  8  are  plotted  as  a  function  of  the  inverse  field.  As  shown 
in  the  figure,  the  fandemental  frequency  is  modulated  by  other  very  weak  harmonics. 

As  it  was  pointed  out  fat  the  previous  section,  this  result  is  consistent  with  the  Hall  data 
which  indicate  that  fat  stage-8  the  transport  mechanisms  are  dominantly  controlled  by 
one  kind  of  charge  carrier. 


Fourier  transforms  of  those  oscillations  are  displayed  in  Fig.(7)  for  stage  8  and  stage 
0.  The  dominant  frequency  of  stage-8  which  is  attributed  to  the  hole  oscillations  occurs 
at  about  SMT.  In  stage-9,  as  shown  in  the  figure,  oscillations  of  different  frequencies  are 
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observed  At  0.55T,  4.4T,  And  8.8T.  Some  high  frequencies  of  insignificant  Amplitudes  are 
Also  obserred  in  stage-9.  It  is  clear  that  some  of  these  frequencies  in  stage-9  are  due  to 
the  HOPG  layers  in  the  syAem.  The  large  frequency  of  stage-5  indicates  a  similar  large 
size  in  the  extremal  Fermi  surface  cross-sectional  area  which  is  a  measure  of  the  carrier 
density  at  the  Fermi  energy. 

CONCLUSION 

The  results  of  galvanomagnetic  measurements  for  stages  5  and  9  of  FeCly  graphite  com¬ 
pounds  indicate  that  the  conduction  in  stage-5  is  dominanted  by  hole  carriers,  whereas 
stage-9  exhibits  a  complicated  field  and  temperature  dependence  similar  to  that  of 
HOPG.  Both  the  Hall  data  and  the  magnetons istance  data  suggest  that  transport 
mechanisms  in  FeCl*  hare  the  behavior  of  two-dimensional  magnetic  systems.  Quan¬ 
tum  oscillations  in  the  magnetorasistaace  are  observed  in  both  stages,  however  the  size 
of  the  oscillations  fat  stage-5  is  very  large  relative  to  that  in  stage-9.  Fourier  tranfbrms 
of  these  data  show  a  dominant  single  frequency  in  stage-5  and  the  modulation  of  several 
frequencies  in  stage-9. 
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Figure  Capdsns 

Figure  1  The  Hall  resistivity  versus  the  applied  field  of  stage-5  at  helium(l),  nitrogen (2), 
and  room(3)  temperatures. 

Figure  2  The  Hall  resistivity  versus  the  applied  field  of  stage-9  at  helium(l),  nitrogen (2), 
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and  room(3)  temperatures. 

Figure  3  The  magnetoresistance  versus  the  applied  field  of  stage-5  at  he- 
lium(l),nitrogen(2),  and  room(3)  temperatures. 

Figure  4  The  magnetorssistance  versus  the  applied  field  of  stage-0  at  he- 
Uum(i),nitrogen(2),  and  room(3)  temperatures. 

Figure  5  The  oee illations  in  the  magnetorssistance  versus  the  field  at  helium  temperature 
of  stage-S(l)  and  stage-0(2) (Stage  5  data  are  multiplied  by  70). 

Figure  0  the  oecilattions  in  the  magnetoresistance  versus  the  inverse  ofield  of  stage- 
5(nonoscillatory  part  is  subtracted). 

Figure  7  Fourier  transform  of  data  obtainsd  from  stage-5  and  stage-0  at  helium  tempera¬ 
ture.  The  top  x-aads  refers  to  stage-0  data,  where  the  bottom  one  refers  to  stage-5  data. 
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CALVANOMAGNETIC  EFFECTS  IN  GRAPHITE  AND  GRAPHITE  INTERCALATED  COMPOUNDS 


A.E.  Ibrahia,  R.  Powers,  and  G.O.  Zimarman*4, 

Phyalca  Department,  Bo » con  University 

INTRODUCTION 

Electronic  properties  of  highly  orieneed  pyrolytle  graphite  (KOPG)  can 
be  draseically  changed  by  eha  intercalation  of  acceptor  or  donor  coapounda. 
The  charge  denalty  transferred  free  or  into  the  graphite  layers  plays  an 
important  role  la  controlling  the  electronic  properties  of  the  graphite 
intercalated  coapounda  (CXC's).  Although  tbs  transport  properties  of  chase 
coapounda  have  been  extensively  studied,  the  roles  of  the  charge  transfer 
between  the  latarealaats  sad  the  graphite  are  only  partly  understood. 
Galvanoaagnetlc  (CM)  aaaaureaants  and  the  related  quantum  oscillations  (QO) 
are  among  those  reliable  techniques  which  provide  information  about  the 
electronic  structure  of  these  compounds.  In  this  work,  wa  have  measured  the 
CM  effects  in  the  atage-2  FeClj-graphite  compound  and  in  HOPC.  The  aaaaure- 
mencs  were  carried  out  la  fields  up  to  20T  at  several  teaparacuraa.  in  high 
magnetic  field  at  hellim  takers  turn,  QO  were  observed  In  the  aagnetorests- 
taaca  (do I  and  in  the  Hall  effect  for  atage-2  tdisreas  HOPC  exihiblts  oscil¬ 
lations  only  in  AD.  These  angne to resistance  oscillations  of  HOPC  are  la 
qualitative  agreement  with  previously  reported  data  [1,2]. 

UPEIXtKNTAL 

FeCli-GIC  series  warm  prepared  ualag  a  standard  two-sone  furnenca 
technique  [3]  where  stage  index  was  controlled  by  the  eeaperature  difference 
between  the  graphitn  host  (HOPC)  and  the  FeClj  powder.  The  sables  ware  la 
eha  form  of  thin  rectangular  plates  of  diaensions  1.3x0.5x0.01  cm3.  The 
five  probes-K  technique  was  used  to  aeaeurn  the  Hall  effect  and  the  uapMt- 
o res is tames.  The  current  was  kept  lower  than  3  KA  to  avoid  thermal  effects.  . 
The  measurements  were  carried  out  at  the  National  Magnet  Laboratory  (MIT). 
Maaauraneaes  have  bean  nads  for  two  spacinena  cut  from  each  sa^la  (HOPC  and 
stage-2)  and  several  scans  were  taken  to  ensure  ths  reproducibility  of  our 
dace. 

RESULTS  AND  ANALYSIS 
A- THE  HALL  EFFECT  DATA: 

The  tOl  resistivity  of  HOPC  verses  the  applied  magnetic  field  at  room, 
liquid  nierogem,  and  liquid  helium  teapera cures  is  shown  in  flg.(l).  At  room 
eeaperature  the  Hall  Voltage  is  negative  at  all  field  values  indicating  that 


Fig.  1  Hall  resistivity  vs.  field'  (H)  Fig. 2  Hall  resistivity  vs.  field  (H) 
of  HOPC  at  He(l),  Hi (2)  and  ta(3)teup.of  stage-2  at  Had)  and  ta<2)  eeap. 
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ch«  majority  of  charga  carriers  art  electrons.  thia  ratult  Is  consistanc  with 
cha  low  field  oaaauraaancs  of  Ollloa  at  al  [<.]  *nd  Ayacha  [5],  howavar  at 
high  fialda  our  data  axlhible  a  saturation  affact.  Tha  alactron  and  hola 
danaltlas  of  HOPG  ara  naarly  aqual,  but  dua  to  tha  diffarancaa  In  thalr 
affactiva  aaasaa  tha  alactrons  ara  cha  dominant  contributors  to  tha  Hall 
Voltaga.  As  cha  flald  bacomas  larga  anough  to  forea  all  carrlars  to  complete 
chair  cyclotron  oscillations,  cha  Hall  Voltaga  saturatas  at  a  nagatlva  value. 
The  striking  result  at  liquid  nlcrogan  Temperature  la  cha  reversal  of  the 
Hall  Voltage  which  becomes  positive  (negative  Hall  Coefficient  fO  at  a  field 
of  about  12T.  Stellar  behavior  was  also  observed  at  low  flelda  (H<  0.1T)  by 
Dillon  at  al  [4] .  The  low  field  data  were  described  as  nob 11a  minority  hola 
behavior  which  also  depends  in  part  on  tha  quality  of  cha  sample  [6] .  In  cha 
high  flald  limit  at  nlcrogan  Temperature  holes  are  tha  dominant  contributors 
to  the  Hall  Voltage.  At  helium  Te^erature  cha  Hall  Voltaga  saturates  at  a 
vary  small  negative  value  similar  to  that  at  room  cetera cure  indicating  a 
dominant  alactron  contribution. 

Our  data  from  cha  Hall  affact  measurements  of  seaga-2  fed 3  CIC  indicate 
differences  between  staga-2  and  HOPG  atamariaad  ae  follows: 

1- Ia  .staga-2  tha  Hall  Voltaga  la  positive  at  all  temperature  and  field  values 
Fig  (2)  shows  tbs  Hall  resistivity  of  stage-2  at  helium  and  room  temperatures 
as  a  function  of  tha  applied  flald.  Aa  shown  in  tha  figure  the  Hall  algnal 

la  smaller  ae  lower  temperature  indicating  a  decrease  la  the  scattering  rata 
of  the  charga  carrlara. 

2-  la  staga-2  cha  aa^Utude  of  the  Hall  signal  is  comparable  to  that  of  tha 
magnatorsalatanca  voltage.  Tha  Hall  data  of  aeage-2,  therefore,  are  more 
reliable  than  those  of  HOPG  in  which  tha  vary  small  Hall  sigial  is  masked  by 
a  relatively  large  magsatoresleemea  voltaga. 

3-  Tha  Hall  Voltage  of  a rage-2  exlhlblcs  shablnhov-de  Haas  (SdH)  quantum 

oscillations  at  low  temperature  and  high  flald.  In  fig  (3)  tha  Hall  resist¬ 
ivity  is  shown  aa  a  function  of  the  applied  field  at  liquid  helium  Ta^er- 
ature  for  both  seaga-2  and  HON.  Aa  shown  in  tha  figure  the  oscillations 
occur  la  seaga-2  la  a  field  am  low  a a  HT,  whereas  HOPG  has  oscillations  of 
much  lower  frequencies  at  lower  fields.  Tha,  eslscaaea  of  the  QO  ae  low 
temperature  errlndea  the  interpretation  of  tha  data  aa  a  one- carrier  model 
and  therefore  cha  Hall  Coefficient  muac  be  given  by  the  equation  of  tha  two 
carrier  modal,  *H  •  1  Pvff  -  nu| 

ec  (Pup  ♦n»i3)2 

whan  P  and  a  denote  the  hola  and  alactron  carrier  danaltlas  Up  and  %  an 
tha  corresponding  nobilities.  Based  on  a  two-carrier  nodal  tha  determination 
of  tha  four  quantities  P,  a.  Up.  and  Up  from  the  CM  data  la  not  possible. 


Pig.  3  Hall  resistivity  vs.  flold  (H)  Pig.  4  Magna  to  resistance  ve.  flald  (B) 
of  HOPG  Cl)  and  stagm-2(2)  at  Be  tamp. of  HOPG  at  Ha(l),  Rl(2)  and  Bm(3)  taap. 


however  ea  described  above  qualitative  information  about  the  charge  carriers 
and  the  carrier  density  can  be  deduced. 

B-MAGRETORESISTANCE  DATA 

The  transverse  magnetorealscance  Ao/o  as  a  function  of  the  applied 
field  at  room,  SI,  and  He  temperatures  for  HOPC  Is  presented  In  Fig  (4).  Our 
data  qualitatively  agree  with  the  high  field  data  of  Lowrev  and  Spain  [7] 

At  all  ceaperature  values  As  increases  continuously  with  no  sign  of  saturation 
up  to  a  field  of  20T.  The  magnetoresistance  is  large  and  significantly  en¬ 
hanced  at  lower  temperatures  indicating  the  dominance  of  Impurities  and  im¬ 
perfection  scattering.  The  magneto resistance  dace  of  stage-4  is  shown  in 
fig  (5).  In  this  figure  the  transverse  component  Ao/o  is  plotted  versus 
the  applied  field  at  room,  Ni,  and  Ha  tenparatures.  At  room  and  Hi  tempera¬ 
tures  As  of  stage-2  sad  that  of  HOPC,  have  a  typical  field  dependent  charac¬ 
ter.  This  similarity  between  tha  data  for  HOPC  and  stage -2  supports  the  two- 
dimensional  modal  of  tha  transport  -nehaniam  in  chase  intercalated  confounds . 
Therefore,  tha  scattering  of  tha  charge  carriers  is  an  in-plane  mechanism 
controlled  by  tha  graphite  lactlea  and  tha  density  of  free  carriers.  The 
increase  of  tha  number  of  free  carriers  due  to  the  charga  transfer  between 
the  graphite  sad  tha  lnterealaat  la  accounted  for  by  tha  decrease  in  tha 
magnitude  of  Ao  of  at age- 2  relative  to  that  of  HOPC.  At  Ha  ceaperature  QO 
for  sc age -2  are  readily  observed  at  fields  of  H>8T.  Tha  comparison  between 
the  oscillations  in  fig  (3)  and  those  of  fig  (4)  indicate  that  tha  SdH 
frequencies  of  scaga-2  are  greater  by  saveral  orders  of  magnitudes  chan  those 
of  HOPC.  This  difference  can  be  intarpracad  as  an  increase  in  tha  Fermi  sur¬ 
face  cross-sectional  areas  as  a  result  of  tha  charge  transfer  between  tha 
graphite  and  tha  Incarcalant.  It  is  Interesting  to  nota,  shown  in  fig  (6), 
tha  damping  of  the  S<*  oscillations  as  tha  ceaperature  increases.  A a  shown 
in  tha  figure  the  oscillations  persist  at  teaparaturma  as  high  as  33K.The 
existence  of  the  oscillations  at  this  teaperaturm,  a  takers  turn  which 
requites  sufficient  perfection  to  satisfy  tha  quantum  condition  u»cT»l, 
indicates  that  tha  FaC^  intercalated  compound  la  an  in  cares  ting  system  for 
further  study. 
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The  behavior  of  a  system  of  dipoles  is  presented.  Here  we  study  a  system  of  dipoles 
located  on  a  honeycomb  lattice  under  the  mean-field  approximation.  The  dipoles  are 
confined  to  rotate  in  the  plane  of  the  lattice.  It  is  found  that  in  zero  external  field,  the 
system  has  a  ground  state  which  is  continuously  degenerate.  This  degeneracy  persists  to 
T  jt  0  in  zero  field.  For  non-zero  planar  field  pointing  in  an  arbitrary  direction  in  the 
plane,  we  compute  numerically  the  magnetization,  susceptibility,  and  the  specific  heat  as 
functions  of  both  temperature  and  field.  Our  results  show  the  occurence  of  cusps  in  these 
quantities,  indicating  the  onset  of  phase  transitions.  On  the  basis  of  these  calculations  a 
phase  diagram  for  the  system  is  then  constructed. 
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1- Introduction 

This  study  is  motivated  by  recent  experimental  measurements  [li  of  the  thermal  and 

magnetic  properties  of  the  FeCh-graphite  intercalated  compounds.  The  intercalate  forms 

monomoiecular  layers  of  FeCh  sandwiched  between  graphite  layers  with  the  Fe3+  ions 

arranged  on  a  honeycomb  lattice.  The  iron  ions  are  in  the  spin-|  state  with  nearest 

neighbor  distance  of  3.5 A9.  There  is  also  evidence  that,  at  low  enough  temperatures,  the 

Fe3+  dipoles  are  confined  to  rotate  in  the  intercalate  plane.  It  is  thus  reasonable  to  assume 

that  dipolar  interactions  play  a  significant  role  in  this  substance,  with  the  intercalate 

« 

described  by  a  system  of  planar  dipoles  arranged  on  the  honeycomb  lattice.  In  this  paper 
we  carry  out  the  theoretical  investigation  of  the  thermal  and  magnetic  properties  of  such  a 
dipole  system  under  the  mean-field  approximation.  Results  of  experimental  measurements 
as  well  as  comparison  of  the  theoretical  predu  .ion  with  experiments  is  given  elsewhere  [1]. 

The  outline  of  this  paper  is  as  follows:  The  problem  is  described,  and  formulated  in 
section  II  using  the  mean-field  approximation.  The  solution  of  a  special  case,  in  which  all 
dipole  moments  are  of  equal  magnitude  and  interact  with  nearest  neighbor  interactions, 
is  analyzed  in  section  III.  The  ground-state  in  zero  field  is  then  studied,  with  interactions 
extending  to  all  neighbors  in  part  a)  of  section  IV.  The  most  general  solution  of  the  mean- 
field  equations,  for  nearest  neighbor  interactions,  are  studied  numerically  in  part  b)  of 
section  IV,  and  in  section  V  we  discuss  the  resulting  phase  diagram. 


II-  Mean-Field  Formulation 

Consider  a  system  of  N  dipoles  on  a  honeycomb  lattice  in  an  external  magnetic  field 
H.  We  assume  the  dipoles  interact  by  means  of  the  dipolar  interaction 

~  «/«/[£*  '  frj  —  3(£t  •  f «> ) ]  (l) 

with  all  dipole  moments  and  the  applied  field  H  confined  to  the  plane  of  the  lattice,  and 
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write  the  Hamiltonian  as 


*=  £  v„- i2) 

<i.j>  i 

Here  the  summation  is  taken  over  ail  pairings  of  dipoles,  fit  is  a  dipole  moment  of 

unit  .magnitude  at  the  ith  lattice  site,  and  is  the  coupling  constant  between  dipoles  /I, 
and  tf} .  The  experimental  situation  corresponds  to  Jt)  >  0. 


The  free  energy  F[p)  of  the  system  is  given  by  the  expression 


0F(p)  =  Tr 


p{3  Uij  -  h  •  fii  +  tnp) 

<*,j>  » 


(3) 


where  3  =  1/kgT,  h  —  3H  and  p  is  the  density  matrix.  In  the  mean-field  analysis  we  seek 


to  minimize  the  free  energy  in  the  subspace 

N 

P-\{pi  (4) 

■  a) 

where  p,  is  the  single  particle  density  matrix  normalised  to 


Tw  =  1. 


Substitution  of  (4)  into  (3)  and  the  use  of  (1)  leads  to 


=  £  k„ 


m.  •  m>  -  3(m,  •  r„)(m;  •  ftJ) 


<*,!> 

-h  •  £  A,  +  Tr(piinpt) 


where  K%]  =  3Jt),  and 


=  i  —  1,2,3, 

is  the  magnetisation  of  the  dipole  fii. 


(5) 


(6) 

(7) 


It  is  straightforward  to  earry  out  the  minimization  dF/ dpi  =  0  subject  to  the  contraint 
(5).  This  procedure  leads  to 


Pi  = 


2x/o(o ,) 


,<*.  a* 


(8) 
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and,  consequently, 


m,  =  m(ax)al 


where 


m(a)  = 


/.(«) 

/o(a)’ 


and  Jn(a)  is  the  modified  Bessel  function  of  order  n.  Here,  the  vectors  d,,  i=l,  2,  N. 
are  determined  from  the  self  consistent  equations 

=  h  —  y  *  \hj  —  3r,j(flj  •  ftJ)  ,  (11) 


where  the  prime  over  the  summation  sign  denotes  the  restriction  of  j#i.  Our  goal  now  is 
to  solve  for  a,  from  (11).  Then  the  per*dipole  magnetization  is  given  by 


l  N 


from  which  we  can  directly  compute  the  magnetic  susceptibility  as 


The  per-site  entropy,  given  by  the  last  term  in  (6),  can  now  be  obtained  by  combining  (7) 


5  =  -ks  a*m(a»)  ~  in(2ffj„(a,)) 


and  (8)  as 


from  which  we  can  compute  the  per-site  specific  heat: 


c  =  r— . 

dT 


So  far  we  have  been  careful  in  the  formulation  to  allow  dipolar  interactions  between 
all  pairs.  For  simplicity  we  now  assume  nearest- neighbor  interactions.  Furthermore,  we 
look  for  solutions  of  (ll)  with  periodic  ordering.  An  honeycomb  lattice  of  N  sites  can 
be  decomposed  into  N/6  nonintersecting  elementary  hexagons,  each  of  which  containing 


6  sites  oriented  and  numbered  as  shown  in  Fig.l.  A  priori,  and  guided  by  our  numerical 
solutions,  we  expect  the  magnetization  of  the  six  sublattices  to  be  different.  Allowing  this 
possibility  then,  for  periodic  ordering.  (11)  reduces  to  six  equations  as  follows: 


3' 


S9  =  h-  K  ^  m(aa») 
3'=v 


ag>  -  3fa0>(as>  •  ra&.) 


a  =  1,2,3 


aff  -  3rft-4( iff  •  fa>a) 


,  a'  =  1\2',3' 


(16) 


aa>  =  h  -  K  ^  m(a^) 

where  K  =  PJnm  Jnn  being  the  nearest  neighbor  interaction,  and  fQ^.,  is  the  unit  vector 
pointing  from  site  a  to  site  /?'.  For  fixed  T  and  H  we  solve  (16)  for  aa  and  aa>.  Then  the 
magnetization  is  given  by  (12).  In  fact,  using  the  identity 

3  -  3 

22  •(A>  fafi.)  =  -A, 

ftal  1 

which  is  valid  for  any  vector  A,  we  find  by  directly  summing  over  (16), 


2 

mi*  +  mj*  +  rfiji  =  —  (ij  +01+03-  3 h) 
3a 

2  _ 

mi  +  rfij  +  m3  =  —  (01*  +  a]*  +■  S31  —  3  h), 
3n 


and,  hence  using  (12), 


3 

l 

Oat 


3' 


(17) 


(18) 


(ID) 


#'*i’ 


in*  Solution  of  the  Special  Caee  of  Equal  Dipole  Moments 


In  this  and  subsequent  sections,  we  present  results  on  numerical  solutions  of  (16)  and 
the  resulting  determination  of  the  magnetisation  and  susceptibility.  For  the  purpose  of  il¬ 
lustration,  however,  we  analyse  in  this  section  the  solution  of  (16)  when  the  magnetisations 
ma  and  ma*  have  equal  magnitude.  This  is,  we  assume 

aa  =  Oa*  =  a.  (20) 
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While  the  assumption  in  (20)  simplifies  the  analysis,  it  must  be  pointed  out  that,  as  will 
be  seen  in  Sec.  IV  below,  solutions  other  than  (20)  do  exist  and  the  analysis  restricted  to 
(20)  is  applicable  only  in  a  limited  region  of  the  phase  space,  including  the  case  of  H  =  0 


Consider  first  the  ground  state  for  H  =  0.  By  symmetry  we  expect  the  dipoles  to  be 
arranged  in  a  symmetric  fashion,  namely,  the  dipoles  mi<  M2,  M3,  and  mi',M2',  M3< ,  making 
120°  with  respect  to  each  other,  although  there  can  be  relative  rotations.  A  straightforward 
calculation  shows  that  the  ground  state  is  obtained  when 


$3  *  B2  -  120°  =  fix  -  240° 


By  =  By  -  120°  =  Bx>  -  240° 


(21) 


where  Bt(B[)  is  the  azimuth  angle  of  the  dipole  at  site  t(i').  [CF.  Fig.  1).  Any  other  relative 
positioning  of  the  two  sets  of  dipoles  yields  a  higher  energy.  When  (21)  is  satisfied,  we  find 
the  per-site  energy 

E  =  ~  JnmCO»(Bt  +  #i*  +  120°).  (22) 


clearly,  the  ground  state  is  highly  degenerate  and  is  attained  by  setting 

Bi  +  Bit  =  240°,  for  Jnn  >  0, 


(23) 


Bx  +  By  =  60°,  for  Jnn  <  0 

Indeed,  numerical  simulations  show  that  the  conditions  (21)  and  (23)  also  lead  to  the 
lowest  energy  when  interactions  are  extend  to  all  neighbors,  with  each  dipole  interacting 
with  every  other  one. 


For  T  #  0,  it  yt  0,  and  using  (20)  and  (10),  we  obtain  from  (18)  the  constraint 


+  fj  +  ij  **  f|i  +  A)'  -4*  O31 


3  h 


'  1  -  §C(«) 


(24) 


where 


C(a)  =  Km(a)/a. 

Note  that  (24)  implies  «i  +  03  +  £3  =  0|«  +  iy  +03-  =0  in  zero  field. 
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(25) 


We  can  further  iterate  (16)  to  obtain  equations  for  aa  and  aQ<  individually.  After 
some  algebra  and  the  use  of  identities  (17)  and 


?a3'  •  fs't  ~  jU  3£a-,), 


we  arrive  at  the  result,  assuming  the  constraint  (20), 


3  3' 


A(a)oa  =  B(a)h  +  D{a)  ^  ^  fai3>{da  •  r30.) 

$'  =  v 


where 

^(®)  =  1  -  \C(a) 

Btd\  -  A(a) 

(  )_  1-1  C(a) 

D(«)  =  ”|C(a)p,  (28) 

and  an  equation  similar  to  (27)  for  aa>  with  the  primed  and  unprimed  subscripts  inter¬ 
changed. 

Each  of  these  two  equations  represents  a  set  of  six  equations  containing  four  unknowns, 
e.g.,  the  parameter  a  and  the  three  directional  angles  of  aa  and,  therefore,  is  generally  a  set 
of  over-determined  equations  which  does  not  possess  a  solution.  However,  when  .4(a)  =  0 


or,  equivalently, 


I  \C(a)  |=  1, 


which  also  happens  to  be  the  characteristic  equation  of  (27).  It  can  be  shown  that  (27) 
reduces  to  a  set  of  two  independent  equations.  While  in  this  case  (29)  and  (27)  are 
insufficient  to  determine  S9  uinquely,  the  overall  constraint  (24)  still  holds.  Thus  from 
(29)  and  (24)  we  And 

9- 

®1  +  «3  +  *3  =  ®1'  +  ®3'  +  ®3'  —  r  h,  Jnn  >  0 

l-  ,  <“> 

—  -h,  Jnn  <  0. 
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Equation  (30)  for  either  Jnn  >  0  or  Jnn  <  0  is  a  set  of  two  independent  equations 
and  is  equivalent  to  (20).  Thus,  for  fixed  T  and  H,  we  use  (29)  to  compute 


a  =  ao(T). 


Then  any  set  of  £ i  and  S\>  whose  magnitudes  are  ao(T)  and  satisfying  (30)  is  a  solution. 
The  constraint  (30)  now  leads  to  the  existence  of  helical  and  nonhelical  phases  similar  to 
those  discussed  in  Ref.[3j.  The  critical  temperature  Tc  is  given  by  ao{Tc)  =  0  or,  using 
(29)  and  (25), 


Tc  = 


The  phase  boundary  between  the  helical  and  nonhelical  phases  is  £i  +  £2  +  <»3  =  a,  or, 


using  (30)  and  (31), 


flo(r)  =  2kT'  ^nn  >  ^ 


Qlf 

*o(T)  =  Jnn  <  0. 


The  border  between  the  paramagnetic  and  nonhelical  phases  is  «t  +  £*  +  £3  =  3£,  or, 


ZH 

ao{T)  —  — Jnn  >  0 


ZH  . 

ao(T)  —  ^2*'  J^n  ^  0- 


Plot  of  the  phase  diagram  is  similar  to  that  of  Fig.2  of  Ref.  [2)  and  will  not  be  repeated 


here. 


IV*  Results  of  the  Numerical  Solution 

In  this  section  we  present  results  of  numerical  solutions  of  (16)  without  the  restriction 
(20)  of  equal  dipole  magnitude. 

a)  Ground  State  Energy  in  Zero  Magnetic  Field 

The  ground  state  calculation  for  zero  magnetic  field  was  performed  on  a  set  of  72 
dipoles.  The  dipoles  were  located  on  a  rectangle  which  measures  9  x  6%/3  where  the 
nearest  neighbor  distance  is  taken  as  one.  In  order  to  eliminate  effects  of  the  boundary, 
this  set  was  surrounded  by  48  similar  dipole  arrays.  The  field  of  all  the  48x72=3456 
dipoles  as  well  as  the  other  71  dipoles  within  the  original  rectangle  was  computed  at  the 
site  of  each  of  the  72  dipoles  in  the  centers!  rectangle. 

For  nearest  neighbor  interactions  we  verified  that,  the  ground  state  energy  is  given 
by  (22),  i.e.  the  assumption  (20)  of  equal  magnetization  magnitudes  was  correct.  We 
have  further  computed  the  zero-field  ground  state  energy  for  interactions  extending  to  all 
neighbors.  This  lead  to  the  same  ground  state  configurations  given  by  (21)  and  (23),  but 
with  a  per-dipole  energy  -2.2269  J„n,  which  was  about  \%  higher  than  the  value  of  -2.25 
Jnn  implied  by  (22)  for  nearest  neighbor  interactions.  Furthermore,  our  numerical  results 
indicated  that  the  ground  state  with  extended-neighbor  interactions  was  again  degenerate 
with  respect  to  relative  rotations  of  the  two  sets  of  dipoles  m,  and  £•'  **  long  as  both  (21) 
and  (23)  were  obeyed. 

b)  Mean- Field  Cakulatksi 

For  general  J?  and  T  ^  0  (16)  was  solved  by  means  of  a  numerical  iterative  method 
starting  from  states  specified  by  (21)  and  (23),  with  some  chosen  initial  values  of  aa- 
Random  initial  orientations  were  also  tried  but  they  did  not  lead  to  new  solutions.  It  turns 
out  that  there  exist  solutions  to  (16)  with  different  final  dipole  configurations  and  slightly 


different  free  energies,  typically  one  tenth  of  kgT.  We  identify  the  one  associated  with  the 
lowest  free  energy  as  the  true  dipole  configuration.  The  condition  for  self-consistency  was 
that  successive  values  of  aa  he  within  0.1 %  of  each  other.  An  accuracy  of  0.01%  was  also 
tried  with  no  significant  change  in  the  values  except  for  a  longer  computational  time.  For 
low  temperature  and  low  field,  convergence  was  achieved  after  as  few  as  5  iterations,  while 
near  the  phase  boundaries,  (see  later  discussion)  it  took  up  to  900  iterations  to  converge. 

For  H=0,  our  numerical  solution  again  verified  that  assumption  (20)  of  equal  magne¬ 
tization  magnitude  is  correct  for  T  ^  0  and  reproduced  the  solution  given  in  section  III. 
Particularly  we  verified  that  oa(r)  =  0  and  aar(T)  =  0  for  T  >Te.  However,  for  H  #  0, 
the  numerical  solutions  yielded  unequal  magnetizations  aQ  and  aa/.  In  general,  magnetic 
fields  applied  along  an  axis  of  symmetry  tended  to  pair  up  the  dipole  magnetisation  val¬ 
ues.  This  pairing  is  shown  in  Table  I.  By  symmetry,  the  behavior  is  identical  modulo  60°, 
although  the  pairing  will  vary.  We  thus  show  only  the  pairing  at  0°  and  90°.  The  pair 
of  dipoles  with  equal  magnitude  orient  themselves  with  positive  and  negative  angles  with 
respect  to  the  field.  This  pairing  was  broken  if  H  was  not  along  a  direction  of  symmetry. 

The  dipole  configurations  corresponding  to  the  lowest  free  energy  are  shown  in  figures 
2  and  3  for  various  temperatures,  and  fields  along  the  horizontal  and  vertical  directions, 
respectively.  In  these  and  subsequent  figures,  T*  is  the  temperature  given  in  units  of 
kgT/Jnn  and  h*  is  the  field  in  units  of  H[Jnn.  One  observes  that  at  low  temperature, 
i.e.  frame  "a*  of  each  figure,  the  dipole  pattern  resembles  one  of  the  zero  field  ground- 
state  patterns  specified  by  (21)  and  (23).  Thus  the  application  of  the  magnetic  field  lifts 
the  degeneracy.  Moreover,  the  configuration  with  the  lowest  free  energy  depends  on  the 
direction  of  the  magnetic  field.  Although  the  difference  in  the  temperature  between  the 
last  two  frames  in  Figs.2  and  3  is  small,  the  patterns  change  drastically.  This  change  is  due 
to  the  180°  rotation  of  dipoles  which  previously  were  antiparallel  to  the  field  and  signifies 
the  occur ance  of  a  new  phase. 
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Once  a  self-consistent  solution  for  (16)  was  obtained,  the  magnetization  was  calcu¬ 
lated  by  means  of  (12).  This  magnetization  is  shown,  normalized  to  M«,  the  saturation 
magnetization,  in  Fig.4  as  a  function  of  the  magnetic  field  at  fixed  temperature,  and  in 
Fig.5  as  a  function  of  temperature  at  fixed  magnetic  field.  The  values  plotted  are  those 
for  the  field  in  the  horizontal  direction.  Magnetization  values  with  fields  applied  in  other 
directions  were  also  calculated.  The  values  thus  obtained  differed  from  each  other  by  only 
4 %.  As  seen  in  Fig.5,  the  magnetization  possesses  a  distinct  cusps  at  low  field.  At  higher 
field  the  cusp  is  replaced  by  an  inflection  point. 

The  susceptibility  computed  from  the  magnetization  by  using  (13)  as  a  function  of 
the  applied  magnetic  field  at  fixed  temperatures  is  shown  in  Fig.  6,  and  as  a  function  of 
the  temperature  at  fixed  field  in  Fig.7.  One  again  observes  in  Fig.7  the  occurrence  of  a 
distinct  cusp  for  low  field,  while  at  higher  field  the  cusp  becomes  a  gradual  maximum. 

The  entropy  computed  by  means  of  (14)  is  shown  as  a  function  of  temperature  at  fixed 
field  in  Fig.8,  while  the  specific  heat  computed  from  the  entropy  by  use  of  (15)  is  shown  in 
Fig.O.  Figure  9  again  exhibits  sharp  cusps  in  the  specific  heat  at  the  same  temperatures  as 
the  cusps  in  the  magnetization  and  magnetic  susceptibility.  There  are  no  distinct  features 
in  the  specific  heat  which  can  be  associated  with  the  inflection  points  in  the  magnetization 
shown  in  Fig.5  or  the  gradual  maxima  in  the  susceptibility  shown  in  Fig.7. 

V-Tlie  Phase  Diagram 

On  the  basis  of  the  above  results,  one  can  now  construct  a  phase  diagram  in  the 
r*  -  h*  plana  by  plotting  the  locus  of  the  susceptibility  maxima.  This  is  shown  in  Fig.  10 
for  the  field  applied  along  the  horizontal  direction.  This  diagram  is  also  valid  to  within 
4%  in  h*  for  fields  applied  in  other  directions  in  the  plane. 

The  diagram  shows  three  distinct  phases.  Phase  l  is  the  ordered  phase  where  in- 


teractions  between  the  dipoles  dominate.  Phase  II  is  a  paramagnetic  phase  where  the 
susceptibility  obeys  Curie's  law.  The  transition  between  phase  I  and  phase  II  is  marked 
by  the  occurrence  of  distinct  cusps  in  the  magnetization,  magnetic  susceptibility  and  the 
specific  heat.  In  a  field  applied  along  an  axis  of  symmetry,  that  transition  is  also  marked 
by  an  Sipping  of  dipoles  which  are  at  an  angle  of  180°  with  the  field  in  phase  I,  to  a  parallel 
orientation  with  the  field  in  phase  II.  This  is  seen  in  the  change  in  configuration  between 
frames  b  and  c  in  figures  2  and  3.  In  these  figures  frame  b  is  in  the  ordered  phase,  while 
frame  c  is  in  the  paramagnetic  phase.  Phase  III  represents  a  magnetic  field  induced  ferro¬ 
magnetic  phase.  In  this  phase*  the  dipoles  are  completely  lined  up  with  the  field.  Starting 
in  phase  III,  as  one  increases  the  temperature  at  a  fixed  field,  there  is  an  abrupt  change 
in  the  alignment  of  the  dipoles  with  the  field.  This  change  is  at  the  transition  to  phase  II. 
In  phase  II  the  average  angle  a  dipole  makes  with  the  applied  field  is  greater  than  1°  near 
the  transition.  The  transition  between  phase  III  and  phase  II  is  marked  by  the  inflection 
points  in  the  magnetisation  of  Fig.5  and  the  maxima  in  the  susceptibility  of  Fig.7. 

Vi-Summary 

A  system  of  dipoles  which  are  confined  to  a  plane  exhibits  many  properties  of  a 
physical  system  with  long  range  interactions  and  an  intrinsic  anisotropy.  For  any  particular 
configuration  at  zero  field  and  temperature,  the  energy  of  the  system  can  be  computed 
exactly  for  extended  neighbor  interactions.  Here  we  have  computed  the  ground-state 
energy  of  such  an  extended  system  on  a  honeycomb  lattice  and  found  its  value  to  be 
within  1%  ot  that  determined  for  nearest  neighbor  interactions  only.  Moreover,  we  found 
that  the  ground-etate  degeneracy  specified  by  (21)  and  (23),  which  was  derived  for  only 
nearest  neighbor  interactions,  also  holds  for  extended  neighbor  interactions. 

We  have  also  studied  the  behavior  of  the  system  at  finite  temperature  and  non-zero 
field  applied  in  an  arbitrary  direction  in  the  plane,  under  the  mean-field  approximation.  In 
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this  approximation  only  nearest  neighbor  interactions  were  taken  into  account.  We  found 
that  in  zero  field  the  degeneracy  persists.  In  non-zero  field,  the  magnetic  field  selects  a 
particular  dipole  configuration  which  in  low  fields  resembles  one  of  the  patterns  described 
by  (21)  and  (23).  Two  such  configurations,  for  fields  along  the  horizontal  and  vertical 
directions,  are  shown  in  figures  2  and  3.  A  set  of  states  of  minimum  free-energy  was  then 
selected  from  which  measurable  quantities,  such  as  the  magnetization,  susceptibility  and 
the  specific  heat  were  computed.  These  quantities  show  distinct  maxima  which  denote 
phase  transitions. 

The  phase  diagram  is  shown  in  Fig. 10.  The  nature  of  the  phases  was  established  by 
the  examination  of  the  dipole  configuration  in  each  phase,  the  values  of  the  mean  field 
parameters,  and  the  values  of  the  free  energy  and  the  measurable  quantities  mentioned 
above.  The  phases  are  shown  in  Fig.10.  Phase  I  is  the  low  temperature  low  field  phase 
where  the  dipole  pattern  resembles  one  of  the  ground  states.  In  this  phase  the  interaction 
between  the  dipoles  predominates.  If  the  field  is  applied  along  an  axis  of  symmetry,  one  set 
of  dipoles  will  line  up  at  180°  to  the  field.  The  transition  from  phase  I  to  phase  II  is  marked 
by  the  flipping  of  the  antialigned  dipoles  by  180°  so  that  they  are  now  parallel  to  the  field. 
Phase  II  is  the  high  temperature  phase  and  is  identified  as  the  paramagnetic  phase.  Phase 
III,  the  high  field  phase,  is  the  field  induced  ferromagnetic  phase.  In  this  phase  the  dipoles 
are  completely  lined  up  with  the  field,  the  dipole-magnetic  field  interaction  dominates,  and 
the  mean  field  adds  to  the  magnetization  of  the  system. 
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Table  I 

Pairing  of  Dipoles  as  a  Function  of  the 
Direction  of  Field. 


Field  direction  refers  to  the  angle  the  field  makes  with  the  horizontal  as  referred 


Figure  Captions 


Figure  1.  The  honeycomb  lattice  with  the  subl&ttice  points  numbered  on  one  of  the 
elementary  hexagons. 

Figure  2.  Dipolar  configuration  for  a  field  h *  =  0.84  applied  in  the  horizontal  direction 
with  a-  T *  —  0.68,  b-  T*  =  1.82  and  c-  T*  =  2.00.  a-  and  b-  are  in  the  ordered  phase, 
phase  I.  c-  is  in  the  paramagnetic  phase,  phase  II,  discussed  in  the  text. 

Figure  3.  Dipolar  configuration  for  a  field  h *  =  0.84  applied  in  the  vertical  direction  with 
a-  T*  =  0.68,  b-  T*  =  1.82  and  c-  T*  =  2.00.  a-  and  b-  are  in  the  ordered  phase,  phase  I. 
c-  is  in  the  paramagnetic  phase,  phase  II,  discussed  in  the  text. 

Figure  4.  Magnetization  as  a  function  of  the  applied  external  magnetic  field,  h*  at  fixed 
temperature.  The  magnetization  is  normalized  to  Af,,  the  saturation  magnetization.The 
magnetization  curves  in  increasing  order  of  magnetization  referred  to  the  right  hand  side 
of  the  graph  are  at  T*=3.41,  2.05,  1.50,  and  0.68  respectively. 

Figure  5.  Magnetization  as  a  function  of  temperature  at  fixed  magnetic  fields.  The  mag¬ 
netization  is  given  in  units  of  A/,,  the  saturation  magnetization.  The  magnetization  curves 
in  increasing  order  of  magnetization  referred  to  the  right  hand  side  of  the  graph  are  at 
h*  =0.084,  0.840,  1.680,  1.932,  2.352,  and  2.520  respectively. 

Figure  6.  The  magnetic  susceptibility  for  the  temperatures  of  Fig.4  as  a  function  of  an 
externally  applied  magnetic  field.  The  three  cusps  which  occur  at  increasing  field  are  at 
T*=2.05,  1.50  and  0.68.  The  line  is  at  T*=3.41. 

Figure  7.  The  magnetic  susceptibility  as  a  function  of  temperature  at  fixed  magnetic  field. 
In  order  of  increasing  susceptibility  as  referred  to  the  right  hand  side  of  the  figure  the 
curves  are  at  h*  =2.520,  2.352,  1.032,  1.680,  0.840  and  0.084. 


Figure  8.  The  entropy  at  constant  applied  magnetic  field  as  a  function  of  temperature.  In 
order  of  decreasing  entropy  as  referred  to  the  right  hand  side  of  the  graph  the  values  are 
at  h*=8.41x  10-5,  0.841,  and  1.68,  respectively. 

Figure  9.  The  specific  heat  as  a  function  of  temperature  at  constant  applied  magnetic 
field.  The  fields  are  those  of  Fig. 8  with  the  addition  of  h*  =  1.26.  Above  T*=3  the  specific 
heat  is  a  monotonically  increasing  function  of  the  field. 

Figure  10.  The  phase  diagram  of  the  system  in  the  T*-h*  plane.  Phase  I  is  the  ordered 
phase,  phase  II  is  the  paramagnetic  phase,  and  phase  III  is  a  field  induced  ferromagnetic 
phase. 
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Abstract 

A  theoretical  investigation  of  a  planar  dipolar  system  on  a  honeycomb  lattice  re¬ 
veals  that  the  ground  state  energy  is  increased  for  finite  systems  over  that  of  an  infinite 
system.  The  infinite  system  has  a  ground  state  degeneracy  with  respect  to  dipole  orien¬ 
tations.  This  degeneracy  is  lifted  for  finite  systems.  A  mean  field  calculation  including 
nearest,  second  and  third  neighbors  shows  that  there  are  three  phases,  an  ordered  phase, 
a  quasiparamagnctic  one  which  goes  over  to  a  paramagnetic  one  at  high  temperature 
and  a  magnetic  field  induced  ferromagnetic  phase.  A  defect  in  the  second  nearest  neigh¬ 
bor  has  the  effect  of  introducing  a  finite  magnetisation  at  zero  field  at  low  temperatures. 

Introduction 

The  attraction  of  a  dipolar  system  lies  in  the  fact  that  the  interaction  between  any 
two  dipoles  can  be  easily  and  exactly  calculated.  Dipolar  forces  are  long  range,  and 
there  is  an  inherent  anisotropy  built  into  the  interaction  because  it  depends  not  only  on 
the  relative  orientation  of  the  dipoles  but  also  on  their  location.  Much  work  was  done  in 
the  1930's  on  the  mathematics  of  three  dimensional  dipole  lattices  [1]  and  later  dipolar 
interactions  were  used  to  explain  the  magnetic  transition  of  Cerium  Magnesium  nitrate, 
a  paramagnetic  salt  which  undergoes  an  antiferromagnetic  transition  in  the  millidegree 
range  [2,3].  The  shape  correction  due  to  dipolar  interactions  was  also  explored  [4-6].  In 
this  paper  we  will  confine  ourselves  to  the  two  dimensional  system  of  dipoles.  By  that 
we  mean  that  the  dipoles  are  located  in  a  plane  and  their  rotation  is  such  that  they  al¬ 
ways  point  parallel  to  the  plane  and  rotate  about  an  axis  normal  to  it.  The  magnetic 
field,  however,  is  three  dimensional. 

One  can  show  that  in  a  perfect  dipolar  honeycomb  lattice,  of  infinite  extent  and  at 
zero  temperature,  the  ground  state  is  infinitely  degenerate  [7].  Fig.  1  shows  one  of  a  de¬ 
generate  configurations  from  which  all  the  others  can  be  derived.  One  can  decompose 
the  lattice  into  elementary  hexagons,  with  the  six  dipoles  numbered  1  ,1’,2  ,2’, 3  and  3’ 
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constituting  one  of  those.  Those  with  ”N”  denote  the  next,  second  and  third  nearest 
neighbor  to  the  initial  hexagon.  Dipoles  with  the  same  number  have  the  same  orienta¬ 
tion.  From  this  configuration  one  obtains  all  the  other  states  with  identical  energies  by 
rotating  the  unprimed  dipoles  by  the  same  arbitrary  angle  clockwise  while  the  primed 
dipoles  are  rotated  by  the  same  angle  in  the  counterclockwise  direction,  or  vice  versa. 

It  is  thus  of  interest  to  explore  if  defects  or  boundaries  favor  one  of  the  ground  states 
over  another  8  .  In  the  same  paper  jf.  a  mean  field  calculation  was  performed  which 
extended  the  behavior  of  the  system  to  finite  temperatures  and  finite  magnetic  fields. 
That  calculation,  which  was  confined  to  the  first  nearest  neighbors  only,  showed  that  the 
system  has  a  critical  point  at  zero  field  where  it  undergoes  a  transition  from  an  ordered 
to  a  disordered  state,  and  that  transition  can  be  traced  as  a  function  of  an  applied  mag¬ 
netic  field.  The  change  of  phase  is  marked  by  susceptibility  and  specific  heat  maxima. 

Section  2  of  this  paper  describes  briefly  the  calculation  method.  In  section  3  we  ex¬ 
tend  our  calculations  of  the  ground  state  to  finite  lattices  of  two  shapes  and  examine  the 
total  energy  as  a  function  of  dipole  orientation.  In  section  4  the  mean  field  calculation  is 
extended  to  the  second  and  third  nearest  neighbors.  This  calculation  is  also  performed 
for  the  lattice  having  a  vacancy  in  a  second  or  third  neighbor.  Section  5  presents  a  dis¬ 
cussion  of  the  results. 


Method  of  Calculation 

The  ground  state  energy  calculations  of  an  infinite  and  tne  finite  systems  was  per¬ 
formed  by  summing  the  interactions  of  each  dipole  with  all  the  others  in  the  lattice,  in 
these  calculations  we  used  the  usual  expression  for  the  dipolar  energy  which  is  given  in 
Eq.(l) 

Uif  a  &  J1’  -  3(ft  '  ^iMft  '  nj 

In  this  expression  i  and  j  denote  the  two  interacting  dipoles,  n  is  the  magnetic  mo¬ 
ment  and  r  is  the  distance  between  the  two  dipoles.  If  the  dipoles  are  identical,  one  can 
rewrite  (1)  in  terms  of  unit  vectors  as  shown  in  Eq.(2). 

U%i  *  J[i*  •  ft  ~  3(ft  '  M(ft  *  *u)l  (2) 

From  now  on  we  use  dimensionless  units  where  J»l,  the  temperature  is  given  in  units 
of  J,  the  distance  between  two  nearest  neighbors,  (NO),  will  be  taken  as  one  and  the 
magnetic  field  will  be  in  units  of  (tiH/J). 


An  infinite  lattice  was  simulated  by  taking  into  account  the  interaction  of  3523 
dipoles  on  a  rectangle  of  72.8  XD  by  62-8  XD  and  taking  the  energy  of  only  the  central 
72  dipoles.  Extending  the  boundary  further  had  no  effect  on  the  central  dipoles  thus 
providing  a  natural  cutoff.  The  ground  state  energy  per  dipole  was  thus  computed  as 
E  = -4.450  . 

The  mean  field  calculation  was  performed  by  selfconsistently  solving  the  equations 
of  Lee  et  al  9j  modified  for  dipolar  interactions.  In  that  method,  the  free  energy,  given 
in  Eq.3  is  minimized. 


F(p)  =  kBTTr  p(K  £  UXJ  -  h‘  £  m, 


In  this  equation  F  denotes  the  free  energy,  p  is  the  density  matrix,  K=J/kflT  and 
h*=Hp/l(0T.  This  minimization  gives  an  expression  for  the  magnetization  of  a  dipole 
given  by  Eq.(4). 

m,  =  m(«,)d,  (4) 


where 


/ 1(« 

*>(«) 


(5) 


and  Li  (a)  is  the  modified  Bessel  function  of  order  n.  The  vectors  a,  are  mean  field  pa¬ 
rameters  which  are  determined  from  the  self  consistent  equations 


where  the  prime  over  the  summation  sign  denotes  the  restriction  that  i#j.  An  iterative 
solution  of  £q.(6)  provides  self  consistent  values  of  a,  from  which  F,  the  free  energy,  M, 
the  magnetisation  and  S,  the  entropy,  can  be  obtained  for  any  temperature  and  applied 
magnetic  field.  The  msgnstir  susceptibility  and  the  specific  heat  can  then  be  obtained 


The  Ground  flats 

The  calculation  of  the  energy  per  dipole  was  carried  out  as  described  in  the  previous 
section.  For  an  infinite  lattice,  the  energy  per  dipole  is  E=-4.450.  As  pointed  out  in  the 
introduction,  that  value  is  invariant  under  the  rotation  of  the  unprimed  dipoles  by  an 
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arbitrary  angle  in  the  counterclockwise  direction  with  a  corresponding  rotation  ot  • 
primed  dipoles  in  the  clockwise  direction  or  vice  versa. 

A  finite  array  raises  the  energy  of  the  system.  Moreover,  that  energy  is.  to  a  small 
degree,  dependent  on  the  orientation  of  the  dipoles.  As  an  example  we  calculated  the 
energies  of  two  such  arrays.  The  first  was  170  dipoles  arranged  on  a  square  with  a  side 
length  of  15.5  \*D.  the  other  was  an  array  of  180  dipoles  contained  within  a  circle  of  17 
ND.  In  both  of  these  arrays  2095  of  the  dipoles  were  on  the  boundary,  where  the  bound¬ 
ary  dipoles  were  defined  as  those  which  were  missing  at  least  one  nearest  neighbor.  The 
average  energy  per  dipole,  E,  is  shown  in  Fig. 2  for  both  arrays  as  a  function  of  spin  ori¬ 
entation.  In  this  simulation  the  ground  state  was  preserved  by  rotating  the  dipoles  as 
previously  described,  such  rotation  also  applied  to  the  dipoles  on  the  boundary.  Since 
the  dipoles  numbered  "2"  in  Fig.l  start  with  a  zero  angle  with  respect  to  the  horizontal 
axis,  we  took  that  dipole  as  an  indicator,  and  thus  the  angle,  in  degrees,  denotes  the  '2'' 
dipole  orientation. 

For  the  square  the  average  minimum  energy  per  dipole  was  -4.160  while  the  maxi¬ 
mum  was  -4.072,  for  the  circle  the  values  were  -4.151  and  -3.078  respectively.  The  overall 
change  in  energy  due  to  orientation  is  about  2.595  for  the  square  and  4.3%  for  the  cir¬ 
cle.  Both  have  the  expected  180°  symmetry.  For  the  finite  size  of  the  array  the  average 
energy  was  raised  from  an  average  of  -4.450  to  -4.120  for  the  square  and  -4.064  for  the 
circle.  Although  there  is  an  overlap  in  the  energies  as  shown  in  Fig.2  it  is  somewhat  sur¬ 
prising  that  the  square  has  a  lower  energy  than  the  circle. 

Mean  Field  Calculation 

In  reference  7  we  calculated  the  behavior  of  the  system  taking  into  account  only  the 
nearest  neighbor  interactions.  Thus  the  interactions  of  the  central  hexagon  containing 
the  dipoles  numbered  1, 1’,  2,  2’,  3  and  3’  and  the  nearest  neighbors  of  each  of  them 
were  taken  into  account,  la  this  approximation,  a  primed  dipole,  as  shown  in  Fig.l.  in¬ 
teracted  only  with  an  unpruned  one,  and  vice  versa.  Since  in  a  honeycomb  lattice,  with 
nearest  neighbor  distance  1  ND  the  second  nearest  neighbors  are  v/3  ND  away  while 
the  third  neighbors  are  2  ND  away,  they  contribute  significantly  to  the  interaction.  The 
calculation  was  thus  extended  to  also  include  second  and  third  nearest  neighbors.  The 
magnetization  calculated  in  this  maimer  is  shown  in  Fig.3  for  a  temperature  of  2  and 
as  a  function  of  the  applied  magnetic  field.  This  figure  shows  our  present  calculation, 
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squares,  as  well  as  our  previous  results  which  contained  only  the  nearest  neighbors. . 

for  comparison. 

One  notes  that  the  nearest  neighbor  calculations  (nn)  give  a  lower  magnetization 
than  those  including  second  and  third  nearest  ones  (nnn).  Fig. 4  shows  the  magnetic 
susceptibility  as  a  functon  of  the  applied  field  which  was  obtained  by  differentiation  of 
the  data  shown  in  the  previous  figure.  As  the  field  increases,  one  sees  a  sharp  and  large 
maximum  in  the  susceptibility,  here  seen  at  h=0.7  and  then  a  less  pronounced  peak  near 
h=l.6.  Those  peaks  indicate  transitions  from  phase  I  to  phase  II,  low  field  high  magni¬ 
tude  maximum,  and  from  phase  II  to  phase  III,  high  field  low  magnitude  peak. 

If  one  maps  those  phases  on  the  temperature-magnetic  field  plane  one  obtains  Fig  5. 
Here  the  boxes  denote  the  phase  I  to  II  transition,  while  the  stars  denote  the  phase  II 
to  III  transition,  as  calculated  for  nearest  neighbors  only.  The  "x"  and  the  respec¬ 
tively,  denote  the  same  transition  for  the  nnn  calculation.  One  notes  that  at  low  fields 
the  nnn  results  raise  the  I  to  II  transition  to  higher  temperatures,  while  lowering  it  for 
fields  above  0.75  with  respect  to  the  nn  results.  The  II  to  m  transition  comes  at  a  lower 
field  for  nnn  at  any  particular  temperature  is  long  as  the  temperature  is  above  0.5. 

Phase  I  is  the  ordered  phase  where  dipolar  interactions  dominate.  Phase  II  has  the 
character  of  a  paramagnetic  phase,  however,  one  can  still  detect  the  results  of  dipolar 
interactions.  Because  of  that  we  will  call  it  the  quasiparamagnetic  phase.  At  zero  exter¬ 
nal  magnetic  field  there  are  no  solutions  of  Eq.(6)  in  phase  II  other  than  a,  =  0,  how¬ 
ever,  if  a  field  is  applied,  the  dipoles  align  making  an  angle  with  the  field  between  3°  and 
10°  with  alternate  dipoles  having  a  positive  and  negative  angle.  Those  angles,  of  course, 
become  smaller  as  the  temperature  and  field  are  increased.  In  phase  III  the  dipoles  are 
fully  alligned  with  the  magnetic  field.  This  may  be  a  field  induced  ferromagnetic  phase. 

The  nnn  calculation  was  repeated  with  one  of  the  second  nearest  neighbors  to  the 
central  hexagon  eliminated.  The  eliminated  dipole  is  indicated  in  Fig.l  by  the  small  box. 
Thus  of  the  18  neighbors  of  the  central  hexagon  one  is  missing.  One  thus  introduces  a 
defect  in  the  lattice;  we  shall  call  this  calculation  nnd.  The  magnetisation  and  suscepti¬ 
bility  given  by  this  calculation  are  shown  in  Fig.3  and  4  respectively  and  are  denoted  by 
"x".  Both  the  overall  magnetization  and  the  susceptibility  of  this  calculation  are  similar 
to  the  nnn  results.  However,  the  magnetization  of  the  nnd  is  generally  greater  in  phase 
I, the  same  in  phase  Q,  and  smaller  in  phase  III  as  compared  to  that  of  nnn.  The  phase 
boundaries  appear  to  be  unaffected  by  the  defect. 
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In  nnd  there  is  a  finite  magnetization  at  zero  and  low  fields.  Once  the  defect  in¬ 
troduced  it  provides  a  focus  for  polarization.  Thus  at  low  fields  the  the  defects  enhance 
the  magnetization.  As  the  field  and  temperature  are  raised,  the  influence  of  the  defect 
becomes  insignificant  and  thus  the  equality  between  the  nnn  and  nnd  results  in  phase 
II.  In  phase  III  the  defect  provides  a  disturbance  in  the  lattice  regularity  and  acts  as  a 
cavity  in  a  ferromagnetic  system,  decreasing  the  overall  magnetization.  Thus  the  lower 
magnetization  in  phase  III.  Figures  6  and  7  show  this  effect  in  greater  detail  for  the  low 
field  region  where  both  the  susceptibility  and  magnetization  are  plotted  for  the  nnn.  box 
and  nnd,  "x",  calculations.  Fig.d  shows  that  in  nnn  the  magnetization  goes  to  zero  in 
zero  field  while  it  is  finite  in  the  nnd  case.  Fig.7,  which  shows  the  susceptibility  with 
0.3  subtracted  from  the  actual  value,  indicates  that  in  nnd  the  susceptibility  has  a  small 
minimum  at  low  fieldj,  while  in  nnn  it  is  a  monotonieally  increasing  function  in  the  low 
field  region. 

Conchialocis 

The  inclusion  of  boundaries  in  a  system  of  dipoles  confined  to  rotate  in  a  plane  and 
located  on  a  honeycomb  lattice  raises  the  average  ground  state  energy.  Moreover,  de¬ 
pending  on  the  shape  of  the  boundaries,  the  average  ground  state  energy  depends  on  the 
dipole  orientation. 

The  inclusion  of  second  and  third  nearest  neighbors  in  the  interaction  alters  the 
phase  boundaries  of  the  system.  Defects  introduce  a  finite  magnetization  at  low  tem¬ 
perature,  even  at  zero  field.  In  low  fields,  the  susceptibility  appears  to  have  a  minimum 
with  the  temperature  held  constant  for  calculations  including  defects,  while  it  is  a  mono- 
tonically  increasing  function  of  the  field  when  no  defects  are  included. 

The  system  studied,  might  be  applicable  to  some  two  dimensional  systems  where 
dipolar  interactions  play  an  important  role.  For  one,  such  systems  migh  be  found  in 
graphite  intercalated  compounds  where  a  magnetic  substance  is  the  in terc slant.  The 
graphite  layers  thus  reduce  the  interaction  between  the  in  terc  slant  layers  creating  a  mag¬ 
netic  two  dimensional  array  within  each  interc slant  layer. 
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Figure  Captions 

Fig.l  The  lattice  with  ground  state  dipole  directions.  >'  denotes  neighbors  of  the  main 
hexagon  taken  into  account  in  the  mean  field  calculation. 

Fig.2  The  energy  per  dipole  of  the  square,”  x”,  and  the  circular,  box.  array  of  dipoles,  as 
a  function  of  the  angle  of  dipole  #2. 

Fig.3  Magnetisation  at  T*2  as  a  function  of  the  applied  field  showing  nn, " ,  nnn,  box. 
and  nnd,  "x",  results. 

Fig.4  The  Magnetic  susceptibility  at  T*2  as  a  function  of  the  applied  field  with  symbols 
as  in  Fig.3. 

Fig.3  The  phase  diagram  of  the  system  with  the  phases  described  in  the  text.  The  boxes 
and  the  denote  the  nn  calculation  while  ”x”  and  ”+"  show  th*  nnn  and  nnd  results. 


Fig.6  The  low  field  behavior  of  the  magnetization  at  T  =  2.  The  box  denotes  the  mm  ra 
culation  while  "x“  denotes  nnd. 


APPENDIX  X 


Self-balancing  raslstanca  bridga 

Alvaro  Kaplan,  QnorgnO.  Zlmmarman,  and  Doan  S.  Edmonds,  Jr. 

Department  of  Phftkt  Batten  Uruventty,  Batten.  Mtaochusttti  022 1 J 
(Received  27  May  1986;  accepted  for  publication  28  July  1986) 

An  automatically  balanced  resistance  bridge  with  bekm-picowatt  dissipation  in  the  measured 
resistance  has  been  developed.  The  bridge,  capable  of  measuring  resistances  between  20  ft  and 
10  left  with  an  accuracy  of  3%,  was  specifically  designed  to  allow  resistance  measurements  to 
be  made  at  low  temperatures  without  introducing  an  appreciable  heat  source  to  the  cryostat  It 
consists  of  a  Wheatstone  bridge  having  a  photoresistor  as  one  of  its  arms.  The  feedback  balance 
is  accomplished  by  means  of  light-emitting  diodes  which  control  the  photoresistor.  A  novel 
feedback  principle  used  in  the  bridge  is  described. 


In  measurements  of  resistance  of  materials  at  temperatures  erable  tradition  as  secondary  thermometers  in  low-tempera- 

below  the  boiling  point  of  He*  (below 4.2  K),  it  is  important  ture  physics,2  and  are  even  now  being  used  in  ultralow-tem- 

that  the  power  dissipated  hi  the  material  whose  resistance  is  perature  applications.  We  have  developed  a  self-contained, 

being  measured  be  small  This  is  necessitated  by  the  fact  that  self-balancing  resistance  bridge  whose  power  dissipation  is 

at  those  low  temperatures,  the  specific  heats  of  substances  less  than  1  pW.  It  is  basically  a  Wheatstone  bridge  with  a 

become  small  and  the  thermal  contact  between  substances  lightsensitive  resistance  as  one  arm.  This  resistance  is  pack- 

end  helium  is  limited  by  the  various  thermal  “resistances,"  aged  with  seven  LED’s  which  shine  on  it  and  are  driven  by 

the  one  most  considered  being  the  Kapitza  resistance.1  Thus  the  control  circuit  output  so  that  at  balance  the  unknown 

if  too  much  power  is  dissipated  in  the  resistive  element  dur-  resistance  Xt  is  given  as  a  function  of  the  control  circuit 

ing  the  measurement  process,  the  temperature  of  the  sub-  output  voltage  V0  by  a  relation  of  the  form 

stance  whose  resistance  is  being  measured  may  bear  little  y  +  BIX  )~U1  ( I ) 

relation  to  the  temperature  of  its  surroundings.  This  point  is  0  "  5 

crucial  if  the  subetiutce  in  question  is  serving  at  a  thermom-  Figure  I  is  a  graph  of  the  actual  results  obtained  with  our 

eter.  setup  and  allows  values  for  the  constants  A  and  B  to  be  deter- 

Carbon  resistors,  carbon  glass  resistors,  phosphor  mined.  Individual  operating  conditions  will  affect  these  val- 

bronze  resistors,  and  o, her  resistive  elements  have  had  a  ven-  ties,  but  the  form  of  Eq.  ( I )  should  hold  in  any  case 
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Fig.  1.  Output  voltage 
plotted  against  sense  resis¬ 
tance  Av  A  linear  plot  is  ob¬ 
tained  by  plotting  K„  against 


l/VRUi  10* 


The  complete  circuit  of  the  self-balancing  bridge  is  given 
in  Fig.  2.  The  ac  excitation  for  the  bridge  is  a  340-Hz  sine 
wave  obtained  from  the  type  8038  oscillator  chip.1  The  fre¬ 
quency  of  this  oscillator  is  given  by 

RCJ~  0.309,  •  (2) 


where  R  —  Ry  =■  R,  is  in  ohms,  C2  is  in  farads,  and  /is  in 
hertz.  Ry  and  R ,  are  made  equal  when  a  symmetrical  square 
wave  is  desired,  and  Eq.  (2)  applies  when  pins  7  and  8  of  the 
8038  are  connected  together.  In  the  circuit  of  Ftg.  2  a  fre¬ 
quency  control  voltage  is  introduced  at  pin  8  to  allow  fine 
tuning  the  8038  to  the  exact  frequency  of  the  tuned  error- 
signal  amplifier.  The  values  shown  give  a  340-Hz  output 
signal  with  R,0  set  near  midrange. 

The  output  signal  from  the  bridge,  when  not  balanced,  is 


amplified  by  a  narrow-band  twin-T  push-pul!  amplifier  built 
for  us  by  the  White  Instrument  Company  under  model  num¬ 
ber  4529.*  The  amplified  340-Hz  sine  wave  is  then  detected 
by  a  full-wave  synchronous  rectifier,  a  phase-sensitive  detec¬ 
tor  using  a  type  4066  quadruple  gate  chip,  whose  output  is 
further  amplified  but  not  filtered,  so  that  we  have  at  this 
point  a  rectified  sine  wave,  either  positive-going  or  negative¬ 
going  depending  on  which  way  the  bridge  is  out  of  balance. 
The  distinguishing  feature  of  our  arrangement  is  that  this 
signal  is  not  used  directly  to  control  the  balancing  of  the 
bridge,  but  is  instead  applied  to  an  operational  amplifier  hav¬ 
ing  no  feedback  and,  therefore,  functioning  as  a  flip-flop 
with  two  stable  states.  When  the  signal  is  positive  going,  the 
flip-flop  spends  most  of  its  time  in  the  state  giving  full  posi¬ 
tive  output,  and  capacitor  C,  charges  through  resistor  R, 
accordingly.  The  resulting  dc  voltage  is  amplified  by  the 
pc 2002  power  amplifier  which  drives  the  LED's  already 
mentioned  and  provides  output  voltage  V0  across  the  cur¬ 
rent-sensing  resistor  R2.  When  the  amplified  error  signal  is 
negative  going,  C,  discharges  through  R,  but,  at  balance, 
flip-flop  A,  flips  back  and  forth  randomly,  spending,  on  the 
average,  as  much  time  in  one  state  as  in  the  other,  so  that  the 
charge  on  C,  remains  constant.  Thus  the  correct  LED  cur¬ 
rent  is  maintained  with  zero  error  signal,  giving  the  system 
the  characteristics  of  a  first-order  servomechanism  rather 
than  a  simple  regulator.  Zener  diode  CR,  and  offset  control 
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Rm  serve  to  asymmetrize  the  system  so  that  the  output  of  the 
/ipc2002  will  not  go  negative.  ( 

The  system  comes  to  balance  in  a  reasonable  time  after  a 
change  occurs  in  R,  and  has  been  used  to  make  temperature 
measurements  to  I  mK  in  the  helium  temperature  range  ( I- 
4  K).  The  system  thus  provides  an  inexpensive,  automated, 
accurate  alternative  to  the  traditional  manually  balanced 
Wheatstone  bridge. 

This  work  was  supported  in  part  by  Air  Force  of  Scien¬ 


tific  Research  under  Grant  No.  AFOSR  82-0286. 
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